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↗
NASA NTRS1995

Supersaturated Electrolyte Solutions: Theory and Experiment

Highly supersaturated electrolyte solutions can be prepared and studied employing an electrodynamic levitator trap (ELT) technique. The ELT technique involves containerless suspension of a microdroplet thus eliminating dust, dirt, and container walls which normally cause heterogeneous nucleation. This allows very high supersaturations to be achieved. A theoretical study of the experimental results obtained for the water activity in microdroplets of various electrolyte solutions is based on the development of the Cahn-Hilliard formalism for electrolyte solutions. In the approach suggested the metastable state for electrolyte solutions is described in terms of the conserved order parameter omega(r,t) associated with fluctuations of the mean solute concentration n(sub 0). Parameters of the corresponding Ginzburg-Landau free energy functional which defines the dynamics of metastable state relaxation are determined and expressed through the experimentally measured quantities. A correspondence of 96-99 % between theory and experiment for all solutions studied was achieved and allowed the determination of an analytical expression for the spinodal concentration n(sub spin), and its calculation for various electrolyte solutions at 298 K. The assumption that subcritical solute clusters consist of the electrically neutral Bjerrum pairs has allowed both analytical and numerical investigation of the number-size N(sub c) of nucleation monomers (aggregates of the Bjerrum pairs) which are elementary units of the solute critical clusters. This has also allowed estimations for the surface tension Alpha, and equilibrium bulk energy Beta per solute molecule in the nucleation monomers. The dependence of these properties on the temperature T and on the solute concentration n(sub 0) through the entire metastable zone (from saturation concentration n(sub sat) to spinodal n(sub spin) is examined. It has been demonstrated that there are the following asymptotics: N(sub c), = I at spinodal concentration and N(sub c) = infinity at saturation.

Izmailov, Alexander F.↗
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↗
NASA NTRS2018

Exploring Li-Air Batteries for High Specific Energy and High Power Applications: A Simulation Study

Commercialization of lithium-air batteries face many challenges, such as electrolyte decomposition, short cycle life, low energy efficiency, low power density, etc. However, commercialization of Li-air batteries for mass sensitive applications such as electric vehicles, portable power source, and drones is more challenging due to additional constraints of safety, electrolyte evaporation, high specific energy requirements, and reliable discharge times. In this presentation, we will present our finite element simulation results comparing Li-O2 and Li-air batteries using power density, energy density, specific power, and discharge times as metrics to evaluate different electrolytes and electrode geometry to reduce total mass and maximize discharge current. We use a finite element model and a discharge product model developed in which is based on porous electrode and concentrated electrolyte theories and the discharge product is modeled using quantum tunneling; for reaction kinetics and oxygen diffusion, an improved model was used. The electrolyte properties such as ion conductivity and ion diffusion were obtained from Molecular Dynamics (MD) simulations while the other parameters for the finite element model were calibrated to match experiments at high discharge current densities (>1.5 mA/cm2). The mass densities of different electrolytes were computed using MD simulations as well. For this presentation, we examine the practical electrochemical mass of a system at different current ratings, the sensitivity of mass to the use of ambient air as compared to pure oxygen as well as the electrolyte, which affects maximum current density and the total mass associated with the electrolyte (which includes the mass of additional components), and, the optimization of battery geometry for total discharge time, average discharge voltage, maximum discharge current, and minimum electrochemical mass.

Mehta, M.↗
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↗
NASA NTRS1990

Theory of electrohydrodynamic instabilities in electrolytic cells

The paper develops the theory of the hydrodynamic stability of an electrolytic cell as a function of the imposed electric current. A new electrohydrodynamic instability is encountered when the current is forced to exceed the Nernst limit. The convection is driven by the volume force exerted by the electric field on space charges in the electrolyte. This intrinsic instability is found to be easily masked by extrinsic convection sources such as gravity or stirring. A linear stability analysis is performed and a dimensionless number Le is derived whose value determines the convection pattern.

Bruinsma, R.↗
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↗
NASA NTRS2019

A Multi-Physics Study on High-Specific Power Li-O2 Batteries for Electric Aircraft

Commercialization of lithium-air batteries faces many challenges, such as electrolyte decomposition, short cycle life, low energy and power density, etc. However, commercialization of Li-O2 batteries for aeronautics is much more challenging due to additional safety constraints on cyclability and performance (high specific power and specific energy). For this presentation, we will discuss inter-related aspects of physics-based modeling of a pack: cell and battery model calibration. In addition, we will evaluate and present optimal battery designs for high discharge current density, high discharge time, and low battery mass using simulation-based optimization.The Finite Element Model (FEM) used to simulate a Li-O2 cell is based on the work of Bevara [1]. The different aspects of the model are based on: porous electrode theory and concentrated electrolyte theory; quantum tunneling model for the resistance of conformal layer of discharge product (Li2O2) [1]; Butler-Volmer kinetics for electrochemical reaction; Fick's diffusion for oxygen transport; and an oxygen dissolution model is applied at the air/electrolyte interface [2]. The electrolyte properties such as ion conductivity, ion diffusion, oxygen diffusion, and mass density of the electrolyte were taken from Molecular Dynamics (MD) simulations [3]; while the other model parameters, which includes mass of cell components, were calibrated to match experiments at high discharge current densities. The cell mass includes the anode, cathode, separator, electrolyte, and other components (such as current collector). This calibrated model is used to perform parametric studies on cathode thickness, porosity, tortuosity, carbon particle size, electrolyte transport and material properties, partial pressure of oxygen, discharge time, and discharge current density to study optimal designs for high specific power and energy. References:1. Bevara, V. & Andrei, P. (2014), J. Electrochem. Soc. 161 (14), A2068-A2079.2.Mehta, M. & Andrei, P. (2015), J. Power Sources. 286, 299-308.3.Liyana-Arachchi, T.; Haskins, J.; Burke, C.; Diederichsen, K.; McCloskey, B.; & Lawson, J. (2018), J. Phys. Chem. B. 122 (36), 8548 - 8559.4.Choi, W.; Kikumoto, H.; Choudhary, R. & Ooka, R. (2018), Applied Energy, 209, 306-321.

Mehta, Mohit↗
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↗
NASA NTRS1970

Solubility of non-polar gases in electrolyte solutions

Solubility theory describes the effects of both concentration and temperature on solute activity coefficients. It predicts the salting-out effect and the decrease in solubility of non-polar gases with increased electrolyte concentration, and can be used to calculate heats of solution, entropies, and partial molal volumes of dissolved gases

Walker, R. L., Jr.↗
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↗