Role of local order in anomalous ion diffusion: Interrogation through tetrahedral entropy of aqueous solvation shells

Characterization of the Coordination and Solvation Dynamics of Solvated Systems─Implications for the Analysis of Molecular Interactions in Solutions and Pure H2O

The characterization of solvation shells of atoms, ions, and molecules in solution is essential to relate solvation properties to chemical phenomena such as complex formation and reactivity. Different definitions of the first-shell coordination sphere from simulation data can lead to potentially conflicting data on the structural properties and associated ligand exchange dynamics. The definition of a solvation shell is typically based on a given threshold distance determined from the respective solute-solvent pair distribution function g(r) (i.e., GC). Alternatively, a nearest neighbor (NN) assignment based on geometric properties of the coordination complex without the need for a predetermined cutoff criterion, such as the relative angular distance (RAD) or the modified Voronoi (MV) tessellation, can be applied. In this study, the effect of different NN algorithms on the coordination number and ligand exchange dynamics evaluated for a series of monatomic ions in aqueous solution, carbon dioxide in aqueous and dichloromethane solutions, and pure liquid water has been investigated. In the case of the monatomic ions, the RAD approach is superior in achieving a well separated definition of the first solvation layer. In contrast, the MV algorithm provides a better separation of the NNs from a molecular point of view, leading to better results in the case of solvated CO2. When analyzing the coordination environment in pure water, the cutoff-based GC framework was found to be the most reliable approach. By comparison of the number of ligand exchange reactions and the associated mean ligand residence times (MRTs) with the properties of the coordination number autocorrelation functions, it is shown that although the average coordination numbers are sensitive to the different definitions of the first solvation shell, highly consistent estimates for the associated MRT of the solvated system are obtained in the majority of cases.

Understanding the Anomalous Diffusion of Water in Aqueous Electrolytes Using Machine Learned Potentials

The diffusivity of water in aqueous cesium iodide solutions is larger than that in neat liquid water and vice versa for sodium chloride solutions. Such peculiar ion-specific behavior, called anomalous diffusion, is not reproduced in typical force field based molecular dynamics (MD) simulations due to inadequate treatment of ion-water interactions. Herein, this hurdle is tackled by using machine learned atomic potentials (MLPs) trained on data from density functional theory calculations. MLP based atomistic MD simulations of aqueous salt solutions reproduce experimentally determined thermodynamic, structural, dynamical, and transport properties, including their varied trends in water diffusivities across salt concentration. This enables an examination of their intermolecular structure to unravel the microscopic underpinnings of the differences in their transport properties. While both ions in CsI solutions contribute to the faster diffusion of water molecules, the competition between the heavy retardation by Na ions and the slight acceleration by Cl ions in NaCl solutions reduces their water diffusivity.

Tailoring water structure with high-tetrahedral-entropy for antifreezing electrolytes and energy storage at -80 °C

One of unsolved puzzles about water lies in how ion-water interplay affects its freezing point. Here, we report the direct link between tetrahedral entropy and the freezing behavior of water in Zn²⁺-based electrolytes by analyzing experimental spectra and molecular simulation results. A higher tetrahedral entropy leads to lower freezing point, and the freezing temperature is directly related to the entropy value. By tailoring the entropy of water using different anions, we develop an ultralow temperature aqueous polyaniline| |Zn battery that exhibits a high capacity (74.17 mAh g^-1) at 1 A g^-1 and -80 °C with ~85% capacity retention after 1200 cycles due to the high electrolyte ionic conductivity (1.12 mS cm^-1). Moreover, an improved cycling life is achieved with ~100% capacity retention after 5000 cycles at -70 °C. The fabricated battery delivers appreciably enhanced performance in terms of frost resistance and stability. This work serves to provide guidance for the design of ultralow temperature aqueous batteries by precisely tuning the water structure within electrolytes.

Diffusions in Aqueous Solutions with Multivalent Cations and Especially in Cationic First Hydration Shell

Compared with univalent cationic diffusion, little is known about the diffusion behavior of multivalent cations let alone the diffusion of water in their first hydration shell. Here, we show that all published translational diffusion coefficients of multivalent cations and water measured at room temperature exhibit the same concentration dependence when plotted as a function of the mass fraction of free water or of hydrated solute. This behavior is held until only hydration and confined water remain in solutions, wherein their concentration dependences become cationic charge number-dependent. We also found that the iceberg model can well describe the diffusions of multivalent cations with decreasing water content until only hydration water is present. However, ¹H-pulsed-field-gradient nuclear magnetic resonance measurements confirmed that ¹H in the first hydration shell diffuses faster than Al³⁺ at room temperature and they have the same diffusion coefficient only with decreasing the temperature down to about 243 K. Therefore, iceberg model only equivalently describes the effect of strong ion-water interaction on multivalent cationic diffusion. These results will also help us reconceive our current understanding of the pathway for hydration water affecting the diffusion behavior of solutes with relatively weak solute-water interactions.

Rigid Cations Induce Enhancement of Microheterogeneity and Exhibit Anomalous Ion Diffusion in Water-Ethanol Mixtures

Because of the amphiphilic nature of ethanol in the aqueous solution, ions cause an interesting microheterogeneity where the water molecules and the hydroxy groups of ethanol preferentially solvate the ions, while the ethyl groups tend to occupy the intervening space. Using computer simulations, we study the dynamics of rigid monovalent cations (Li⁺, Na⁺, K⁺, and Cs⁺) in aqueous ethanol solutions with chloride as the counterion. We vary both the size of the ions and the composition of the mixture to explore size- and composition-dependent ion diffusion. The relative stability of enhanced microheterogeneous configurations makes ion diffusion slower than what would be surmised by using the bulk properties of the mixture, using the Stokes-Einstein relation. We study the structure through partial radial distribution functions and the stability through coordination number fluctuations. The ion diffusion coefficient exhibits sharp re-entrant behavior when plotted against viscosity varied by composition. Our studies reveal multiple anomalous features of ion motion in this mixture. We formulate a mode-coupling theory (MCT) that takes into account the interaction between different dynamical components; MCT can incorporate the effects of heterogeneous dynamics and nonlinearity in composition dependence that arise from the feedback between mutually dependent ion-solvent dynamics.

Molecular Mechanism of Acceleration and Retardation of Collective Orientation Relaxation of Water Molecules in Aqueous Solutions

The collective orientation relaxation (COR) of water molecules in aqueous solutions is faster or slower with an increase in the concentration of the solutions than that in pure water; for example, acceleration (deceleration) of the COR is observed in a solution of sodium chloride (tetramethylammonium chloride) with increasing concentration. However, the molecular mechanism of the solution and concentration dependence of the relaxation time of the COR has not yet been clarified. We theoretically investigate the concentration dependence of the COR of water molecules in solutions of tetramethylammonium chloride (TMACl), guanidinium chloride (GdmCl), and sodium chloride (NaCl). Based on the Mori-Zwanzig equation, we identify two opposing factors that determine the COR of water molecules in any aqueous solution: the correlation of dipole moments and the single-molecule orientation relaxation. We reveal the molecular mechanism of the concentration dependence of the relaxation time of the COR in the TMACl, GdmCl, and NaCl solutions in terms of these two factors.