Search NASA⌕ Search

SEARCH · Search NASA

Results for “electrolyte theory”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Multiple Operando Fields Can Identify a Predictive Mass Transport Theory in Electrolytes

An electrolyte transport theory connects its transport properties, evolution of spatiotemporal fields (e.g., concentration), and corresponding macroscopic current and voltage responses. Given this interconnection, the transport properties are typically inferred by analyzing the macroscopic response through the lens of a chosen electrolyte transport theory. Unfortunately, the same macroscopic measurements can be analyzed by using different theories to arrive at seemingly dissimilar properties that are inconsistent with each other. Here, we offer a resolution to this dilemma by analyzing multiple spatiotemporal (i.e., operando) fields for a given electrolyte. We show that predictive transport theory can analyze operando fields to estimate underlying transport properties and subsequently predict another operando field. A theory passing this test identifies meaningful transport properties such that the electrolyte behavior can be accurately predicted over a wide range of excitations, and is critical to property-based electrolyte screening and electrolyte discovery efforts.

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

A Systematic Way to Extend the Debye–Hückel Theory beyond Dilute Electrolyte Solutions

An extended Debye–Hückel theory with fourth order gradient term is developed for electrolyte solutions; namely, the electric potential φ(r) of the bulk electrolyte solution can be described by ∇ 2 φ(r) = κ 2 φ(r) + L Q 2 ∇ 4 φ(r), where the parameters κ and L Q are chosen to reproduce the first two roots of the dielectric response function of the bulk solution. Three boundary conditions for solving the electric potential problem are proposed based upon the continuity conditions of involving functions at the dielectric boundary, with which a boundary element method for the electric potential of a solute with a general geometrical shape and charge distribution is derived. Solutions for the electric potential of a spherical ion and a diatomic molecule are found and used to calculate their electrostatic solvation energies. In this work, the validity of the theory is successfully demonstrated when applied to binary as well as multicomponent primitive models of electrolyte solutions.

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

Beyond the Debye–Hückel limit: Toward a general theory for concentrated electrolytes

The phenomenon of underscreening in concentrated electrolyte solutions leads to a larger decay length of the charge–charge correlation than the prediction of Debye–Hückel (DH) theory and has found a resurgence of both theoretical and experimental interest in the chemical physics community. To systematically understand and investigate this phenomenon in electrolytes requires a theory of concentrated electrolytes to describe charge–charge correlations beyond the DH theory. We review the theories of electrolytes that can transition from the DH limit to concentrations where charge correlations dominate, giving rise to underscreening and the associated Kirkwood Transitions (KTs). In this perspective, we provide a conceptual approach to a theoretical formulation of electrolyte solutions that exploits the competition between molecular-informed short-range (SR) and long-range interactions. We demonstrate that all deviations from the DH limit for real electrolyte solutions can be expressed through a single function ΣQ that can be determined both theoretically and numerically. Importantly, ΣQ can be directly related to the details of SR interactions and, therefore, can be used as a tool to understand how differences in representations of interaction can influence collective effects. The precise function form of ΣQ can be inferred through a Gaussian field theory of both the number and charge densities. The resulting formulation is validated by experiment and can accurately describe the collective phenomenon of screening in concentrated bulk electrolytes. Importantly, the Gaussian field theory predictions of the screening lengths appear to be less than ∼1 nm at concentrations above KTs.

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

Complete characterization of a lithium battery electrolyte using a combination of electrophoretic NMR and electrochemical methods

Improving transport properties of electrolytes is important for developing lithium-ion batteries for future energy storage applications. In Newman's concentrated solution theory, electrolytes are characterized by three transport parameters, conductivity, diffusion coefficient, and transference number, in addition to the thermodynamic factor. In this work, these parameters are all determined for an exemplar liquid electrolyte, lithium bis(trifluoromethanesulfonyl)imide mixed in tetraethylene glycol dimethyl ether, using electrochemical methods. The intrinsic coupling between parameters obtained by electrochemical methods results in large error bars in the transference number that obscure the transport behavior of the electrolyte. Here, we use electrophoretic NMR (eNMR) to measure the electric-field-induced ion and solvent velocities to obtain the transference number directly, which enables determination of the thermodynamic factor with greater certainty. In conclusion, our work indicates that the combination of eNMR and electrochemical methods provides a robust approach for complete characterization of battery electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2021

Mass action model of solution activity via speciation by solvation and ion pairing equilibria

Abstract Solutes and their concentrations influence many natural and anthropogenic solution processes. Electrolyte and solution models are used to quantify and predict such behavior. Here we present a mechanistic solution model based on mass action equilibria. Solvation and ion pairing are used to model speciated solute and solvent concentrations such that they correlate to a solution’s vapor pressure (solvent activity) according to Raoult’s law from dilute conditions to saturation. This model introduces a hydration equilibrium constant ( K ha ) that is used with either an ion dissociation constant ( K id ) or a hydration modifier ( m ) with an experimentally determined ion dissociation constant, as adjustable parameters to fit vapor–liquid equilibrium data. The modeled solvation equilibria are accompanied by molecular dynamics (MD) studies that support a decline in the observed degree of solvation with increased concentration. MD calculations indicate this finding is a combination of a solvent that solvates multiple solutes, and changes in a solute’s solvation sphere, with the dominant factor changing with concentration. This speciation-based solution model is lateral to established electrostatics-based electrolyte theories. With its basis in mass action, the model can directly relate experimental data to the modeled solute and solvent speciated concentrations and structures.

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

Solvation Environment and Interface Dynamics of Li 2 B 12 H 12 and Li 2 B 12 F 12 Electrolytes Uncovered by Theory and Operando Optical and FTIR Spectroelectrochemistry

In this work, we evaluated two closo-borate salts (Li 2 B 12 H 12 and Li 2 B 12 F 12 ) in propylene carbonate from theoretical and experimental perspectives to understand how the coordination environment influences their spectroscopic and electrochemical properties. The coordination environments of the closo-borate salts were modeled via density functional theory (DFT) and molecular dynamics (MD). Vibrational spectra calculated from the predicted coordination environments are in agreement with experimentally measured steady-state FTIR data. This theoretical investigation also suggested that Li 2 B 12 F 12 would possess a higher ionic conductivity than Li 2 B 12 H 12 , which was corroborated experimentally. Additionally, an electrochemical cell was designed and fabricated that enabled operando optical and FTIR spectroelectrochemical (OP-IR-SEC) measurements. This allowed for the simultaneous measurement of the relative changes of species at a lithium electrode-liquid electrolyte interface and the visualization of lithium plating at the electrode surface. This technique could provide new chemical insights and potentially link optical changes at the electrode–electrolyte interface to specific chemical species in similar electrochemical systems. The Li 2 B 12 F 12 electrolyte was found to have a higher thermal stability, which may find utility in applications for batteries that are subject to high-temperature conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comparing Experimentally-Measured Sand’s Times with Concentrated Solution Theory Predictions in a Polymer Electrolyte

We compare the electrochemically measured Sand’s time, the time required for the cell potential to diverge when the applied current density exceeds the limiting current, with theoretical predictions for a 0.47 M poly(ethylene oxide) (5 kg mol −1 )/LiTFSI electrolyte. The theoretical predictions are made using concentrated solution theory which accounts for both concentration polarization and polymer motion, using independently measured parameters that depend on concentration, c : conductivity ( κ ), salt diffusion coefficient ( D ), cationic transference number with respect to the solvent velocity ( t + 0 ), thermodynamic factor 1 + dln f ± dln c , and partial molar volume of the salt ( V ̅ ); f ± is the mean molar activity coefficient of the salt. We find quantitative agreement between experimental data and theoretical predictions. We derive a generalized analytical expression for Sand’s time for electrolytes based on dilute solution theory. This expression correctly predicts the divergence of the Sand’s time at the limiting current, in agreement with experimental data and concentrated solution theory predictions. When the applied current is large compared to the limiting current, the analytical expression approaches the standard expression for Sand’s time used in the literature.

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

Elucidating and contrasting the mechanisms for Mg and Ca sulfate ion-pair formation with multi-level embedded quantum mechanics/molecular dynamics simulations

Here, solutions and minerals containing sulfate (SO 4 2- ), and Ca 2+ and Mg 2+ cations, are ubiquitous throughout the lithosphere and are significant components of seawater, thus presenting a prototypical system for the study of strong electrolytes and crystal nucleation mechanisms. However, despite their relative abundance, key questions remain unanswered about the most fundamental atomic-level steps of their mineralization pathways and aqueous dynamics. Here, we carry out enhanced sampling multi-level molecular dynamics (MD) embedded correlated wavefunction theory simulations to elucidate ion-pairing mechanisms for Mg–SO 4 and Ca–SO 4 in concentrated aqueous solution, accurately capturing effects arising from both structural dynamics and electron exchange–correlation. We predict contact-ion-pair formation to be barrierless and highly exoergic for Ca–SO 4 , in agreement with its minimal solubility, whereas for Mg–SO 4 , solvent-shared and contact ion pairs have similar free energies, qualitatively consistent with its higher solubility. Finally, we demonstrate that brief high-temperature pre-equilibration may be utilized to accelerate convergence of free energies in blue-moon-ensemble enhanced-sampling MD.

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

Identification of efficient and stable solid-solid electrochemical interfaces for Li-S batteries. Final Report

The objective of this work is to develop a better understanding of the interactions between battery components that drive the observed battery performance. Using computational tools at the atomistic level, we elucidate molecular, crystal, and amorphous structures, interfacial structures and reactions, and electrochemical behavior. We focus on Li metal anodes, liquid electrolytes, solid electrolytes, S cathodes, and metal-oxide cathodes, and their corresponding interfaces. Within liquid electrolytes, we study their bulk structure and dynamics, as a function of salt concentration. The analysis allows to understand the formation of localized high concentrated regions as the solvent is substituted by a weakly interacting diluent species. Analysis of the electron distribution yields the reactive behavior at interfaces. Moreover, we employ tools that allow the characterization of ion and electron transport at interfaces, during the processes of solvation and desolvation occurring while cycling. For Sulfur cathodes we elucidate synthesis processes and the role of carbon-sulfur and nitrogen-sulfur interactions not only to provide electronic conductivity but also to define alternative reaction pathways during discharge of Li-S batteries. Full-cell simulations offer new insights regarding the integration of all the battery components. Our studies also include solid electrolytes and their interfaces with Li metal, S-C, and metal-oxide electrodes.

25 ENERGY STORAGE↗
DOE OSTI2021

Deep potential generation scheme and simulation protocol for the Li 10 GeP 2 S 12 -type superionic conductors

We report solid-state electrolyte materials with superior lithium ionic conductivities are vital to the next-generation Li-ion batteries. Molecular dynamics could provide atomic scale information to understand the diffusion process of Li-ion in these superionic conductor materials. Here, we implement the deep potential generator to set up an efficient protocol to automatically generate interatomic potentials for Li 10 GeP 2 S 12 -type solid-state electrolyte materials (Li 10 GeP 2 S 12 , Li 10 SiP 2 S 12 , and Li 10 SnP 2 S 12 ). The reliability and accuracy of the fast interatomic potentials are validated. With the potentials, we extend the simulation of the diffusion process to a wide temperature range (300 K–1000 K) and systems with large size (~1000 atoms). Important technical aspects such as the statistical error and size effect are carefully investigated, and benchmark tests including the effect of density functional, thermal expansion, and configurational disorder are performed. The computed data that consider these factors agree well with the experimental results, and we find that the three structures show different behaviors with respect to configurational disorder. Our work paves the way for further research on computation screening of solid-state electrolyte materials.

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

Infrared spectral signatures of interfacial water at TiO 2 –electrolyte interfaces from deep potential molecular dynamics

Vibrational spectroscopy is a powerful tool for probing water at oxide–electrolyte interfaces, but its molecular interpretation can be challenging. Here, we employ deep potential long-range molecular dynamics simulations with layer-resolved spectral analysis to investigate the microscopic origins of the infrared (IR) response of water at the interface with anatase TiO 2 (101), a prototypical oxide surface. The calculated interfacial spectra exhibit characteristic modifications compared to bulk water IR spectra, including enhanced intensities, a red shifted and broadened stretching band, and a higher-frequency shoulder, in qualitative agreement with experiments. Spectral decomposition shows that these signatures originate mainly from the first interfacial water layer, dominated by surface-bound H 2 O at Ti 5C sites, with secondary contributions from the second layer. A moderate salt concentration (0.4 M NaCl) leaves both the interfacial structure and the spectra essentially unchanged, while tuning the pH strongly modulates the spectral intensity. We establish a scaling relation linking the spectral intensity to the surface water dissociation fraction and the dipole moment, both governed by interfacial electric fields. These findings provide a microscopic framework for interpreting IR spectra of oxide–electrolyte interfaces.

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

Localized high concentration electrolytes decomposition under electron-rich environments

Localized high concentration electrolytes have been proposed as an effective route to construct stable solid–electrolyte interphase (SEI) layers near Li-metal anodes. However, there is still a limited understanding of decomposition mechanisms of electrolyte components during SEI formation. In this work, we investigate reactivities of lithium bis(fluorosulfonyl)imide (LiFSI, salt), 1,2-dimethoxyethane (DME, solvent) and tris(2,2,2-trifluoroethyl)orthoformate (TFEO, diluent) species in DME+TFEO mixed solvents and 1M LiFSI/DME/TFEO solutions. By supplying excess of electrons into the simulation cell, LiFSI is initially reduced via a 4-electron charge transferring reaction yielding F - and N(SO 2 ) 2 3- . The local solvation environment has little effect on the subsequent TFEO reaction, which typically requires 6 |e| to decompose into F - , HCOO - , CH2CF and - OCH 2 CF 3 . Besides, TFEO dehydrogenation reaction mechanism under the attack of anions is also identified. Unlike salt and diluent, DME shows good stability with any excess of electrons. Here, the energetics of most relevant reactions are characterized. Most reactions are thermodynamically favorable with low activation barriers.

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

Uncovering the binding nature of thiocyanate in contact ion pairs with lithium ions

Ion pair formation is a fundamental molecular process that occurs in a wide variety of systems, including electrolytes, biological systems, and materials. In solution, the thiocyanate (SCN–) anion interacts with cations to form contact ion pairs (CIPs). Due to its ambidentate nature, thiocyanate can bind through either its sulfur or nitrogen atoms, depending on the solvent. This study focuses on the binding nature of thiocyanate with lithium ions as a function of the solvents using FTIR, 2D infrared spectroscopy (2DIR) spectroscopies, and theoretical calculations. The study reveals that the SCN– binding mode (S or N end) in CIPs can be identified through 2DIR spectroscopy but not by linear IR spectroscopy. Linear IR spectroscopy shows that the CN stretch frequencies are too close to one another to separate N- and S-bound CIPs. Moreover, the IR spectrum shows that the S–C stretch presents different frequencies for the salt in different solvents, but it is related to the anion speciation rather than to its binding mode. A similar trend is observed for the anion bend. 2DIR spectra show different dynamics for N-bound and S-bound thiocyanate. In particular, the frequency–frequency correlation function (FFCF) dynamics extracted from the 2DIR spectra have a single picosecond exponential decay for N-bound thiocyanate and a biexponential decay for S-bound thiocyanate, consistent with the binding mode of the anion. Lastly, it is also observed that the binding mode also affects the line shape parameters, probably due to the different molecular mechanisms of the FFCF for N- and S-bound CIPs.

25 ENERGY STORAGE↗
DOE OSTI2023

Molecular dynamics study on the interfacial properties of mixtures of monomers of polyvinylpyrrolidone (PVP)-based battery binders on graphene and graphite surfaces

This study investigates the behavior of two different mixtures of monomers of polyvinylpyrrolidone (PVP)-based battery binders, polyvinylpyrrolidone:polyvinylidene difluoride (PVP:PVDF) and polyvinylpyrrolidone:polyacrylic acid (PVP:PAA), at graphene and graphite interfaces using classical molecular dynamics simulations. The aim is to identify the best performing monomer binder blend and carbon-based material for the design of battery-optimized energy devices. The PVP:PAA monomer binder blend and graphite are found to have the best interaction energies, densification upon adsorption, and more ordered structure. The adsorption of both monomer binder blends is strongly guided by the higher affinity of PVP and PAA monomeric molecules for the surfaces compared to PVDF. The structure of adsorbed layers of PVP:PVDF monomer binder blend on graphene and graphite develops more quickly than PVP:PAA, indicating faster kinetics. This study complements a previous density functional theory study recently reported by our group and contributes to a better understanding of the nanoscopic features of relevant interfacial regions involving mixtures of monomers of PVP-based battery binders and different carbon-based materials. In conclusion, the effect of a blend of commonly used monomer binders on carbon-based materials is essential for obtaining tightly bound anode and cathode active materials in lithium-ion batteries, which is crucial for designing battery-optimized energy devices.

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

Influence of surface chemistry on Li nucleation energetics on graphene-based surfaces

Lithium metal is a promising high-capacity anode material for solid-state batteries, but it typically suffers from poor cyclability. Carbon scaffold hosts have the potential to improve this performance due to their high electronic conductivity and large surface area, which facilitates lithium-ion adsorption and desorption. Scaffold surface chemistry is known to significantly influence performance outcomes, but the details of these interactions are not fully understood. Here, this study employs first-principles simulations to explore lithium transport and nucleation on graphene anodes with various surface chemistries. Using enhanced sampling techniques, ab initio molecular dynamics, and density functional theory calculations, we find that although surface chemistry has a minimal impact on lithium interfacial transport, it influences surface nucleation significantly. Both heteroatom dopants and intrinsic defects lower the nucleation barrier, creating a more favorable environment for lithium nucleation compared to pristine graphene. In addition, our results reveal a complex interplay between surface lithium concentration, lithium transport, and nucleation kinetics. These findings highlight the potential of surface modifications to precisely control nucleation processes on carbon-based anodes and provide design guidance for reducing dendrite formation and improving the cycle life of solid-state batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Defect Thermodynamics and Transport Properties of Proton Conducting Perovskite Electrode and Electrolyte Materials Evaluated Based on Density Functional Theory Modeling

Both electron-rich and electron-poor perovskite oxides have been used in solid oxide cell applications as electrode and electrolyte materials. The rich oxygen defect chemistry and its coupling to temperature, hydrogen-steam or oxygen-steam gas pressure, or to the applied potentials creates enormous complexities for modeling performance and degradation of the materials. Herein, density functional theory-based thermodynamic modeling was carried out to describe the defect chemistry and transport properties of the proton-conducting electrolyte BaZr1-xYxO3-δ (x≤0.1) and of the triple-conducting perovskite (La,Ba)(Fe,M)O3-δ (M=Y and Zr). The defect thermodynamics of intrinsic point defects and the hydrogen-related defect reactions were solved in integrated defect models and further used to predict the Brouwer diagram and the transport properties of the functional perovskites. For the electron-poor electrolytes BaZr0.9Y0.1O3-δ, the developed model has been used to describe the experimental transport properties in the SOC operating conditions. Specifically, the roles played by the acceptor-bound holes and the intrinsic and hydrogen point defects upon the conductivities of holes, protons, and oxygen vacancies under the hydrogen-rich and oxygen-rich conditions at various humidity levels were demonstrated. A defect modeling tool was also developed for the triple-conducting perovskite (La,Ba)(Fe,M)O3-δ (M=Y and Zr) to examine magnetic effects and hydride defects in defect equilibria.

defect thermodynamics↗