Fast Structural Dynamics in Concentrated HCl Solutions: From Proton Hopping to the Bulk Viscosity

Concentrated acid solutions, particularly HCl, have been studied extensively to examine the proton hopping and infrared spectral signatures of hydronium ions. Much less attention has been given to the structural dynamics of concentrated HCl solutions. Here, we apply optical heterodyne detected-optical Kerr effect (OHD-OKE) measurements to examine HCl concentration-dependent dynamics from moderate (0.8 m) to very high (15.5 m) concentrations and compare the results to the dynamics of NaCl solutions, as Na+ is similar in size to the hydronium cation. Both HCl and NaCl OHD-OKE signals decay as triexponentials at all concentrations, in contrast to pure water, which decays as a biexponential. Two remarkable features of the HCl dynamics are the following: (1) the bulk viscosity is linearly related to the slowest decay constant, t3, and (2) the concentration-dependent proton hopping times, determined by ab initio MD simulations and 2D IR chemical exchange experiments, both obtained from the literature, fall on the same line as the slowest structural dynamics relaxation time, t3, within experimental error. The structural dynamics of hydronium/chloride/water clusters, with relaxation times t3, are responsible for the concentration dependence of microscopic property of proton hopping and the macroscopic bulk viscosity. The slowest time constant (t3), which does not have a counterpart in pure water, is 3 ps at 0.8 m and increases by a factor of ∼2 by 15.5 m. The two fastest HCl decay constants, t1 and t2, are similar to those of pure water and increase mildly with the concentration.

Comparative Investigation of the Microstructure of MgCl2 Aqueous Solutions Using Different X-ray Scattering Sources, Raman Spectroscopy, and Atomistic Simulations

Aqueous solutions of magnesium chloride (MgCl2(aq)) are often used to test advances in the theory of electrolyte solutions because they are considered an ideal strong 2:1 electrolyte. However, there is evidence that some ion association occurs in these solutions, even at low concentrations. Even a small ion-pairing constant can have a significant impact on the chemical speciation of ions, so it is important to determine whether ion pairing actually occurs. In this study, MgCl2(aq) with concentrations ranging from 1 to 35% was studied using three methods: X-ray scattering (XRS) with the Shanghai Synchrotron Radiation Facility (SSRF) and silver-anode laboratory sources, Raman spectroscopy, and molecular dynamics (MD) simulations with the COMPASS-II and Madrid force fields. XRS results were analyzed in the framework of PDF theory to obtain the reduced structure function F(Q) and the reduced pair distribution function G(r). The F(Q) values from synchrotron radiation and laboratory sources both showed that the tetrahedral hydrogen bonds in bulk water were destroyed with the increased MgCl2 concentration. The results of G(r) indicated that the main peaks centered at 2.05 and 2.80 Å can be ascribed to the interactions of Mg-O and O-O, respectively. The peak at 3.10 Å is attributed to the combined effect of O-O and Cl-O. By comparing the structural information on MgCl2 solution obtained from the two light sources, it was found that both SSRF and silver-anode laboratory sources can reflect the above-mentioned structural information on MgCl2 solution. The radial distribution function (RDF) obtained from MD simulations of MgCl2 solutions assigned the peaks at 2.0, 2.8, and 3.2 Å to the Mg-O, O-O, and Cl-O interatomic pairs, respectively. The decrease in the O-O coordination number confirms that the hydrogen-bonding network of water is disrupted by increasing MgCl2 observed by X-ray scattering. The proportion of Mg-Cl contact ion pairs gradually increases with MgCl2 concentration as does the coordination number. Raman spectroscopy results show that the bond type changes from double donor double acceptor (DDAA) to single donor-single acceptor (DA) with increasing concentration, providing explicit details of the hydrogen-bond evolution in the aqueous solution.

Water-in-Salt: Fast Dynamics, Structure, Thermodynamics, and Bulk Properties

We present concentration-dependent dynamics of highly concentrated LiBr solutions and LiCl temperature-dependent dynamics for two high concentrations and compare the results to those of prior LiCl concentration-dependent data. The dynamical data are obtained using ultrafast optical heterodyne-detected optical Kerr effect (OHD-OKE). The OHD-OKE decays are composed of two pairs of biexponentials, i.e., tetra-exponentials. The fastest decay (t1) is the same as pure water's at all concentrations within error, while the second component (t2) slows slightly with concentration. The slower components (t3 and t4), not present in pure water, slow substantially, and their contributions to the decays increase significantly with increasing concentration, similar to LiCl solutions. Simulations of LiCl solutions from the literature show that the slow components arise from large ion/water clusters, while the fast components are from ion/water structures that are not part of large clusters. Temperature-dependent studies (15-95 °C) of two high LiCl concentrations show that decreasing the temperature is equivalent to increasing the room temperature concentration. The LiBr and LiCl concentration dependences and the two LiCl concentrations' temperature dependences all have bulk viscosities that are linearly dependent on τcslow, the correlation time of the slow dynamics (weighted averages of t3 and t4). Remarkably, all four viscosity vs 1/τCslow plots fall on the same line. Application of transition state theory to the temperature-dependent data yields the activation enthalpies and entropies for the dynamics of the large ion/water clusters, which underpin the bulk viscosity.

Lifting Hofmeister's Curse: Impact of Cations on Diffusion, Hydrogen Bonding, and Clustering of Water

Water plays a role in the stability, reactivity, and dynamics of the solutes that it contains. The presence of ions alters this capacity by changing the dynamics and structure of water. However, our understanding of how and to what extent this occurs is still incomplete. Here, a study of the low-frequency Raman spectra of aqueous solutions of various cations by using optical Kerr-effect spectroscopy is presented. This technique allows for the measurement of the changes that ions cause in both the diffusive dynamics and the vibrations of the hydrogen-bond structure of water. It is found that when salts are added, some of the water molecules become part of the ion solvation layers, while the rest retain the same diffusional properties as those of pure water. The slowing of the dynamics of the water molecules in the solvation shell of each ion was found to depend on its charge density at infinite dilution conditions and on its position in the Hofmeister series at higher concentrations. It is also observed that all cations weaken the hydrogen-bond structure of the solution and that this weakening depends only on the size of the cation. Finally, evidence is found that ions tend to form amorphous aggregates, even at very dilute concentrations. This work provides a novel approach to water dynamics that can be used to better study the mechanisms of solute nucleation and crystallization, the structural stability of biomolecules, and the dynamic properties of complex solutions, such as water-in-salt electrolytes.

Structure and dynamics of the Li+ ion in water, methanol and acetonitrile solvents: ab initio molecular dynamics simulations

Fundamental understanding of the structure and dynamics of the Li+ ion in solution is of utmost importance in different fields of science and technology, especially in the field of ion batteries. In view of this, ab initio molecular dynamics (AIMD) simulations of the LiCl salt in water, methanol and acetonitrile were performed to elucidate structural parameters such as radial distribution function and coordination number, and dynamical properties like diffusion coefficient, limiting ion conductivity and hydrogen bond correlation function. In the present AIMD simulation, one LiCl in water is equivalent to 0.8 M, which is close to the concentration of the lithium salt used in the Li-ion battery. The first sphere of coordination number of the Li+ ion was reaffirmed to be 4. The radial distribution function for different pairs of atoms is seen to be in good agreement with the experimental results. The calculated potential of mean force indicates the stronger interaction of the Li+ ion with methanol over water followed by acetonitrile. The dynamical parameters convey quite high diffusion and limiting ionic conductivity of the Li+ ion in acetonitrile compared to that in water and methanol which has been attributed to the transport of the Li-Cl ion pair in a non-dissociated form in acetonitrile. The AIMD results were found to be in accordance with the experimental findings, i.e. the limiting ion conductivity was found to follow the order acetonitrile > methanol > water. This study shows the importance of atomistic level simulations in evaluating the structural and dynamical parameters and in implementing the results for predicting and synthesizing better next generation solvents for lithium ion batteries (LIBs).

All-Day Multicyclic Atmospheric Water Harvesting Enabled by Polyelectrolyte Hydrogel with Hybrid Desorption Mode

Sorption-based atmospheric water harvesting (AWH) is a promising approach for mitigating worldwide water scarcity. However, reliable water supply driven by sustainable energy regardless of diurnal variation and weather remains a long-standing challenge. To address this issue, a polyelectrolyte hydrogel sorbent with an optimal hybrid-desorption multicyclic-operation strategy is proposed, achieving all-day AWH and a significant increase in daily water production. The polyelectrolyte hydrogel possesses a large interior osmotic pressure of 659 atm, which refreshes sorption sites by continuously migrating the sorbed water within its interior, and thus enhancing sorption kinetics. The charged polymeric chains coordinate with hygroscopic salt ions, anchoring the salts and preventing agglomeration and leakage, thereby enhancing cyclic stability. The hybrid desorption mode, which couples solar energy and simulated waste heat, introduces a uniform and adjustable sorbent temperature for achieving all-day ultrafast water release. With rapid sorption-desorption kinetics, an optimization model suggests that eight moisture capture-release cycles are capable of achieving high water yield of 2410 mLwater kgsorbent -1 day-1 , up to 3.5 times that of single-cyclic non-hybrid modes. The polyelectrolyte hydrogel sorbent and the coupling with sustainable energy driven desorption mode pave the way for the next-generation AWH systems, significantly bringing freshwater on a multi-kilogram scale closer.

Chain-like Structures Facilitate Li+ Transport in Concentrated Aqueous Electrolytes: Insights from Ultrafast Infrared Spectroscopy and Molecular Dynamics Simulations

Highly concentrated aqueous electrolytes have attracted attention due to their unique applications in lithium ion batteries (LIBs). However, the solvation structure and transport mechanism of Li+ cations at concentrated concentrations remain largely unexplored. To address this gap in knowledge, we employ ultrafast infrared spectroscopy and molecular dynamics (MD) simulations to reveal the dynamic and spatial structural heterogeneity in aqueous lithium chloride (LiCl) solutions. The coupling between the reorientation dynamics of the extrinsic probe and the macroscopic viscosity in aqueous LiCl solutions was analyzed using the Stokes-Einstein-Debye (SED) equations. MD simulations reveal that the Cl- and Li+ form chain-like structures through electrostatic interactions, supporting the vehicular migration of Li+ through the chain-like structure. The concentration dependent conductivity of the LiCl solution is well reproduced, where Li(H2O)2+ and Li(H2O)3+ are the dominant species that contribute to the conduction of Li+. This study is expected to establish correlations between ion pair structures and macroscopic properties.

Ion/Water Network Structural Dynamics in Highly Concentrated Lithium Chloride and Lithium Bromide Solutions Probed with Ultrafast Infrared Spectroscopy

The structural dynamics of highly concentrated LiCl and LiBr aqueous solutions were observed from 1-4 to 1-16 water molecules per ion pair using ultrafast polarization-selective pump-probe (PSPP) experiments on the OD stretch of dilute HOD. At these high salt concentrations, an extended ion/water network exists with complex structural dynamics. Population decays from PSPP experiments highlight two distinct water components. From the frequency-dependent amplitudes of the decays, the spectra of hydroxyls bound to halides and to water oxygens are obtained, which are not observable in the FT-IR spectra. PSPP experiments also measure frequency-dependent water orientational relaxation. At short times, wobbling dynamics within a restricted angular cone occurs. At high concentrations, the cone angles are dependent on frequency (hydrogen bond strength), but at higher water concentrations (>10 waters per ion pair), there is no frequency dependence. The average cone angle increases as the ion concentration decreases. The slow time constant for complete HOD orientational relaxation is independent of concentration but slower in LiCl than in LiBr. Comparison to structural MD simulations of LiCl from the literature indicates that the loss of the cone angle wavelength dependence and the increase in the cone angles as the concentration decreases occur as the prevalence of large ion/water clusters gives way to contact ion pairs.

Dynamics of Concentrated Aqueous Lithium Chloride Solutions Investigated with Optical Kerr Effect Experiments

We report the dynamics of concentrated lithium chloride aqueous solutions over a range of moderate to high concentrations. Concentrations (1-29 to 1-3.3 LiCl-water) were studied in which, at the highest concentrations, there are far too few water molecules to solvate the ions. The measurements were made with optically heterodyne-detected optical Kerr effect experiments, a non-resonant technique able to observe dynamics over a wide range of time scales and signal amplitudes. While the pure water decay is a biexponential, the LiCl-water decays are tetra-exponentials at all concentrations. The faster two decays arise from water dynamics, while the slower two decays reflect the dynamics of the ion-water network. The fastest decay (t₁) is the same as pure water at all concentrations. The second decay (t₂) is also the same as that of pure water at the lower concentrations, and then, it slows with increasing concentration. The slower dynamics (t₃ and t₄), which do not have counterparts in pure water, arise from ion-water complexes and, at the highest concentrations, an extended ion-water network. Comparisons are made between the concentration dependence of the observed dynamics and simulations of structural changes from the literature, which enable the assignment of dynamics to specific ion-water structures. The concentration dependences of the bulk viscosity and the ion-water network dynamics are directly correlated. The correlation provides an atomistic-level understanding of the viscosity.

Water Dynamics and Structure of Highly Concentrated LiCl Solutions Investigated Using Ultrafast Infrared Spectroscopy

In highly concentrated salt solutions, the water hydrogen bond (H-bond) network is completely disrupted by the presence of ions. Water is forced to restructure as dictated by the water-ion and ion-ion interactions. Using ultrafast polarization-selective pump-probe (PSPP) spectroscopy measurements of the OD stretch of dilute HOD, we demonstrate that the limited water-water H-bonding present in concentrated lithium chloride solutions (up to four waters per ion pair) is, on average, stronger than that occurring in bulk water. Furthermore, information on the orientational dynamics and the angular restriction of water H-bonded to both water oxygens and chloride anions was obtained through analysis of the frequency-dependent anisotropy decays. It was found that, when the salt concentration increased, the water showed increasing restriction and slowing at frequencies correlated with strong H-bonding. The angular restriction of the water molecules and strengthening of water-water H-bonds are due to the formation of a water-ion network not present in bulk water and dilute salt solutions. The structural evolution of the ionic medium was also observed through spectral diffusion of the OD stretch using 2D IR spectroscopy. Compared to bulk water, there is significant slowing of the biexponential spectral diffusion dynamics. The slowest component of the spectral diffusion (13 ps) is virtually identical to the time for complete reorientation of HOD measured with the PSPP experiments. This result suggests that the slowest component of the spectral diffusion reflects rearrangement of water molecules in the water-ion network.

A Numerical Study on Chloride Diffusion in Cracked Concrete

The cracks in concrete are a fast transport path for chlorides and influence the service life of concrete structures in chloride environments. This study aimed to reveal the effect of crack geometry on chloride diffusion in cracked concrete. The chloride diffusion process in cracked concrete was simulated with the finite difference method by solving Fick’s law. The results showed that the apparent chloride diffusivity was lower in more tortuous cracks, and the cracks with more narrow points also showed lower apparent chloride diffusivity. For tortuous cracks, a higher crack width meant relatively more straight cracks, and consequently, higher apparent chloride diffusivity, while a lower crack width resulted in more tortuous cracks and lower apparent chloride diffusivity. The crack depth showed a more significant influence on the chloride penetration depth in cracked concrete than crack geometry did. Compared with rectangular and V-shaped cracks, the chloride diffusion process in cracked concrete with a tortuous crack was slower at the early immersion age. At the same crack depth, the crack geometry showed a marginal influence on the chloride penetration depth in cracked concrete during long-term immersion.

On the transferability of ion parameters to the TIP4P/2005 water model using molecular dynamics simulations

Countless molecular dynamics studies have relied on available ion and water force field parameters to model aqueous electrolyte solutions. The TIP4P/2005 model has proven itself to be among the best rigid water force fields, whereas many of the most successful ion parameters were optimized in combination with SPC/E, TIP3P, or TIP4P/Ew water. Many researchers have combined these ions with TIP4P/2005, hoping to leverage the strengths of both parameter sets. To assess if this widely used approach is justified and to provide a guide in selecting ion parameters, we investigated the transferability of various commonly used monovalent and multivalent ion parameters to the TIP4P/2005 water model. The transferability is evaluated in terms of ion hydration free energy, hydration radius, coordination number, and self-diffusion coefficient at infinite dilution. For selected ion parameters, we also investigated density, ion pairing, chemical potential, and mean ionic activity coefficients at finite concentrations. We found that not all ions are equally transferable to TIP4P/2005 without compromising their performance. In particular, ions optimized for TIP3P water were found to be poorly transferable to TIP4P/2005, whereas ions optimized for TIP4P/Ew water provided nearly perfect transferability. The latter ions also showed good overall agreement with experimental values. The one exception is that no combination of ion parameters and water model considered here was found to accurately reproduce experimental self-diffusion coefficients. Additionally, we found that cations optimized for SPC/E and TIP3P water displayed consistent underpredictions in the hydration free energy, whereas anions consistently overpredicted the hydration free energy.