Ion aggregation in high salt solutions: Ion network versus ion cluster

Water-Ion Interaction Determines the Mobility of Ions in Highly Concentrated Aqueous Electrolytes

Solvation engineering plays a critical role in tailoring the performance of batteries, particularly through the use of highly concentrated electrolytes, which offer heterogeneous solvation structures of mobile ions with distinct electrochemical properties. In this study, we employed spectroscopic techniques and molecular dynamics simulations to investigate mixed-cation (Li+/K+) acetate aqueous electrolytes. Our research unravels the pivotal role of water in facilitating ion transport within a highly viscous medium. Notably, Li+ cations primarily form ion aggregates, predominantly interacting with acetate anions, while K+ cations emerge as the principal charge carriers, which is attributed to their strong interaction with water molecules. Intriguingly, even at a concentration as high as 40 m, a substantial amount of water molecules persistently engages in hydrogen bonding with one another, creating mobile regions rich in K+ ions. Our observations of a redshift of the OH stretching band of water suggest that the strength of the hydrogen bond alone cannot account for the expansion of the electrochemical stability window. These findings offer valuable insights into the cation transfer mechanism, shedding light on the contribution of water-bound cations to both the ion conductivity and the electrochemical stability window of aqueous electrolytes for rechargeable batteries. Our comprehensive molecular-level understanding of the interplay between cations and water provides a foundation for future advances in solvation engineering, leading to the development of high-performance batteries with improved energy storage and safety profiles.

From Networked to Isolated: Observing Water Hydrogen Bonds in Concentrated Electrolytes with Two-Dimensional Infrared Spectroscopy

Superconcentrated electrolytes have emerged as a promising class of materials for energy storage devices, with evidence that high voltage performance is possible even with water as the solvent. Here, we study the changes in the water hydrogen bonding network induced by the dissolution of lithium bis(trifluoromethane sulfonyl)imide (LiTFSI) in concentrations ranging from the dilute to the superconcentrated regimes. Using time-resolved two-dimensional infrared spectroscopy, we observe the progressive disruption of the water-water hydrogen bond network and the appearance of isolated water molecules interacting only with ions, which can be identified and spectroscopically isolated through the intermolecular cross-peaks between the water and the TFSI^- ions. Analyzing the vibrational relaxation of excitations of the H₂O stretching mode, we observe a transition in the dominant relaxation path as the bulk-like water vanishes and is replaced by ion-solvation water with the rapid single-step relaxation of delocalized stretching vibrations into the low frequency modes being replaced by multistep relaxation through the intramolecular H₂O bend and into the TFSI^- high frequency modes prior to relaxing to the low frequency structural degrees of freedom. These results definitively demonstrate the absence of vibrationally bulk-like water in the presence of high concentrations of LiTFSI and especially in the superconcentrated regime, while additionally revealing aspects of the water hydrogen bond network that have been difficult to discern from the vibrational spectroscopy of the neat liquid.

Applying Classical, Ab Initio, and Machine-Learning Molecular Dynamics Simulations to the Liquid Electrolyte for Rechargeable Batteries

Rechargeable batteries have become indispensable implements in our daily life and are considered a promising technology to construct sustainable energy systems in the future. The liquid electrolyte is one of the most important parts of a battery and is extremely critical in stabilizing the electrode-electrolyte interfaces and constructing safe and long-life-span batteries. Tremendous efforts have been devoted to developing new electrolyte solvents, salts, additives, and recipes, where molecular dynamics (MD) simulations play an increasingly important role in exploring electrolyte structures, physicochemical properties such as ionic conductivity, and interfacial reaction mechanisms. This review affords an overview of applying MD simulations in the study of liquid electrolytes for rechargeable batteries. First, the fundamentals and recent theoretical progress in three-class MD simulations are summarized, including classical, ab initio, and machine-learning MD simulations (section 2). Next, the application of MD simulations to the exploration of liquid electrolytes, including probing bulk and interfacial structures (section 3), deriving macroscopic properties such as ionic conductivity and dielectric constant of electrolytes (section 4), and revealing the electrode-electrolyte interfacial reaction mechanisms (section 5), are sequentially presented. Finally, a general conclusion and an insightful perspective on current challenges and future directions in applying MD simulations to liquid electrolytes are provided. Machine-learning technologies are highlighted to figure out these challenging issues facing MD simulations and electrolyte research and promote the rational design of advanced electrolytes for next-generation rechargeable batteries.

Towards understanding specific ion effects in aqueous media using thermodiffusion

Specific ion effects play an important role in scientific and technological processes. According to Hofmeister, the influence on the hydrogen bond network depends on the ion and leads to a specific order of the ions. Also thermodiffusion the mass transport caused by a temperature gradient is very sensitive to changes of the hydrogen bond network leading to a ranking according to hydrophilicity of the salt. Hence, we investigate various salt solutions in order to compare with the Hofmeister concept. We have studied three different sodium salts in water as a function of temperature (25-45[Formula: see text]C) and concentration (0.5-5 mol kg[Formula: see text]) using Thermal Diffusion Forced Rayleigh Scattering (TDFRS). The three anions studied, carbonate, acetate and thiocyanate, span the entire range of the Hofmeister series from hydrophilic to hydrophobic. We compare the results with the recent measurements of the corresponding potassium salts to see to what extent the cation changes the thermodiffusion of the salt.

Effects of molecular shape on alcohol aggregation and water hydrogen bond network behavior in butanol isomer solutions

Despite butanol isomers such as n-butanol, sec-butanol, isobutanol and tert-butanol having the same chemical formula, their liquid-liquid phase diagrams are distinct. That is, tert-butanol is miscible in water at all concentrations, while the other three butanol isomers are partially miscible under ambient conditions. The molecular shape of tert-butanol is close to globular and differs from the other three butanol molecules with a relatively long carbon chain. By performing molecular dynamics simulations and graph theoretical analysis of the four water-butanol isomer mixtures at varying concentrations, we show how distinct butanol aggregates are formed which depend upon the molecular shape and affect the water H-bond network structure and phase diagram in the binary liquid. The three butanol isomers of n-butanol, sec-butanol and isobutanol at concentrated solutions form chain-like alcohol aggregates, but tert-butanol forms small aggregates due to the distinct packing behavior caused by its globular molecular shape. By employing the graph theoretical analysis such as the degree distribution and the eigenvalue spectrum from the adjacency matrix in the graphical representation of the alcohol H-bond network, we show that the tert-butanol aggregates have a different morphological structure from that of the other three butanol isomers in aqueous solution. The graph theoretically distinct butanol aggregates are categorized into two groups, water-compatible and water-incompatible, depending upon the interaction between the alcohol and water molecules. Based upon our observations, we propose that the water-incompatible networks of n-butanol, sec-butanol and isobutanol aggregates do not change the water structure significantly, forming two separate liquid phases that are alcohol-rich and water-rich. However, the water-compatible network of tert-butanol aggregates has a considerable interaction with the water molecules and causes significant disruption of the water H-bond network, forming a homogeneous solution. Understanding the alcohol aggregation behavior and water structure in butanol-water mixtures provides a critical clue in appreciating fundamental issues such as miscibility and phase separation in aqueous solution systems.

Molecular heterogeneity in aqueous cosolvent systems

Aqueous cosolvent systems (ACoSs) are mixtures of small polar molecules such as amides, alcohols, dimethyl sulfoxide, or ions in water. These liquids have been the focus of fundamental studies due to their complex intermolecular interactions as well as their broad applications in chemistry, medicine, and materials science. ACoSs are fully miscible at the macroscopic level but exhibit nanometer-scale spatial heterogeneity. ACoSs have recently received renewed attention within the chemical physics community as model systems to explore the relationship between intermolecular interactions and microscopic liquid-liquid phase separation. In this perspective, we provide an overview of ACoS spatial segregation, dynamic heterogeneity, and multiscale relaxation dynamics. We describe emerging approaches to characterize liquid microstructure, H-bond networks, and dynamics using modern experimental tools combined with molecular dynamics simulations and network-based analysis techniques.

Insights into the hydrogen bond network topology of phosphoric acid and water systems

Phosphoric acid and its mixtures with water are some of the best proton conducting materials known to science. Although the proton conductivity in pure phosphoric acid decreases upon external doping with excess H+ or OH-, the addition of water improves it substantially. A number of experimental and theoretical studies indicate that these systems form a very special case of hydrogen bond networks which not only facilitate fast proton transport but also show a number of other interesting properties such as glass forming ability. In this work, we present the molecular dynamics simulation results of the H3PO4-H2O system over the entire concentration range. The hydrogen bond networks were analyzed in terms of conventional microscopic as well as topological properties based on graph and network theory. The results show that the hydrogen bond network of H3PO4 is fundamentally different from that of H2O. On average, each phosphoric acid molecule tends to form more and stronger hydrogen bonds than water which leads to a much more connected and clustered network showing small-world properties which are absent in pure water. Moreover, these hydrogen bond network properties persist in the H3PO4-H2O mixtures as well, even at relatively high water contents. Finally, many of the physical properties such as molecular diffusion coefficients seem to be also intimately related to the network topological properties and follow similar trends with respect to system content. These results strongly indicate that many important properties such as proton transport in phosphoric acid and its aqueous systems are fundamentally related to their hydrogen bond network topology and might hold the key for their ultimate molecular understanding.

Hydration breaking and chemical ordering in a levitated NaCl solution droplet beyond the metastable zone width limit: evidence for the early stage of two-step nucleation

For over two decades, NaCl nucleation from a supersaturated aqueous solution has been predicted to occur via a two-step nucleation (TSN) mechanism, i.e., two sequential events, the formation of locally dense liquid regions followed by structural ordering. However, the formation of dense liquid regions in the very early stage of TSN has never been experimentally observed. By using a state-of-the-art technique, a combination of electrostatic levitation (ESL) and in situ synchrotron X-ray and Raman scatterings, we find experimental evidence that indicates the formation of dense liquid regions in NaCl bulk solution at an unprecedentedly high level of supersaturation (S = 2.31). As supersaturation increases, evolution of ion clusters leads to chemical ordering, but no topological ordering, which is a precursor for forming the dense disordered regions of ion clusters at the early stage of TSN. Moreover, as the ion clusters proceed to evolve under highly supersaturated conditions, we observe the breakage of the water hydration structure indicating the stability limit of the dense liquid regions, and thus leading to nucleation. The evolution of solute clusters and breakage of hydration in highly supersaturated NaCl bulk solution will provide new insights into the detailed mechanism of TSN for many other aqueous solutions.

Specific Host-Guest Interactions in the Crown Ether Complexes with K+ and NH4+ Revealed from the Vibrational Relaxation Dynamics of the Counteranion

The specific host-guest interactions in the corresponding complexes of K⁺ and NH₄⁺ with typical crown ethers were investigated by using FTIR and ultrafast IR spectroscopies. The counteranions, i.e., SCN^-, were employed as a local vibrational probe to report the structural dynamics of the complexation. It was found that the vibrational relaxation dynamics of the SCN^- was strongly affected by the cations confined in the cavities of the crown ethers. The time constant of the vibrational population decay of SCN^- in the complex of NH₄⁺ with the 18-crown-6 was determined to be 6 ± 2 ps, which is ∼30 times faster than that in the complex of K⁺ with the crown ethers. Control experiments showed that the vibrational population decay of SCN^- depended on the size of the cavities of the crown ethers. A theoretical calculation further indicated that the nitrogen atom of SCN^- showed preferential coordination to the K⁺ ions hosted by the crown ethers, while the NH₄⁺ can form hydrogen bonds with the oxygen atoms in the studied crown ethers. The geometric constraints formed in the complex of crown ethers can cause a specific interaction between the NH₄⁺ and SCN^-, which can facilitate the intermolecular vibrational energy redistribution of the SCN^-.

Persistent homology analysis of ion aggregations and hydrogen-bonding networks

Despite the great advancement of experimental tools and theoretical models, a quantitative characterization of the microscopic structures of ion aggregates and their associated water hydrogen-bonding networks still remains a challenging problem. In this paper, a newly-invented mathematical method called persistent homology is introduced, for the first time, to quantitatively analyze the intrinsic topological properties of ion aggregation systems and hydrogen-bonding networks. The two most distinguishable properties of persistent homology analysis of assembly systems are as follows. First, it does not require a predefined bond length to construct the ion or hydrogen-bonding network. Persistent homology results are determined by the morphological structure of the data only. Second, it can directly measure the size of circles or holes in ion aggregates and hydrogen-bonding networks. To validate our model, we consider two well-studied systems, i.e., NaCl and KSCN solutions, generated from molecular dynamics simulations. They are believed to represent two morphological types of aggregation, i.e., local clusters and extended ion networks. It has been found that the two aggregation types have distinguishable topological features and can be characterized by our topological model very well. Further, we construct two types of networks, i.e., O-networks and H2O-networks, for analyzing the topological properties of hydrogen-bonding networks. It is found that for both models, KSCN systems demonstrate much more dramatic variations in their local circle structures with a concentration increase. A consistent increase of large-sized local circle structures is observed and the sizes of these circles become more and more diverse. In contrast, NaCl systems show no obvious increase of large-sized circles. Instead a consistent decline of the average size of the circle structures is observed and the sizes of these circles become more and more uniform with a concentration increase. As far as we know, these unique intrinsic topological features in ion aggregation systems have never been pointed out before. More importantly, our models can be directly used to quantitatively analyze the intrinsic topological invariants, including circles, loops, holes, and cavities, of any network-like structures, such as nanomaterials, colloidal systems, biomolecular assemblies, among others. These topological invariants cannot be described by traditional graph and network models.

The effect of Hofmeister anions on water structure at protein surfaces

To understand the effects of specific ions on protein-water interactions and the thermodynamic stability of proteins in salt solutions, we use a molecular dynamics (MD) simulation to examine the water structure, orientational distribution, and dynamics near the surface of ubiquitin. In particular, we consider NaCl, NaBF₄, NaSCN, and NaClO₄ salt solutions containing ubiquitin, where the anions of the latter three salts are well-known chaotropic ions in the Hofmiester anion series. The number of hydrogen bonds (H-bonds) per water molecule is found to decrease significantly at the ubiquitin-water interface, indicating a significant disruption of the water H-bonding network. The distribution of the water H-bond numbers near the protein surface is modulated by dissolved ions, and the extent of the ion effect on the H-bonding network structure follows the order of the Hofmeister anion series, while there are no specific ion effects on water properties at distances larger than 5 Å from the protein surface. From detailed analyses of the surface area, volume, and root-mean-square deviation (RMSD) of ubiquitin, we show that changes in the properties of the protein could originate from the disruption of the water H-bond network induced by ions with a higher affinity for the protein surface instead of direct protein residue-ion interactions. An interesting observation made here is that the orientational distribution of water molecules at the protein-water interface is close to random, but there is a slight preference for interfacial water molecules with a straddle structure within 2.5 Å of the protein surface, where one of the two OH groups points away from the protein surface and the other points toward the surface. In addition, comparing the MD simulation results for ubiquitin solutions with dissolved NaSCN and KSCN, we show that Na⁺ affects the water H-bonding structure at the protein surface more than K⁺. It is clear that the H-bonding network structure of water more than one water layer away from the protein surface is not distinguishably different from that of neat water. We thus anticipate that the present work will provide insights into the scale of specific ion effects on the H-bonding structure and orientational distribution of water in the vicinity of protein surfaces in aqueous solutions.

Dynamic heterogeneity in aqueous ionic solutions

It is well known that supercooled liquids have heterogeneous dynamics, but it is still unclear whether dynamic heterogeneity also exists in aqueous ionic solutions at room or even higher temperatures.

Molecular mechanism of water reorientational slowing down in concentrated ionic solutions

Water dynamics in concentrated ionic solutions plays an important role in a number of material and energy conversion processes such as the charge transfer at the electrolyte-electrode interface in aqueous rechargeable ion batteries. One long-standing puzzle is that all electrolytes, regardless of their "structure-making/breaking" nature, make water rotate slower at high concentrations. To understand this effect, we present a theoretical simulation study of the reorientational motion of water molecules in different ionic solutions. Using an extended Ivanov model, water rotation is decomposed into contributions from large-amplitude angular jumps and a slower frame motion which was studied in a coarse-grained manner. Bearing a certain resemblance to water rotation near large biological molecules, the general deceleration is found to be largely due to the coupling of the slow, collective component of water rotation with the motion of large hydrated ion clusters ubiquitously existing in the concentrated ionic solutions. This finding is at variance with the intuitive expectation that the slowing down is caused by the change in fast, single-molecular water hydrogen bond switching adjacent to the ions.

Role of Solvent Water in the Temperature-Induced Self-Assembly of a Triblock Copolymer

Water-soluble triblock copolymers have received much attention in industrial applications and scientific fields. We here show that femtosecond mid-IR pump-probe spectroscopy is useful to study the role of water in the temperature-induced self-assembly of triblock copolymers. Our experimental results suggest two distinct subpopulations of water molecules: those that interact with other water molecules and those involved in the hydration of a triblock copolymer surface. We find that the vibrational dynamics of bulk-like water is not affected by either micellation or gelation of triblock copolymers. The increased population of water interacting with ether oxygen atoms of the copolymer during the unimer to micelle phase transition is important evidence for the entropic role of water in temperature-induced micelle formation at a low copolymer concentration. In contrast, at the critical gelation temperature and beyond, the population of surface-associated water molecules interacting with ether oxygen atoms decreases, which indicates important enthalpic control by water. The present study on the roles of water in the two different phase transitions of triblock copolymers sheds new light on the underlying mechanisms of temperature-induced self-aggregation behaviors of amphiphiles that are ubiquitous in nature.

Ultrafast fluxional exchange dynamics in electrolyte solvation sheath of lithium ion battery

Lithium cation is the charge carrier in lithium-ion battery. Electrolyte solution in lithium-ion battery is usually based on mixed solvents consisting of polar carbonates with different aliphatic chains. Despite various experimental evidences indicating that lithium ion forms a rigid and stable solvation sheath through electrostatic interactions with polar carbonates, both the lithium solvation structure and more importantly fluctuation dynamics and functional role of carbonate solvent molecules have not been fully elucidated yet with femtosecond vibrational spectroscopic methods. Here we investigate the ultrafast carbonate solvent exchange dynamics around lithium ions in electrolyte solutions with coherent two-dimensional infrared spectroscopy and find that the time constants of the formation and dissociation of lithium-ion···carbonate complex in solvation sheaths are on a picosecond timescale. We anticipate that such ultrafast microscopic fluxional processes in lithium-solvent complexes could provide an important clue to understanding macroscopic mobility of lithium cation in lithium-ion battery on a molecular level.

Simultaneous spectral and temporal analyses of kinetic energies in nonequilibrium systems: theory and application to vibrational relaxation of O-D stretch mode of HOD in water

A time series of kinetic energies (KE) from classical molecular dynamics (MD) simulation contains fundamental information on system dynamics. It can also be analyzed in the frequency domain through Fourier transformation (FT) of velocity correlation functions, providing energy content of different spectral regions. By limiting the FT time span, we have previously shown that spectral resolution of KE evolution is possible in the nonequilibrium situations [Jeon and Cho, J. Chem. Phys. 2011, 135, 214504]. In this paper, we refine the method by employing the concept of instantaneous power spectra, extending it to reflect an instantaneous time-correlation of velocities with those in the future as well as with those in the past, and present a new method to obtain the instantaneous spectral density of KE (iKESD). This approach enables the simultaneous spectral and temporal resolution of KE with unlimited time precision. We discuss the formal and novel properties of the new iKESD approaches and how to optimize computational methods and determine parameters for practical applications. The method is specifically applied to the nonequilibrium MD simulation of vibrational relaxation of the OD stretch mode in a hydrated HOD molecule by employing a hybrid quantum mechanical/molecular mechanical (QM/MM) potential. We directly compare the computational results with the OD band population relaxation time profiles extracted from the IR pump-probe measurements for 5% HOD in water. The calculated iKESD yields the OD bond relaxation time scale ∼30% larger than the experimental value, and this decay is largely frequency-independent if the classical anharmonicity is accounted for. From the integrated iKESD over intra- and intermolecular bands, the major energy transfer pathways were found to involve the HOD bending mode in the subps range, then the internal modes of the solvent until 5 ps after excitation, and eventually the solvent intermolecular modes. Also, strong hydrogen-bonding of HOD is found to significantly hinder the initial intramolecular energy transfer process.