SEARCH · Engineering Papers
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.
The composition of the lithium cation (Li + ) solvation shell in mixed linear and cyclic carbonate-based electrolytes has been re-examined using Born–Oppenheimer molecular dynamics and Li + (EC) n (DMC) m cluster calculations.
There is increasing evidence that the solvation structure of the active components in a liquid electrolyte solution strongly impacts the performance in electrochemical applications. In this work, the nanoscale solvation structures and dynamics of Mg(BH 4 ) 2 and Mg(TFSI) 2 dissolved in diglyme (DGM) at various concentrations and ratios of Mg(BH 4 ) 2 /Mg(TFSI) 2 were investigated using a combination of natural abundance 25 Mg NMR, quantum chemistry calculations of 25 Mg NMR chemical shifts, classical molecular dynamics (MD) calculations, and electrochemical performance tests. By mixing two competing Mg salts, we were able to reduce the strong covalent interactions between Mg 2+ and BH 4 – anions. A small increase is observed in the coordination number of Mg-TFSI and a significant increase in the interaction of Mg 2+ ions with glymes. Through a combination of NMR, DFT and MD simulations, various stable species around 1 nm in size were detected in the mixed salt solution, which play key roles in the enhanced electrochemical performance of the mixed electrolyte. It is established that for the neat Mg(TFSI) 2 in DGM electrolyte at dilute concentrations the TFSI - is fully dissociated from Mg 2+ . At higher concentrations, Mg 2+ and TFSI - are only partially dissociated as contact ion pairs are formed. In contrast, at 0.01 M Mg(BH 4 ) 2 (saturated concentration) in DGM, the first solvation shell of a Mg 2+ ion contains two BH 4 - anions and one DGM molecule, while the second solvation shell consists of five to six DGM molecules. An exchange mechanism between the solvation structures in the combined electrolyte containing both Mg(BH 4 ) 2 and Mg(TFSI) 2 in DGM was found to result in the observation of a single 25 Mg NMR peak. This exchange is responsible for an increase in uncoordinated anions, as well as improved stability and ionic conductivity as compared to single anion solution. Solvent molecule rearrangement and direct Mg-ion exchange between the basic solvation structures are hypothesized as likely reasons for the exchange. We elucidate that the solvent rearrangement is energetically much more favorable than direct Mg-ion hopping and is thus suggested as the dominant exchange mechanism.
Explore the source record for details and available documents.
We develop efficient handling of solvation forces in the multiscale method of multiple time step molecular dynamics (MTS-MD) of a biomolecule steered by the solvation free energy (effective solvation forces) obtained from the 3D-RISM-KH molecular theory of solvation (three-dimensional reference interaction site model complemented with the Kovalenko-Hirata closure approximation). To reduce the computational expenses, we calculate the effective solvation forces acting on the biomolecule by using advanced solvation force extrapolation (ASFE) at inner time steps while converging the 3D-RISM-KH integral equations only at large outer time steps. The idea of ASFE consists in developing a discrete non-Eckart rotational transformation of atomic coordinates that minimizes the distances between the atomic positions of the biomolecule at different time moments. The effective solvation forces for the biomolecule in a current conformation at an inner time step are then extrapolated in the transformed subspace of those at outer time steps by using a modified least square fit approach applied to a relatively small number of the best force-coordinate pairs. The latter are selected from an extended set collecting the effective solvation forces obtained from 3D-RISM-KH at outer time steps over a broad time interval. The MTS-MD integration with effective solvation forces obtained by converging 3D-RISM-KH at outer time steps and applying ASFE at inner time steps is stabilized by employing the optimized isokinetic Nosé-Hoover chain (OIN) ensemble. Compared to the previous extrapolation schemes used in combination with the Langevin thermostat, the ASFE approach substantially improves the accuracy of evaluation of effective solvation forces and in combination with the OIN thermostat enables a dramatic increase of outer time steps. We demonstrate on a fully flexible model of alanine dipeptide in aqueous solution that the MTS-MD/OIN/ASFE/3D-RISM-KH multiscale method of molecular dynamics steered by effective solvation forces allows huge outer time steps up to tens of picoseconds without affecting the equilibrium and conformational properties, and thus provides a 100- to 500-fold effective speedup in comparison to conventional MD with explicit solvent. With the statistical-mechanical 3D-RISM-KH account for effective solvation forces, the method provides efficient sampling of biomolecular processes with slow and/or rare solvation events such as conformational transitions of hydrated alanine dipeptide with the mean life times ranging from 30 ps up to 10 ns for “flip-flop” conformations, and is particularly beneficial for biomolecular systems with exchange and localization of solvent and ions, ligand binding, and molecular recognition.
Continuum solvation models enable efficient first principles calculations of chemical reactions in solution, but require extensive parametrization and fitting for each solvent and class of solute systems. Here, we examine the assumptions of continuum solvation models in detail and replace empirical terms with physical models in order to construct a minimally-empirical solvation model. Specifically, we derive solvent radii from the nonlocal dielectric response of the solvent from ab initio calculations, construct a closed-form and parameter-free weighted-density approximation for the free energy of the cavity formation, and employ a pair-potential approximation for the dispersion energy. We show that the resulting model with a single solvent-independent parameter: the electron density threshold (n{sub c}), and a single solvent-dependent parameter: the dispersion scale factor (s{sub 6}), reproduces solvation energies of organic molecules in water, chloroform, and carbon tetrachloride with RMS errors of 1.1, 0.6 and 0.5 kcal/mol, respectively. We additionally show that fitting the solvent-dependent s{sub 6} parameter to the solvation energy of a single non-polar molecule does not substantially increase these errors. Parametrization of this model for other solvents, therefore, requires minimal effort and is possible without extensive databases of experimental solvation free energies.
Continuum solvation models enable electronic structure calculations of systems in liquid environments, but because of the large number of empirical parameters, they are limited to the class of systems in their fit set (typically organic molecules). Here, we derive a solvation model with no empirical parameters for the dielectric response by taking the linear response limit of a classical density functional for molecular liquids. This model directly incorporates the nonlocal dielectric response of the liquid using an angular momentum expansion, and with a single fit parameter for dispersion contributions it predicts solvation energies of neutral molecules with a RMS error of 1.3 kcal/mol in water and 0.8 kcal/mol in chloroform and carbon tetrachloride. We show that this model is more accurate for strongly polar and charged systems than previous solvation models because of the parameter-free electric response, and demonstrate its suitability for ab initio solvation, including self-consistent solvation in quantum Monte Carlo calculations.
Relative solvating power, that is, the ratio of the coordination ratios between a solvent and the reference solvent, was used to probe the quantitative structure–activity relationship of electrolyte solvents and the lithium polysulfide (LiPS) dissolution in lithium–sulfur batteries. Internally referenced diffusion-ordered nuclear magnetic resonance spectroscopy (IR-DOSY) was used to determine the diffusion coefficient and coordination ratio, from which the relative solvating power can be easily measured. The higher the relative solvating power of an ethereal solvent, the more severe will be the LiPS dissolution and the lower the coulombic efficiency of the lithium–sulfur battery. A linear relationship was established between the logarithm of relative solvating power of a solvent and the degree of LiPS dissolution, rendering relative solvating power an important parameter in choosing the electrolyte solvent for lithium–sulfur batteries.
Kinetic models based on Fermi's Golden Rule are commonly employed to understand photoinduced electron transfer dynamics at molecule-semiconductor interfaces. Implicit in such second-order perturbative descriptions is the assumption that nuclear relaxation of the photoexcited electron donor is fast compared to electron injection into the semiconductor. This approximation breaks down in systems where electron transfer transitions occur on 100-fs time scale. Here, we present a fourth-order perturbative model that captures the interplay between time-coincident electron transfer and nuclear relaxation processes initiated by light absorption. The model consists of a fairly small number of parameters, which can be derived from standard spectroscopic measurements (e.g., linear absorbance, fluorescence) and/or first-principles electronic structure calculations. Insights provided by the model are illustrated for a two-level donor molecule coupled to both (i) a single acceptor level and (ii) a density of states (DOS) calculated for TiO{sub 2} using a first-principles electronic structure theory. These numerical calculations show that second-order kinetic theories fail to capture basic physical effects when the DOS exhibits narrow maxima near the energy of the molecular excited state. Overall, we conclude that the present fourth-order rate formula constitutes a rigorous and intuitive framework for understanding photoinduced electron transfer dynamics that occur on the 100-fs time scale.
Electrolyte is very critical to the performance of the high-voltage lithium (Li) metal battery (LMB), which is one of the most attractive candidates for the next generation high-density energy storage systems. Electrolyte formulation and structure determine the physical properties of the electrolytes and their interfacial chemistries on the electrode surfaces. Localized high-concentration electrolytes (LHCEs) outperform state-of-the-art carbonate electrolytes in many aspects in LMBs due to their unique solvation structure. Types of fluorinated co-solvents used in LHCEs are investigated here in searching for the most suitable diluent for high concentration electrolytes (HCEs). Non-solvating solvents (including fluorinated ethers, fluorinated borate and fluorinated orthoformate) added in HCEs enable the formation of LHCEs with high-concentration solvation structures. However, low solvating fluorinated carbonate will coordinate with Li+ ions and form a second solvation shell or a pseudo-LHCE which diminishes the benefits of LHCE. In addition, it is evident that the diluent has significant influence on the electrode/electrolyte interphases (EEIs) beyond retaining the high-concentration solvation structures. Diluent molecules surrounding the high-concentration clusters could accelerate or decelerate the anion decomposition through co-participation of diluent decomposition in the EEI formation. The varied interphase features lead to significantly different battery-performance. This study points out the importance of diluents and their synergetic effects with the conductive salt and the solvating solvent in designing LHCEs. These systematic comparisons and fundamental insights into LHCEs using different types of fluorinated solvents can guide further development of advanced electrolytes for high-voltage LMBs.
Li-S batteries are a promising next-generation battery technology. Due to the formation of soluble polysulfides during cell operation, the electrolyte composition of the cell plays an active role in directing the formation and speciation of the soluble lithium polysulfides. Recently, new classes of electrolytes termed "solvates" that contain stoichiometric quantities of salt and solvent and form a liquid at room temperature have been explored due to their sparingly solvating properties with respect to polysulfides. The viscosity of the solvate electrolytes is understandably high limiting their viability; however, hydrofluoroether cosolvents, thought to be inert to the solvate structure itself, can be introduced to reduce viscosity and enhance diffusion. Nazar and co-workers previously reported that addition of 1,1,2,2-tetrafluoroethyl 2,2,3,3-tetrafluoropropyl ether (TTE) to the LiTFSI in acetonitrile solvate, (MeCN) 2 -LiTFSI, results in enhanced capacity retention compared to the neat solvate. Here, we evaluate the effect of TTE addition on both the electrochemical behavior of the Li-S cell and the solvation structure of the (MeCN) 2 -LiTFSI electrolyte. Contrary to previous suggestions, Raman and NMR spectroscopy coupled with ab initio molecular dynamics simulations show that TTE coordinates to Li + at the expense of MeCN coordination, thereby producing a higher content of free MeCN, a good polysulfide solvent, in the electrolyte. Furthermore, the electrolytes containing a higher free MeCN content facilitate faster polysulfide formation kinetics during the electrochemical reduction of S in a Li-S cell likely as a result of the solvation power of the free MeCN.
Abstract We demonstrate that the solvation‐layer interface condition (SLIC) continuum dielectric model for molecular electrostatics, combined with a simple solvent‐accessible‐surface‐area (SASA)‐proportional model for nonpolar solvent effects, accurately predicts solvation entropies of neutral and charged small molecules. The SLIC/SASA model has only seven fitting parameters in total and achieves this accuracy using a training set with only 20 compounds. Despite this simplicity, solvation free energies and entropies are nearly as accurate as those predicted by the more sophisticated Langevin dipoles solvation model. Surprisingly, the model automatically reproduces the negligible contribution of electrostatics to the solvation of hydrophobic compounds. Opportunities for improvement include nonpolar solvation, anion solvation entropies, and heat capacities. More molecular realism may be needed for these quantities. To enable a future, explicit‐solvent‐based assessment of the SLIC/SASA implicit‐solvent model, we predict solvation entropies for the Mobley test set, which are available as Supporting Information.
The correlation between Li + solvation and interphasial chemistry on anodes firmly established in Li-ion batteries, the effect of cation–solvent interaction has gone beyond bulk thermodynamic and transport properties and become an essential element that determines the reversibility of electrochemistry and kinetics of Li-ion intercalation chemistries. Now, most studies are dedicated to the solvation of Li + , and the solvation of anions in carbonate-based electrolytes and its possible effect on the electrochemical stability of such electrolytes remains little understood. Moreover, as a mirror effort to prior Li + solvation studies, this work focuses on the interactions between carbonate-based solvents and two anions (hexafluorophosphate, PF 6– , and tetrafluoroborate, BF 4– ) that are most frequently used in Li-ion batteries. The possible correlation between such interaction and the interphasial chemistry on cathode surface is also explored.
Proton solvation properties and transport mechanisms have been studied in hydrated Nafion using the self-consistent multistate empirical valence bond (SCI-MS-EVB) method that includes the effects excess proton charge defect delocalization and Grotthuss proton hopping. It was found that sulfonate groups influence excess proton solvation, as well as the proton hydration structure, by stabilizing a more Zundel-like (H 5 O 2 + ) structure in their first solvation shells. Hydrate proton-related hydrogen bond networks were observed to be more stable than networks with water alone. Diffusion rates, Arrhenius activation energies, and transport pathways were calculated and analyzed to characterize the nature of the proton transport. Diffusion rate analysis suggests that a proton-hopping mechanism dominates the proton transport for the studied water loading levels and that there is a clear degree of anticorrelation with the vehicular transport. The activation energy drops quickly with an increasing water content when the water loading level is smaller than ~10 H 2 O/SO 3 – , which is consistent with experimental observations. The sulfonate groups were also found to affect the proton hopping directions. The temperature and water content effects on the proton transport pathways were also investigated.
Electrolyte design has long followed a solvation-first paradigm that prioritizes solvents capable of maximizing salt dissociation and ionic conductivity, while treating interfacial degradation and rheological limitations as secondary constraints. Although this approach has enabled significant progress in lithium-ion and sodium-ion batteries, it inherently favors solvent-dominated solvation structures that destabilize reactive metal interfaces. Weakly solvating electrolytes (WSEs) offer a fundamentally different strategy. By using solvents with intrinsically low donor strength and minimal electrostatic affinity for cations, WSEs suppress cation–solvent coordination and promote anion-rich solvation shells without relying on salt superconcentration. This shift lowers desolvation barriers, redirects interfacial decomposition pathways, and supports the formation of inorganic-rich, stable interphases. In this review, we discuss the molecular and solvation criteria that distinguish WSEs from conventional and concentrated or locally concentrated electrolyte systems, examine their implementation across different chemistries, and identify unresolved design challenges. WSEs are presented not as a niche formulation, but as a solvation-centric framework for rethinking electrolyte function in metal battery technologies.
We evaluate hydrofluoroether (HFE) cosolvents with varying degrees of fluorination in the acetonitrile-based solvate electrolyte to determine the effect of the HFE structure on the electrochemical performance of the Li-S battery. Solvates or sparingly solvating electrolytes are an interesting electrolyte choice for the Li-S battery due to their low polysulfide solubility. The solvate electrolyte with a stoichiometric ratio of LiTFSI salt in acetonitrile, (MeCN)(2)-LiTFSI, exhibits limited polysulfide solubility due to the high concentration of LiTFSI. We demonstrate that the addition of highly fluorinated HFEs to the solvate yields better capacity retention compared to that of less fluorinated HFE cosolvents. Raman and NMR spectroscopy coupled with ab initio molecular dynamics simulations show that HFEs exhibiting a higher degree of fluorination coordinate to Li+ at the expense of MeCN coordination, resulting in higher free MeCN content in solution. However, the polysulfide solubility remains low, and no crossover of polysulfides from the S cathode to the Li anode is observed.
Functional electrolyte is the key to stabilizing the highly reductive lithium (Li) metal anode (LMA) and high voltage cathode for long life, high energy-density rechargeable Li metal batteries (LMBs). However, fundamental knowledge of the interaction principles between reactive electrodes and electrolytes is still limited. Recently localized high-concentration electrolytes (LHCEs) are emerging as promising electrolyte design strategies for LMBs. They can also serve as an ideal platform for understanding the reactivity characteristics of the inner solvation sheath on electrode surfaces due to their unique solvation structures. Here, we study the effects of a series of LHCEs with model electrolyte solvents (carbonate, sulfone, phosphate and ether) in high voltage LMBs. Varied electrode stabilities exhibited in different LHCEs indicate the intricate synergies between the salt and the solvent on electrode surfaces. Experimental and theoretical analyses reveal an intriguing general rule that the strong interactions between the salt and the solvent in the inner solvation sheath promote their intermolecular proton/charge transfer reactions, which dictates the properties of the electrode/electrolyte interphases and thus the battery performances.
We report molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. The SLIC/CDC multiscale model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). The SLIC/CDC model predicts Gibbs energies of solvation for a database of 500 solutes in water with average accuracy better than 1 kcal/mol both for experimental measurements and for explicit-solvent molecular dynamics simulations. The separate SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent accessible surface area (SASA) nonpolar energy mode. The SLIC/SASA model predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems and better than 1.6 kcal/mol average accuracy in ionic liquids. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.