SEARCH · Search NASA

Results for “diffusion processes”

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 937 records · Page 52

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Multi-scale model based design of membrane reactor/separator processes for intensified hydrogen production through the water gas shift reaction

This paper aims to first quantify the impact of various diffusion models (Maxwell-Stefan, Wilke, Dusty-Gas) on the predictions of a multi-scale membrane reactor/separator mathematical model, and to then demonstrate this model's use for the design and process intensification of membrane reactor/separator systems for hydrogen production. This multi-scale model captures velocity, temperature and species' concentration profiles along the catalyst pellet's radial direction, and along the reactor's axial direction, by solving the momentum, energy, and species transport equations, accounting for convection, conduction, reaction, and diffusion mechanisms. In the first part of work, the effect of pellet-scale design parameters (mean pore diameter, volumetric porosity, tortuosity factor, etc.) and various species' flux models on the model predictions is studied. In the second part, the study focuses on the comparison, in terms of their process intensification characteristics, of various hydrogen production processes. These include a conventional high-temperature shift reactor (HTSR)/low-temperature shift reactor (LTSR) sequence, a novel HTSR/membrane separator (MS)/LTSR/MS sequence, and a process that involves low-temperature shift membrane reactors-LTSMR in a series.

08 HYDROGEN↗

DOE OSTI2020

Influence of Coating Method on the Performance of Roll-to-Roll Coated PEM Fuel Cell Catalyst Layers

As production of fuel cell vehicles and power systems increases there is a need to understand the process science of fuel cell manufacturing. Industrial production of polymer electrolyte membrane fuel cell (PEMFC) electrodes will rely on continuous roll-to-roll (R2R) coating methods to meet the rates needed for mass production. There are a variety of coating methods that could be used to coat catalyst layers. The different physics of each method may result in different morphologies and performance. Here we present the results of a study comparing fuel cell catalyst layers coated using two R2R coating methods: slot-die and gravure. These two methods were selected because the physics of applying the liquid to the substrate is very different between the two cases. Pt/C catalysts layers were coated directly onto carbon-fiber gas diffusion media to create gas diffusion electrodes (GDEs), without the use of a decal transfer process. These GDEs were assembled into membrane electrode assemblies (MEAs) by hot pressing to Nafion membranes. In situ performance and other electrochemical diagnostics were used to determine the influence of coating method on catalyst layer electrochemical properties. Electron microscopy was utilized to understand the influence of coating method on electrode morphology and its influence on electrochemical properties. We also explored the influence of solvent ratio to understand the potential coupling between ink formulation and coating methodology.

DIRECT ENERGY CONVERSION,MATERIALS SCIENCE↗

DOE OSTI2024

Phase field-volumetric lattice Boltzmann model of ion uptake in porous nuclear waste form materials under continuous flow

The flow field within the mesopores of sorbent particles plays a crucial role in radionuclide diffusion and ion uptake kinetics, thus, impacting the overall performance of porous nuclear waste form materials. To fundamentally understand the influence of microstructures and material properties on the radionuclide absorption and retention processes requires a coupled multi-physics model that considers the advection and diffusion within the flow field, the reaction at liquid-solid interfaces, and finally, the solid-state diffusion within a complex nanoporous medium. Here, this study employs the volumetric lattice Boltzmann method (VLBM) to accurately and efficiently calculate the steady state velocity field inside the mesopores of sorbent particles. The obtained velocity field is then utilized to calculate the advection of ions in the steady flow. A phase field (PF) model of ion uptake is used to describe the reaction occurring at the solid-liquid interface and diffusion inside the porous medium. The integrated PF-VLBM model is verified in terms of the mass conservation and numerical efficiency and validated qualitatively with experimental observation data. Then, it is applied to study the influence of thermodynamic and kinetic properties, as well as flow field conditions on the ion uptake kinetics. The numerical results demonstrate that the ion uptake kinetics in porous particles has three distinct stages, which is in agreement with the observations in continuous flow experiments. In the first stage, the kinetics is predominantly controlled by the flow field and ion diffusivity in the liquid phase. The kinetics in the second stage is primarily governed by ion diffusivity in the solid phase. In the third stage the system reaches a dynamic equilibrium with a net zero uptake flux at the interface. It is also found that porous structures significantly affect the efficiency and capacity of ion uptake. The simulation results can help to understand the physics behind the observed ion uptake kinetics in experiments and to facilitate the development of constitutive equations that can account for heterogeneous microstructures in engineering performance codes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2018

Computational Studies of Solid Stoppers for Generating Intense Exotic Radioactive Ion Beams at FRIB

Low-emittance beams of rare isotopes at low energies (less than 3 MeV/u) are required for many high precision measurements in nuclear science research. Such beams will be produced at the Facility for Rare Isotope Beams (FRIB) by first slowing down and thermalizing or stopping fast fragmentation beams, and subsequently reaccelerating them to the required energies for experiments. The use of a solid stopper has the potential to provide high-intensity reaccelerated beams of many species. We performed simulations of the processes critical for generating beams in solid stoppers, including beam heating, external heating, diffusive release, effusive transport, and surface sticking. Using a variety of software tools, we investigated the viability of numerous solid stopper system designs for FRIB. We present our computational approach, the physical processes we modeled, the codes we used, and our results. We conclude with a proof-of-concept solid stopper system that releases ions sufficiently quickly for experiments with unstable reaccelerated beams at FRIB.

47 OTHER INSTRUMENTATION↗

DOE OSTI2018

Characterization of Grain‐Size Distribution, Thermal Conductivity, and Gas Diffusivity in Variably Saturated Binary Sand Mixtures

Core Ideas Parameterized models are relevant to applications where different porous media mixtures are used. Models include a grain‐size distribution function to describe bimodal behavior for binary mixtures. Improved model describes observed thermal conductivity–saturation relations for binary mixtures. Combined Buckingham–Penman model used to describe observed gas diffusivity–air content relations. This work is relevant for proper simulation of mixed porous media. Characterization of differently textured porous materials, as well as different volumetric porous media mixtures, in relation to mass and heat transport is vital for many engineering and research applications. Functional relations describing physical (e.g., grain‐size distribution, total porosity), thermal, and gas diffusion properties of porous media and mixtures are necessary to optimize the design of porous systems that involve heat and gas transport processes. However, only a limited number of studies provide characterization of soil physical, thermal, and gas diffusion properties and the functional relationships of these properties under varying soil water contents, especially for soil mixtures, complicating optimization efforts. To better understand how mixing controls the physical, thermal, and gas diffusion properties of porous media, a set of laboratory experiments was performed using five volumetric mixtures of coarse‐ and fine‐grained sand particles. For each mixture, the grain‐size distribution (GSD), thermal conductivity, and gas diffusivity were obtained and parameterized using existing and suggested parametric models. Results show that the extended, two‐region Rosin–Rammler particle‐size distribution model proposed in this study could successfully describe the bimodal behavior of the GSD of binary mixtures. Further, the modified Côté and Konrad thermal conductivity model adequately described the thermal conductivity–water saturation relations observed in different mixtures. The proposed simple soil‐gas diffusivity descriptive model parameterized the upper limit, average, and lower limit behavior in gas diffusivity–air content relations in apparently texture‐invariant gas diffusivity data. Results further show a close analogy between gas diffusivity and thermal conductivity and their variation with saturation across different binary mixtures. Overall, the results of the study provide useful numerical insight into the physical, thermal, and gas transport characteristics of binary mixtures, with wide implications for future engineering and research applications that involve multicomponent porous systems.

Deepagoda, T. K. K. Chamindu↗

DOE OSTI2025

Understanding Oxide–Metal Interactions During Hot Isostatic Pressing to Diffusion Bond Aluminum Alloy 6061 Plates

The interaction between Mg, Si, and Al 2 O 3 during hot isostatic pressing diffusion bonding of aluminum alloy 6061 (AA6061) plates was investigated through thermodynamic calculations and experimental microstructural characterization. Thermodynamic calculations as functions of temperature, pressure, and composition revealed that the interaction among Mg, Si, and Al 2 O 3 yields Mg 2 Si and either MgO + Al or MgAl 2 O 4 + Al, facilitating the reduction of Al 2 O 3 and allowing Al/Al metallic bonds to form. Total pressure variation had a negligible influence on the oxygen partial pressure, and consequently, the reaction product formation. Oxygen partial pressure variation as a function of temperature and initial amount of Al 2 O 3 determined the formation of either MgO or MgAl 2 O 4 . Experimental Hot Isostatic Pressure (HIP) bonding at 723 K and 833 K under a constant pressure of 1017 atm documented the cooling rate-dependent formation of β–Mg 2 Si precipitates. High-resolution transmission electron microscopy imaging and selected area electron diffraction patterns verified the formation of β–Mg 2 Si and MgO at the interface but did not detect MgAl 2 O 4 . In conclusion, findings from this study clarify the role of thermochemical interactions in oxide disruption and bonding mechanisms during HIP diffusion bonding of AA6061 and provide guidance for optimizing joining processes for monolithic nuclear fuel assemblies.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Interfacial Void Formation and Self-Healing in Oxide Scales on Al-containing High-Entropy Alloy

The exceptional high-temperature oxidation resistance of Al-containing high-entropy alloys (HEAs) is often attributed to the formation of a protective α-Al 2 O 3 scale. However, the dynamic, atomic-scale mechanisms governing the stability of this scale—including interfacial void formation and the often-postulated but rarely visualized “self-healing” capacity—remain poorly understood. Herein, we reveal the complex evolution of the triple-layer oxide scale on an Al 10 CoCrFeNi HEA through combined electron microscopy and diffraction study. We show that interfacial voids are an inherent consequence of the scaling process, originating from two distinct mechanisms: the Kirkendall effect at the interface between the γ-Al 2 O 3 /α-Al 2 O 3 and alloy driven by cationic diffusion imbalance and volumetric contraction due to phase transformations at the spinel/Cr 2 O 3 interface. Crucially, we provide microstructural evidence consistent with an intrinsic self-healing response. This process is driven by coupled inward diffusion of oxygen and outward diffusion of metal cations, leading to the in-situ formation of transient θ-Al 2 O 3 and spinel phases that partially fill and seal the voids. Here, these results provide atomic-scale insights into the phase evolution, defect formation, and self-repair of oxide scales in HEAs—highlighting pathways to enhance their oxidation resistance in extreme environments.

High-entropy alloy↗

DOE OSTI2024

Dynamics of High Molecular Weight Cylindrical and Lamellar Block Copolymers with X‐ray Photon Correlation Spectroscopy

Here, the structure and dynamics of polystyrene (PS)-b-poly(ethylene oxide) block copolymers (BCPs) are studied. The BCPs exhibit microphase-separated cylindrical and lamellar morphologies. Structural dynamics are measured with X-ray photon correlation spectroscopy in the small-angle regime. Morphologies and domain sizes are evaluated using small-angle X-ray scattering (SAXS), scanning electron microscopy, and atomic force microscopy. Different solvent processing conditions are investigated. Grain sizes evaluated using SAXS are found to depend on processing only for the rubbery majority BCP. The structural relaxation times are examined as a function of PS volume fraction, temperature, morphology, and structural sizes. Well above the glass transition temperature (T g ) of PS, all samples exhibit stretched autocorrelation decays and diffusive dynamics. Near T g of PS, the dynamics of all samples are anomalous with compressed autocorrelation decays and hyperdiffusive dynamics. This transition occurs at 153 °C or 1.13 T g of PS. In the diffusive regime (at high temperature), structural relaxation times are dependent on the processing method. Near PS T g (at low temperature), structural relaxation times scale with the PS volume fraction. Structural relaxation times do not correlate with grain size, indicating that the out-of-equilibrium state of PS dominates the structural dynamics of these strongly phase-segregated BCPs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Nature of Moisture-Induced fogging defects in scintillator plastic

Polyvinyl toluene (PVT) scintillator plastic may degrade in the field due to inward water diffusion at elevated temperatures, which exceeds the saturation limit at lower temperatures, and after long periods of time (e.g., many years), leads to the formation of disk-like defects that attenuate scintillation light, leading to detector degradation. In this work, using fractography and high-magnification optical and electron microscopy to characterize the water-induced defects, a model of the fogging process is hypothesized as follows: excess water present at low temperatures diffuses to spheroids to minimize contact with the hydrophobic polymer. Through the hydrophobic effect, the entropy of water increases by forming nanoclusters which minimizes the contact between the water and the PVT. As the water nanoclusters grow, they break and fold the polymer into densely-packed crystalline regions creating more space for water within the spheroids. The polymer outside the spheroid resists the shrinkage which builds up tension within the spheroid. Once the tensile stress exceeds the yield strength of the plastic, the spheroid is torn in half resulting in a defect. Excess water then drains into the cavities along the disk thereby further increasing its entropy. Slower cooling (over 1 day) leads to larger spheroids and hence, larger “permanent” defects. Freezing causes some defects to further grow due to the expansion of water to ice. These findings imply that the remaining lifetime of scintillator plastic in the field could be predicted using temperature and humidity data thereby mitigating security risks of degraded radiation portal monitors.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Bio‐Imaging quorum sensing signal molecules in a soil‐mimic gel

Rhizosphere, the narrow region surrounding plant roots directly influenced by root exudates and the root associated microbiome, plays an important role in plant productivity and the rhizosphere is well known for stimulating microbial metabolic activities. How microbial communities interact to form stable, metabolically interconnected functional communities is an area of intense interest. The question remains on how microbially produced secreted molecules that function as intercellular communication signals shape the structure and function of microbial communities. Our goal is to detect and quantify signal molecules produced by microbes or plants in the root‐soil environment, spatially and temporarily. As a proof‐of‐concept, we are imaging diffusible extracellular microbial metabolites involved in a bacteria‐bacteria communication process called quorum‐sensing. Quorum‐sensing relies on the accumulation of high concentrations of signal molecules in the environment to control bacterial gene expression, influencing rhizosphere colonization and plant health. Local concentration of quorum sensing molecules are determined using aptamer based sensors. Aptamers that can specifically bind the desired signal molecules are selected through SELEX and immobilized on the surface of nano‐porous membranes. Binding of the signal molecules and aptamer covered surface results in changes in surface charge distribution and steric hinderance and thus modifying the transmembrane ionic transport. Changes in transmembrane impedance can be measured through electrochemical impedance spectroscopy methods to monitor local concentrations of signal molecules. Sensor responses were determined for different concentrations of signal molecules and results showed detection of C4‐HSL in a soil‐mimic solution with K D of 10 nM. We demonstrate that the quorum‐sensing signal molecule C4‐homoserine lactone and (C4‐HSL) can rapidly diffuse in a soil‐mimic gel, providing evidence for our use of a soil‐mimic for further aptasensor development and validation. The sensors were then inserted into soil mimic gel for monitoring C4‐HSL diffusion and the resulted impedance changes were used to determine the C4‐HSL concentration at different positions in the gel. Measurements of local C4‐HSL concentrations variation and numerical solution of diffusion equation were used to create 4D images of C4‐HSL molecule diffusion in the soil‐mimic gel.

Jiang, Nianyu↗

DOE OSTI2023

Impact of oxidation and carbonation on the release rates of iodine, selenium, technetium, and nitrogen from a cementitious waste form

Evaluation of the long-term retention mechanisms and potential release rates for the primary constituents of potential concern (COPCs) (i.e., Tc, I, Se, and nitrate) is necessary to determine if Cast Stone, a radioactive waste form, can meet performance objectives under near-surface disposal scenarios. In this work, a mineral and parameter set accounting for the solubility of I and Se in Cast Stone was developed based on pH-dependent and monolithic diffusion leaching test results, to extend a geochemical speciation model previously developed. The impact of oxidation and carbonation as environmental aging processes on the retention properties of Cast Stone for primary COPCs was systematically estimated. Physically, the effective diffusion coefficients of 4 COPCs in Cast Stone were increased after carbonation and/or oxidation, reflecting an increase in permeability to diffusion. Chemically, i) pH & pe conditions in the original Cast Stone were favorable for the stabilization of Tc, but not for I, Se, and N; ii) oxidation (with/without carbonation) of Cast Stone changed the pe & pH conditions to be detrimental for Tc stabilization; and iii) carbonation (with/without oxidation) of Cast Stone modified the pH & pe conditions to be beneficial for the stabilization of I (in system with Ag added) and Se.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗

DOE OSTI2017

Matching time and spatial scales of rapid solidification: Dynamic TEM experiments coupled to CALPHAD-informed phase-field simulations

A combination of dynamic transmission electron microscopy (DTEM) experiments and CALPHAD-informed phase-field simulations was used to study rapid solidification in Cu–Ni thin-film alloys. Experiments—conducted in the DTEM—consisted of in situ laser melting and determination of the solidification kinetics by monitoring the solid–liquid interface and the overall microstructure evolution (time-resolved measurements) during the solidification process. Modelling of the Cu–Ni alloy microstructure evolution was based on a phase-field model that included realistic Gibbs energies and diffusion coefficients from the CALPHAD framework (thermodynamic and mobility databases). DTEM and post mortem experiments highlighted the formation of microsegregation-free columnar grains with interface velocities varying from ~0.1 to ~0.6 m s –1 . After an 'incubation' time, the velocity of the planar solid–liquid interface accelerated until solidification was complete. In addition, a decrease of the temperature gradient induced a decrease in the interface velocity. The modelling strategy permitted the simulation (in 1D and 2D) of the solidification process from the initially diffusion-controlled to the nearly partitionless regimes. Lastly, results of DTEM experiments and phase-field simulations (grain morphology, solute distribution, and solid–liquid interface velocity) were consistent at similar time (μs) and spatial scales (μm).

36 MATERIALS SCIENCE↗

DOE OSTI2017

Mathematical simulation of heat and mass transfer in a cylindrical channel versus the tangential momentum accommodation coefficient

The process of heat and mass transfer in a long cylindrical channel has been considered in terms of the mirror–diffuse model of the Maxwell boundary condition. The Hazen–Williams equation is used as a basic equation of the process kinetics. A constant temperature gradient is maintained in the channel. The heat and mass fluxes through the cross section of the channel versus the tangential momentum accommodation coefficient have been calculated in a wide range of the Knudsen number. The heat flux profiles have been constructed. A comparison with relevant published data has been carried out.

42 ENGINEERING↗

DOE OSTI2020

Tunable chemical complexity to control atomic diffusion in alloys

Abstract In this paper we report a new fundamental understanding of chemically-biased diffusion in Ni–Fe random alloys that is tuned/controlled by the intrinsic quantifiable chemical complexity. Development of radiation-tolerant alloys has been a long-standing challenge. Here we show how intrinsic chemical complexity can be utilized to guide the atomic diffusion and suppress radiation damage. The influence of chemical complexity is shown by the example of interstitial atom (IA) diffusion that is the most important defect in radiation effects. We use μs-scale molecular dynamics to reveal sluggish diffusion and percolation of IAs in concentrated Ni–Fe alloys. We develop a mean field diffusion model to take into account the effect of migrating defect energy properties on diffusion percolation, which is verified by a new kinetic Monte Carlo approach addressing detailed processes. We demonstrate that the local variations in the ground state energy of IA configurations in alloys, reflecting the chemical difference between alloying components, drives the percolation effects for atomic diffusion. Percolation, chemically-biased and sluggish diffusion are phenomena that are directly related to the chemical complexity intrinsically to multicomponent alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Electron Injection via Modified Diffusive Shock Acceleration in High-Mach-number Collisionless Shocks

The ability of collisionless shocks to efficiently accelerate nonthermal electrons via diffusive shock acceleration (DSA) is thought to require an injection mechanism capable of preaccelerating electrons to high enough energy where they can start crossing the shock front potential. We propose, and show via fully kinetic plasma simulations, that in high-Mach-number shocks electrons can be effectively injected by scattering in kinetic-scale magnetic turbulence produced near the shock transition by the ion Weibel, or current filamentation, instability. We describe this process as a modified DSA mechanism where initially thermal electrons experience the flow velocity gradient in the shock transition and are accelerated via a first-order Fermi process as they scatter back and forth. The electron energization rate, diffusion coefficient, and acceleration time obtained in the model are consistent with particle-in-cell simulations and with the results of recent laboratory experiments where nonthermal electron acceleration was observed. This injection model represents a natural extension of DSA and could account for electron injection in high-Mach-number astrophysical shocks, such as those associated with young supernova remnants and accretion shocks in galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Formation, Structural Variety, and Impact of Antiphase Boundaries on Li Diffusion in LiCoO 2 Thin-Film Cathodes

Anti-phase boundaries (APBs) are observed in as synthesized and processed LiCoO 2 , which is used widely as a cathode material. Using a combination of scanning transmission electron microscopy and ab initio simulations, we investigate APB formation, structure, stability, and impact on Li ion diffusion. It is found that APB low-coordinated sites give rise to cation defects which, in turn, provide a variety of APB structures. Li diffusion along an APB can enter a correlated mode, leading to a ~40% decrease in activation energy with respect to that for the uncorrelated hopping of Li ions. We propose that APBs function as additional mass transfer channels that couple in-plane Li ion diffusion pathways, thus, facilitating Li transfer from one two-dimensional basin to another, potentially enabling new energy storage architectures.

Ong, Phuong-Vu↗

DOE OSTI2019

Dynamical properties of a room temperature ionic liquid: Using molecular dynamics simulations to implement a dynamic ion cage model [plus supplemental information]

The transport behavior of ionic liquids (ILs) is pivotal for a variety of applications, especially when ILs are used as electrolytes. Many aspects of the transport dynamics of ILs remain to be understood. Here, a common ionic liquid, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (BmimNTf 2 ), was studied with molecular dynamics simulations. The results show that BmimNTf 2 displays typical structural relaxation, subdiffusive behavior, and a breakdown of the Stokes-Einstein diffusion relation as in glass-forming liquids. In addition, the simulations show that the translational dynamics, reorientation dynamics, and structural relaxation dynamics are well described by the Vogel-Fulcher-Tammann equation like fragile glass forming liquids. Building on previous work that employed ion cage models, it was found that the diffusion dynamics of the cations and anions were well described by a hopping process random walk where the step time is the ion cage lifetime obtained from the cage correlation function. Detailed analysis of the ion cage structures indicated that the electrostatic potential energy of the ion cage dominates the diffusion dynamics of the caged ion. The ion orientational relaxation dynamics showed that ion reorientation is a necessary step for ion cage restructuring. The dynamic ion cage model description of ion diffusion presented here may have implications for designing ILs to control their transport behavior.

36 MATERIALS SCIENCE↗