Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ionic diffusivity”

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.

431 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The “good,” the “bad,” and the “hidden” in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl 2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca 2+ –Cl - interactions were primarily “hidden” in the Ca-NDIS signal due to overlapping with Ca 2+ –O w and Ca 2+ –H w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models “on average,” the question “which model is the best” has not been answered but replaced by the question “which model is better for a given research.” An overall “good” model can still be inappropriate or, in some instances, “bad” and, unfortunately, produce erroneous results. Additionally, the accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

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

NMR Relaxometry and Diffusometry Analysis of Dynamics in Ionic Liquids and Ionogels for Use in Lithium-Ion Batteries

Solvation and transport properties of methlyviologen dichloride in mixtures of choline chloride (ChCl) and ethylene glycol (EG) including the deep eutectic solvent (DES) ethaline were studied through physical property measurements, cyclic voltammetry, and Raman spectroscopy. The transport properties of the solvate structures in ethaline and ChCl:EG mixtures were compared to a well-studied ionic liquid (IL), 1-methyl-1-propylpyrrolidinium bis(trifluoromethylsulfonyl)imide, [PYR13][TFSI] through the lens of hole theory. Exchanging the chloride anion of ethaline with bis(trifuoromethylsufonyl)imide) anion, [TFSI], is found to increase the fluidity of the solvent and promoted stronger solute-solvent interactions. Raman spectroscopy shows that methylviologen is solvated by EG in solutions with ChTFSI but interstitially accommodated in void spaces in both ChCl solutions and [PYR13][TFSI]. Complex solvents such as DESs are regarded as “designer solvents” and we demonstrate that the physical properties and solvation characteristics of these fluids can be varied by anion exchange of hydrogen bond acceptors.

25 ENERGY STORAGE↗
DOE OSTI2020

From Salt in Solution to Solely Ions: Solvation of Methyl Viologen in Deep Eutectic Solvents and Ionic Liquids

Solvation and transport properties of methlyviologen dichloride in mixtures of choline chloride (ChCl) and ethylene glycol (EG) including the deep eutectic solvent (DES) ethaline were studied through physical property measurements, cyclic voltammetry, and Raman spectroscopy. The transport properties of the solvate structures in ethaline and ChCl:EG mixtures were compared to a well-studied ionic liquid (IL), 1-methyl-1-propylpyrrolidinium bis(trifluoromethylsulfonyl)imide, [PYR 13 ][TFSI] through the lens of hole theory. Here by exchanging the chloride anion of ethaline with bis(trifuoromethylsufonyl)imide) anion, [TFSI], is found to increase the fluidity of the solvent and promoted stronger solute-solvent interactions. Raman spectroscopy shows that methylviologen is solvated by EG in solutions with ChTFSI but interstitially accommodated in void spaces in both ChCl solutions and [PYR 13 ][TFSI]. Here, complex solvents such as DESs are regarded as “designer solvents” and we demonstrate that the physical properties and solvation characteristics of these fluids can be varied by anion exchange of hydrogen bond acceptors.

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

Chromium Tolerant, Highly Active and Stable Electrocatalytic Internal Surface Coating for Cathode of Commercial SOFCs (Final Report)

This project is aimed to develop a chromium (Cr) tolerant, highly active, and stable coating layer on the internal surfaces of the porous composite cathode from commercially available SOFCs. Such coating layer was developed using the additive manufacturing process of Atomic Layer Deposition (ALD) and has been applied on the cathode consisting of either an electronic conductor of LaxSr 1-x MnyO 3-δ (LSM) or mixed ionic and electronic conducting La x Sr 1-x Co y Fe 1-y O 3-δ (LSCF). PI's work has demonstrated that the internal surface of cathode from the commercial cells, can be further tailored using ALD coating to dramatically enhance the cell performance. For instance, ALD layer consisting heterostructured nano composite of nano-Pt and nano-(Mn 0.8 Co 0.2 ) 3 O 4 oxide on the internal surface of porous LSM/YSZ cathode from SOFCs, has resulted in the large reduction of the cell polarizations resistance by up to 53%, and enormous increase of power density over 370%. For the cells with LSCF/Sm 2 O 3 doped CeO 2 (SDC) cathode, the conformal layer of nano-composite consisting of superjacent CoOx and subjacent minimum amount of Pt nano-grains has resulted in the power density enhancement by 126% for the large scale industry tubular cells at 750°C, and both the performance enhancement and nanostructure of the ALD layer are stable over ~ 2000 h continuous operation performed at industry test station. In the meanwhile, those ALD coating layer developed by PI's work is also inherently Cr-tolerant, and could act as physical barrier for preventing Cr diffusion into the cathode backbone, so as to mitigate the Cr poisoning effect on the cathode. In this project, the impact of Cr on the performance of those ALD coated commercial cells has been evaluated. Based on evolution of the cell performance, the ALD coating layer chemistry and ALD coating layer thickness has been optimized to maximize the overall Cr tolerance, cell power density and cell longevity. Different ALD coating has been applied onto the internal surface of LSM/YSZ and LSCF/SDC backbone respectively. The architecture/scaffold structures on the internal surface of different cathode, designed by this project, was catalogued and analyzed using High Resolution Transmission Electron Microscopy (HRTEM), and cell power/durability performance are assured via comprehensive electrochemical performance testing in the industry operation relevant conditions. The impact of the electrochemical operation current density, the water humidity, the cell operation temperature, and cell operation duration on the Cr tolerance of ALD coated cells has been systematically investigated. There is completely different nanostructure degradation mechanisms between LSM and LSCF cells induced by Cr contamination. For the LSCF/SDC baseline cell, With the Cr source, there is no apparent Sr surface segregation phase even for the baseline cell operated for 3000 h at 750 °C. With the Cr source, there is significant amorphous (SrCr)Ox phase accumulated in the original pore region. For the commercial baseline cells, Cr contaminants on the LSM electrode severely impacted the entire cell's electrochemical performance and nanostructure degradation. Those degradations include (1). Peak power density loss of 64 % after 109 h of operation. The dramatic increase in Rp (2). They are cracking at LSM/SSZ interface, LSM grains. SSZ remains intact but with (CrMn)Ox. By contrast, ALD coating (MnCo)Ox/Pt dramatically improves the Cr resistance, as follows (1). ALD-coated cell with a power density is 280-380 % of the baseline cell, depending on the ALD layer thickness. (2). For a cell with a 20 nm thick ALD layer, there is a large performance enhancement (> 200 % power density) induced by ALD coating of Cr-tolerant Mn 0.8 Co 0.2 Ox. (3). For a cell with a 20 nm thick ALD layer, after 168 h at 750 °C power density of the ALD-coated cell is ~ 600% of that baseline cell upon operation with Cr contamination for 109 h. The ALD coating on the internal surface of cathode developed by this project integrated multi-functions. Those multi-functions include (1). Dramatically improving the cell power density for the commercial cells; (2). Dramatically improving contamination resistance of the cathode, for being an excellent protection coating layer sealing off Cr contamination. (3). Dramatically increasing the cell longevity by potentially preventing the microstructure evolution and grain coarsening of the cathode. Overall, this project will provide a simple solution to simultaneously enhance power density and increase the reliability, robustness, and endurance of commercial SOFCs, over the entire operating temperature range of 650-800 °C. For the inherently functional SOFC, the ALD coating of LSM based cathode mitigate the Cr-contamination. Power density of ALD-coated cell is ~ 600% of that baseline cell upon operation with Cr contamination. In addition to SOFCs, the novel on-demand design approach and creation of multifunctional heterogeneous architecture on the electrode surface presented in this work opens further research for their application in other types of fuel cells, batteries, and sensors for which electrochemical reactions on the surface are similarly critical.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nanoconfinement of Carbon Dioxide within Interfacial Aqueous/Ionic Liquid Systems

Nanoporous, gas-selective membranes have shown encouraging results for the removal of CO 2 from flue gas, yet the optimal design for such membranes is often unknown. Therefore, we used molecular dynamics simulations to elucidate the behavior of CO 2 within aqueous and ionic liquid (IL) systems ([EMIM][TFSI] and [OMIM][TFSI]), both confined individually and as an interfacial aqueous/IL system. We found that within aqueous systems the mobility of CO 2 is reduced due to interactions between the CO 2 oxygens and hydroxyl groups on the pore surface. Within the IL systems, we found that confinement has a greater effect on the [EMIM][TFSI] system as opposed to the [OMIM][TFSI] system. Paradoxically, the larger and more asymmetrical [OMIM]+ molecule undergoes less efficient packing, resulting in fewer confinement effects. Finally, free energy surfaces of the nanoconfined aqueous/IL interface demonstrate that CO 2 will transfer spontaneously from the aqueous to the IL phase.

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

Dynamical correlation between ions in liquid seen in real space and time

Controlling ionic transport in liquids is anticipated to provide new scientific and technological opportunities, but it requires accurate knowledge of atomic-scale dynamics of ions beyond the hydrodynamic description. Atomic dynamics in liquids is characterized by strong and dynamical correlations among atoms, which render a conventional approach to describing the dynamics in reciprocal space challenging. We propose an alternative approach to describe liquid dynamics in real space and time using the time-delayed two-body correlation function, Van Hove correlation function, and discuss recent results on the local dynamics in water and electrolytes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Simulations with machine learning potentials identify the ion conduction mechanism mediating non-Arrhenius behavior in LGPS

Abstract Li 10 Ge(PS 6 ) 2 (LGPS) is a highly concentrated solid electrolyte, in which Coulombic repulsion between neighboring cations is hypothesized as the underlying reason for concerted ion hopping, a mechanism common among superionic conductors such as Li 7 La 3 Zr 2 O 12 (LLZO) and Li 1.3 Al 0.3 Ti 1.7 (PO 4 ) 3 (LATP). While first principles simulations using molecular dynamics (MD) provide insight into the Li + transport mechanism, historically, there has been a gap in the temperature ranges studied in simulations and experiments. Here, we used a neural network potential trained on density functional theory (DFT) simulations, to run up to 40-nanosecond long MD simulations at DFT-like accuracy to characterize the ion conduction mechanisms across a range of temperatures that includes previous simulations and experimental studies. We have confirmed a Li + sublattice phase transition in LGPS around 400 K, below which the ab -plane diffusivity D a b ∗ is drastically reduced. Concomitant with the sublattice phase transition near 400 K, there is less cation-cation (cross) correlation, as characterized by Haven ratios closer to 1, and the vibrations in the system are more harmonic at lower temperature. Intuitively, at high temperature, the collection of vibrational modes may be sufficient to drive concerted ion hops. However, near room temperature, the vibrational modes available may be insufficient to overcome electrostatic repulsion, thus resulting in less correlated ion motion and comparatively slower ion conduction. Such phenomena of a sublattice phase transition, below which concerted hopping plays a less significant role, may be extended to other highly concentrated solid electrolytes such as LLZO and LATP.

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

Metallization of colloidal crystals

Colloidal crystals formed by size-asymmetric binary particles coassemble into a wide variety of colloidal compounds with lattices akin to ionic crystals. Recently, a transition from a compound phase with a sublattice of small particles to a metal-like phase in which the small particles are delocalized has been predicted computationally and observed experimentally. In this colloidal metallic phase, the small particles roam the crystal maintaining the integrity of the lattice of large particles, as electrons do in metals. A similar transition also occurs in superionic crystals, termed sublattice melting. Here, we use energetic principles and a generalized molecular dynamics model of a binary system of functionalized nanoparticles to analyze the transition to sublattice delocalization in different coassembled crystal phases as a function of temperature (T), number of grafted chains on the small particles, and number ratio between the small and large particles n s : n l . Here we find that n s : n l is the primary determinant of crystal type due to energetic interactions and interstitial site filling, while the number of grafted chains per small particle determines the stability of these crystals. We observe first-order sublattice delocalization transitions as T increases, in which the host lattice transforms from low- to high-symmetry crystal structures, including A20 → bct → bcc, A d → bct → bcc, and bcc → bcc/fcc → fcc transitions and lattices. Analogous sublattice transitions driven primarily by lattice vibrations have been seen in some atomic materials exhibiting an insulator-metal transition also referred to as metallization. We also find minima in the lattice vibrations and diffusion coefficient of small particles as a function of n s : n l , indicating enhanced stability of certain crystal structures for n s : n l values that form compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Determination of Molecular Structure and Dynamics of Molten Salts by Advanced Neutron and X-ray Scattering Measurements and Computer Modeling

The design and development of fully functional Molten-Salt Reactors (MSR) require detailed knowledge of the molten salt properties in order to understand and predict the salt’s behavior. Fundamental properties of interest include molecular structure, speciation, and dynamics (such as diffusion coefficients) of salt components and dissolved corrosion and fission products. Computer modeling is necessary to predict changes in physical and chemical properties due to irradiation, burning of dissolved fuel, and corrosion. The modeling requires experimental data, and advanced neutron and x-ray scattering and spectroscopy provide the most reliable and direct determination of the structure (Pair-Distribution Functions, PDF), and dynamics of ions in the melt. This project dealt with both fluoride and chloride salts. The PDFs have been measured by a combination of neutron and x-ray diffraction. We utilized the techniques of isotope substitutions, a very powerful tool available for neutron-scattering, to extract the details of the liquid structure. Although similar measurements have been done before, modern advanced neutron and x-ray-scattering techniques allow collecting the data at much higher resolution and in a wider range of temperatures. Importantly, we were among the first to study fluoride salts by neutron scattering. The importance of impurities and their effects on salt properties have become apparent recently and so new methods of salt purification were developed. We took advantage of these developments to produce reliable data, which have been used for computer simulations of both clean salts and those with added fission and corrosion products most relevant for MSRs. Ab initio molecular dynamics simulations have been performed to understand the multi-component liquid solution, in particular solubility of impurities and thermodynamic interactions in relation to the ionic-cluster structure of the fluid. We applied machine learning to regress from the simulation and experimental data in order to develop a fast-acting model that can handle molten salt with an arbitrary (≥ 10) number of chemical elements and be able to predict chemical potential as a function of composition and temperature. This project resulted in a number of experimental and computer-simulation publications, a patent application, and numerous conference presentations (American Physical Society, American Chemical Society, and The Electrochemical Society among others). Multiple students and postdocs participated and collaborated on aspects of this project. This project seeded new collaborations between MIT and other institutions, such as the University of Massachusetts Lowell, the University of Illinois Urbana-Champaign, the University of California Berkeley, and Oak Ridge and Los Alamos National Labs. As such, this project has had a broad and lasting impact beyond its original scientific scope.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code was efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators.

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

Investigating intermolecular interactions among CO 2 , water and PEEK-ionene membrane using cryo ToF-SIMS and isotopic labeling

Cryogenic time-of-flight secondary ion mass spectrometry (cryo ToF-SIMS) has emerged as a powerful tool for investigating molecular interactions, speciation, and dynamics in materials for CO 2 capture. In this study, we apply cryo ToF-SIMS to probe interactions between CO 2 , water, and PEEK-ionene membranes—a promising material for direct CO 2 capture due to its selectivity, durability, and efficiency. Despite this potential, the mechanisms governing CO 2 diffusion and the influence of water vapor on CO 2 behavior remain unclear. To address this, we loaded PEEK-ionene membranes with 13 CO 2 and D 2 O and employed cryo ToF-SIMS to visualize the 3D distribution of CO 2 and water within the membrane. While prior studies suggest that 13 CO 2 is absorbed under ambient conditions, our cryo ToF-SIMS analysis revealed no enhancement of the 13 C/ 12 C ratio, suggesting weak CO 2 -membrane interactions. As a result, CO 2 vaporizes even at low temperatures (−140°C) under vacuum conditions. In contrast, D 2 O displayed a relatively homogeneous distribution in the membrane, suggesting stronger water-membrane interactions via hydrogen bonding (18–20 kJ/mol). Interestingly, CO 2 was not detected in D 2 O-loaded membranes, indicating minimal interference from water vapor on CO 2 diffusion. As a comparison, the cryo ToF-SIMS data show that 13 CO 2 can readily react with a basic Na 2 CO 3 aqueous solution to form NaH 13 CO 3 . These findings demonstrate cryo ToF-SIMS as a critical technique for understanding gas-water-membrane interactions, offering insights for membrane functionalization to improve CO 2 capture efficiency.

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

Critical Role of Anion–Solvent Interactions for Dynamics of Solvent-in-Salt Solutions

Most polar solvent molecules are unstable toward electrode materials used in Li-based batteries. Solid electrolytes and ionic liquids are far more stable; however, they have relatively low conductivity, and therefore electrical energy storage devices based on them would suffer from low power. Solvent-in-salt (SIS) systems combine chemical stability with relatively high conductivity. Here, we show how the nature of the employed anion affects the structure and dynamics of SIS systems. The transport of ions in lithium bis(fluorosulfonyl)imide (Li-FSI) systems was determined to be always faster than that in lithium bis(trifluoromethanesulfonyl)imide (Li-TFSI) systems. Moreover, we found that viscosity does not solely control conductivity and that the lower conductivity of TFSI – solutions is related to their stronger interaction with the solvent. This restricts solvent dynamics and slows down ion motions compared to that of FSI – . Interestingly, the TFSI–solvent interaction also leads to better charge separation (weaker ion–ion correlations) and a higher transference number for Li. Finally, our results suggest that the ability to tune the solvent network formed around the anions may further improve electrolyte conductivity and Li transference number for safer and more efficient energy storage devices.

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

Theory of Cation Solvation in the Helmholtz Layer of Li-Ion Battery Electrolytes

The solvation environments of Li + in conventional nonaqueous battery electrolytes, such as LiPF 6 in mixtures of ethylene carbaronate (EC) and ethyl methyl carbonate (EMC), are often used to rationalize transport properties and solid electrolyte interphase (SEI) formation. Solvation environments in the compact electrical double layer (EDL) next to the electrode, also known as the Helmholtz layer, determine (partially) what species can react to form the SEI, with bulk solvation environments often being used as a proxy. Here, we develop and test a theory of cation solvation in the Helmholtz layer of nonaqueous Li-ion battery electrolytes. First, we validate the theory against bulk and diffuse EDL atomistic molecular dynamics (MD) simulations of LiPF 6 EC/EMC mixtures as a function of surface charge, where we find the theory can qualitatively capture the solvation environments. Next, we turn to the Helmholtz layer, where we find the main effect of the solvation structures next to the electrode is an apparent reduction in the number of binding sites between Li + and the solvents, again where we find reasonable agreement with our developed theory. Finally, by solving a simplified version of the theory, we find that the probability of Li + binding to each solvent remains equal to the bulk probability, suggesting that the bulk solvation environments are a reasonable place to start when understanding battery electrolytes. Our developed formalism can be parametrized from bulk MD simulations and used to predict the solvation environments in the Helmholtz layer through reducing the number of available coordination sites, which can be used to determine what could react and form the SEI.

Helmholtz↗
DOE OSTI2022

Lithium Solvation and Mobility in Ionic Liquid Electrolytes with Asymmetric Sulfonyl-Cyano Anion

The solvation structure and transport properties of Li + in ionic liquid (IL) electrolytes based on n-methyl-n-butylpyrrolidinium cyano(trifluoromethanesulfonyl)imide [PYR 14 ][CTFSI] and [Li][CTFSI] (0 ≤ x Li ≤ 0.7) were studied by Raman and Nuclear Magnetic Resonance (NMR) diffusometry, and molecular dynamics (MD) simulations. At x Li < 0.3, Li + coordination is dominated by the cyano group. As x Li is increased, free cyano-sites become limited, resulting in increased coordination via the sulfonyl group. Here, the 1:1 mixture of the symmetric anions bis(trifluoromethanesulfonyl)imide ([TFSI]) and dicyanamide ([DCA]) results in similar physical properties as the IL with [CTFSI]. However, anion asymmetry is shown to increase Li-salt solubility and promote Li+ transference. The lifetimes of Li + -cyano coordination for [CTFSI] are calculated to be shorter than those for [DCA], indicating that the competition from the sulfonyl group weakens its solvation with Li + . This resulted in higher Li + transference for the electrolyte with [CTFSI]. In relation to the utility of these electrolytes in energy storage, the Li–LiFePO 4 half cells assembled with IL electrolyte (x Li = 0.3, 0.5, and 0.7) demonstrated a nominal capacity of 140 mAh/g at 0.1C rate and 90 °C where the cell with x Li = 0.7 IL electrolyte demonstrated 61% capacity retention after 100 cycles and superior rate capability owing to increased electrochemical stability.

25 ENERGY STORAGE↗
DOE OSTI2024

Review of the second charged-particle transport coefficient code comparison workshop

We report the results of the second charged-particle transport coefficient code comparison workshop, which was held in Livermore, California on 24–27 July 2023. This workshop gathered theoretical, computational, and experimental scientists to assess the state of computational and experimental techniques for understanding charged-particle transport coefficients relevant to high-energy-density plasma science. Data for electronic and ionic transport coefficients, namely, the direct current electrical conductivity, electron thermal conductivity, ion shear viscosity, and ion thermal conductivity were computed and compared for multiple plasma conditions. Additional comparisons were carried out for electron–ion properties such as the electron–ion equilibration time and alpha particle stopping power. Overall, 39 participants submitted calculated results from 18 independent approaches, spanning methods from parameterized semi-empirical models to time-dependent density functional theory. In the cases studied here, we find significant differences—several orders of magnitude—between approaches, particularly at lower temperatures, and smaller differences—roughly a factor of five—among first-principles models. We investigate the origins of these differences through comparisons of underlying predictions of ionic and electronic structure. The results of this workshop help to identify plasma conditions where computationally inexpensive approaches are accurate, where computationally expensive models are required, and where experimental measurements will have high impact.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Pauling’s rules for oxide-based minerals: A re-examination based on quantum mechanical constraints and modern applications of bond-valence theory to Earth materials

Since their introduction in 1929, Pauling’s five rules have been used by scientists from many disciplines to rationalize and predict stable arrangements of atoms and coordination polyhedra in crystalline solids; amorphous materials such as silicate glasses and melts; nanomaterials, poorly crystalline solids; aqueous cation and anion complexes; and sorption complexes at mineral-aqueous solution interfaces. The predictive power of these simple yet powerful rules was challenged recently by George et al. (2020), who performed a statistical analysis of the performance of Pauling’s five rules for about 5000 oxide crystal structures. They concluded that only 13% of the oxides satisfy the last four rules simultaneously and that the second rule has the most exceptions. They also found that Pauling’s first rule is satisfied for only 66% of the coordination environments tested and concluded that no simple rule linking ionic radius to coordination environment will be predictive due to the variable quality of univalent radii. We address these concerns and discuss quantum mechanical calculations that complement Pauling’s rules, particularly his first (radius sum and radius ratio rule) and second (electrostatic valence rule) rules. We also present a more realistic view of the bonded radii of atoms, derived by determining the local minimum in the electron density distribution measured along trajectories between bonded atoms known as bond paths, i.e., the bond critical point (rc). Electron density at the bond critical point is a quantum mechanical observable that correlates well with Pauling bond strength. Moreover, a metal atom in a polyhedron has as many bonded radii as it has bonded interactions, resulting in metal and O atoms that may not be spherical. O atoms, for example, are not spherical in many oxide-based crystal structures. Instead, the electron density of a bonded oxygen is often highly distorted or polarized, with its bonded radius decreasing systematically from ~1.38 Å when bonded to highly electropositive atoms like sodium to 0.64 Å when bonded to highly electronegative atoms like nitrogen. Bonded radii determined for metal atoms match the Shannon (1976) radii for more electropositive atoms, but the match decreases systematically as the electronegativities of the M atoms increase. As a result, significant departures from the radius ratio rule in the analysis by George et al. (2020) is not surprising. We offer a modified, more fundamental version of Pauling’s first rule and demonstrate that the second rule has a one-to-one connection between the electron density accumulated between the bonded atoms at the bond critical point and the Pauling bond strength of the bonded interaction. Pauling’s second rule implicitly assumes that bond strength is invariant with bond length for a given pair of bonded atoms. Many studies have since shown that this is not the case, and Brown and Shannon (1973) developed an equation and a set of parameters to describe the relation between bond length and bond strength, now redefined as bond valence to avoid confusion with Pauling bond-strength. Brown (1980) used the valence-sum rule, together with the path rule and the valence-matching principle, as the three axioms of bond-valence theory (BVT), a powerful method for understanding many otherwise elusive aspects of crystals and also their participation in dynamic processes. We show how a priori bond-valence calculations can predict unstrained bond-lengths and how bond-valence mapping can locate low-Z atoms in a crystal structure (e.g., Li) or examine possible diffusion pathways for atoms through crystal structures. In addition, we briefly discuss Pauling’s third, fourth, and fifth rules, the first two of which concern the sharing of polyhedron elements (edges and faces) and the common instability associated with structures in which a polyhedron shares an edge or face with another polyhedron and contains high-valence cations. The olivine [α-(MgxFe1–x)2SiO4] crystal structure is used to illustrate the distortions from hexagonal close-packing of O atoms caused by metal-metal repulsion across shared polyhedron edges. We conclude by discussing several applications of BVT to Earth materials, including the use of BVT to: (1) locate H+ ions in crystal structures, including the location of protons in the crystal structures of nominally anhydrous minerals in Earth’s mantle; (2) determine how strongly bonded (usually anionic) structural units interact with weakly bonded (usually cationic) interstitial complexes in complex uranyl-oxide and uranyl-oxysalt minerals using the valence-matching principle; (3) calculate Lewis acid strengths of cations and Lewis base strengths of anions; (4) determine how (H2O) groups can function as bond-valence transformers by dividing one bond into two bonds of half the bond valence; (5) help characterize products of sorption reactions of aqueous cations (e.g., Co2+ and Pb2+) and oxyanions [e.g., selenate (Se6+O4)2- and selenite (Se4+O3)2-] at mineral-aqueous solution interfaces and the important role of protons in these reactions; and (6) help characterize the local coordination environments of highly charged cations (e.g., Zr4+, Ti4+, U4+, U5+, and U6+) in silicate glasses and melts.

Geochemistry & Geophysics↗
DOE OSTI2025

Tuning Solvation Dynamics of Electrolytes at Their Eutectic Point Through Halide Identity

Deep eutectic solvents (DESs) are regarded as highly promising solvent systems for redox flow batteries. DESs, composed of choline halides (ChX, X = F − , Cl − , Br − , I − ) and ethylene glycol (EG), exhibit distinct physicochemical properties at their eutectic points, including halide-dependent phase behavior, viscosity, polarity, conductivity, and solvation dynamics. In this study, we investigate the effects of the halide identity on the solvation properties of ChX:EG mixtures at varying mol % of ChX salt content. The solvatochromic polarity based on ET(30) measurements indicates higher polarity for larger halides (I − > Br − ) than for smaller halides (Cl − > F − ), which exhibit larger compensating solvation shells. The ionic conductivity follows the trend of the solvent fluidity (the inverse of the viscosity), namely ChCl > ChBr > ChI > ChF, influenced by the ion mobility and solvodynamic radii. Measurements of the liquidus temperatures (T L ) reveal that the system with ChCl exhibits the deepest eutectic point (at ~20 mol % ChCl), while ChBr and ChI have shallower minima at ~10 mol % ChBr and ~3 mol % ChI, respectively. ChF does not display a eutectic transition but instead appears to readily supercool at salt concentrations above 30 mol % ChF. Consistent with the phase transition measurements, femtosecond transient absorption spectroscopy shows that in the ChCl system, the solvation dynamics become faster with an increasing salt concentration up to ~16.67 mol %, after which the dynamics slow down with further increases in the salt content. The ChF-based system exhibits similar behavior, though with slower dynamics. In contrast, the solvation dynamics of the systems containing ChBr and ChI monotonously slow down with an increasing salt concentration, in agreement with the phase transition measurements, which show that the eutectic points occur at low salt concentrations. These measurements suggest that the solvent composition and, in particular, the identity of the halide anion play a significant role in the solvation behavior of these ethylene-glycol-based DESs, offering a foundation for tuning the DES properties for specific applications.

Biochemistry & Molecular Biology↗