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

Dynamic Ion Correlations and Ion-Pair Lifetimes in Aqueous Alkali Metal Chloride Electrolytes

Electrolytes are central to many technological applications, as well as life itself. The behavior and properties of electrolytes are often described in terms of ion pairs, whereby ions associate as either contact ion pairs (in which ions are "touching") solvent-separated ion pairs (in which ions' solvent shells overlap) or solvent-solvent-separated ion pairs (in which ions' solvent shells are distinct). However, this paradigm is generally restricted to statistically averaged descriptions of solution structure and ignores temporal behavior. Here we elucidate the time-resolved dynamics of these ion-ion interactions in aqueous metal chloride electrolytes using the partial van Hove correlation function, based on polarizable molecular dynamics simulations. Our results show that the existence and persistence of ion pairs in aqueous metal chloride electrolytes should not be assumed a priori, but in fact are ion specific features of the solution with lifetimes on subpicosecond time scales.

Unveiling the Structural and Dynamic Characteristics of Concentrated LiNO3 Aqueous Solutions through Ultrafast Infrared Spectroscopy and Molecular Dynamics Simulations

Highly concentrated aqueous electrolytes have attracted a significant amount of attention for their potential applications in lithium-ion batteries. Nevertheless, a comprehensive understanding of the Li+ solvation structure and its migration within electrolyte solutions remains elusive. This study employs linear vibrational spectroscopy, ultrafast infrared spectroscopy, and molecular dynamics (MD) simulations to elucidate the structural dynamics in LiNO3 solutions by using intrinsic and extrinsic vibrational probes. The N-O stretching vibrations of NO3- exhibit a distinct spectral splitting, attributed to its asymmetric interaction with the surrounding solvation structure. Analysis of the vibrational relaxation dynamics of intrinsic and extrinsic probes, in combination with MD simulations, reveals cage-like networks formed through electrostatic interactions between Li+ and NO3-. This microscopic heterogeneity is reflected in the intertwined arrangement of ions and water molecules. Furthermore, both vehicular transport and structural diffusion assisted by solvent rearrangement for Li+ were analyzed, which are closely linked with the bulk concentration.

Accelerated Bayesian Inference for Molecular Simulations using Local Gaussian Process Surrogate Models

While Bayesian inference is the gold standard for uncertainty quantification and propagation, its use within physical chemistry encounters formidable computational barriers. These bottlenecks are magnified for modeling data with many independent variables, such as X-ray/neutron scattering patterns and electromagnetic spectra. To address this challenge, we employ local Gaussian process (LGP) surrogate models to accelerate Bayesian optimization over these complex thermophysical properties. The time-complexity of the LGPs scales linearly in the number of independent variables, in stark contrast to the computationally expensive cubic scaling of conventional Gaussian processes. To illustrate the method, we trained a LGP surrogate model on the radial distribution function of liquid neon and observed a 1,760,000-fold speed-up compared to molecular dynamics simulation, beating a conventional GP by three orders-of-magnitude. We conclude that LGPs are robust and efficient surrogate models poised to expand the application of Bayesian inference in molecular simulations to a broad spectrum of experimental data.

Dynamics of Formamide-Water Mixtures Investigated by Linear and Nonlinear Infrared Spectroscopy

Formamide (FA) exhibits complete miscibility with water, offering a simplified model for exploring the solvation dynamics of peptide linkages in biophysical processes. Its liquid state demonstrates a three-dimensional hydrogen bonding network akin to water, reflecting solvent-like behavior. Analyzing the microscopic structure and dynamics of FA-water mixtures is expected to provide crucial insights into hydrogen bonding dynamics─a key aspect of various biophysical phenomena. This study is focused on the dynamics of FA-water mixtures using linear and femtosecond infrared spectroscopies. By using the intrinsic OD stretch and extrinsic probe SCN-, the local vibrational behaviors across various FA-water compositions were systematically investigated. The vibrational relaxation of OD stretch revealed a negligible impact of FA addition on the vibrational lifetime of water molecules, underscoring the mixture's water-like behavior. However, the reorientational dynamics of OD stretch slowed with increasing FA mole fraction (XFA), plateauing beyond XFA > 0.5. This suggests a correlation between OD's reorientational time and the strength of the hydrogen bond network, likely tied to the solution's changing dielectric constant. Conversely, the vibrational relaxation dynamics of SCN- was strongly correlated with XFA, highlighting a competition between water and FA molecules in solvating SCN-. Moreover, a linear relationship between rising viscosity and the prolonged correlation time of SCN-'s slow dynamics indicates that the solution's macroscopic viscosity is dictated by the extended structures formed between FA and water molecules. The relation between the reorientation dynamics of the SCN- and the macroscopic viscosity in aqueous FA-water mixture solutions was analyzed by using the Stokes-Einstein-Debye equations. The direct viscosity-diffusion coupling is observed, which can be attributed to the homogeneous dynamics feature in FA-water mixture solutions. The inclusion of these intrinsic and extrinsic probes not only enhances the comprehensiveness of our analysis but also provides valuable insights into various aspects of the dynamics within the FA-water system. This investigation sheds light on the fundamental dynamics of FA-water mixtures, emphasizing their molecular-level homogeneity in this binary mixture 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.

Regulating polystyrene glass transition temperature by varying the hydration levels of aromatic ring/Li+ interaction

Polymer properties can be altered via lithium ion doping, whereby adsorbed Li+ binds with H2O within the polymer chain. However, direct spectroscopic evidence of the tightness of Li+/H2O binding in the solid state is limited, and the impact of Li+ on polymer sidechain packing is rarely reported. Here, we investigate a polystyrene/H2O/LiCl system using solid-state NMR, from which we determined a dipolar coupling of 11.4 kHz between adsorbed Li+ and H2O protons. This coupling corroborates a model whereby Li+ interacts with the oxygen atom in H2O via charge affinity, which we believe is the main driving force of Li+ binding. We demonstrated the impact of hydrated Li+ on sidechain packing and dynamics in polystyrene using proton-detected solid-state NMR. Experimental data and density functional theory (DFT) simulations revealed that the addition of Li+ and the increase in the hydration levels of Li+, coupled with aromatic ring binding, change the energy barrier of sidechain packing and dynamics and, consequently, changes the glass transition temperature of polystyrene.

A Nanoconfined Water-Ion Coordination Network for Flexible Energy-Dissipation Devices

Water-ion interaction in a nanoconfined environment that deeply constrains spatial freedoms of local atomistic motion with unconventional coupling mechanisms beyond that in a free, bulk state is essential to spark designs of a broad spectrum of nanofluidic devices with unique properties and functionalities. Here, it is reported that the interaction between ions and water molecules in a hydrophobic nanopore forms a coordination network with an interaction density that is nearly fourfold that of the bulk counterpart. Such strong interaction facilitates the connectivity of the water-ion network and is uncovered by corroborating the formation of ion clusters and the reduction of particle dynamics. A liquid-nanopore energy-dissipation system is designed and demonstrated in both molecular simulations and experiments that the formed coordination network controls the outflow of confined electrolytes along with a pressure reduction, capable of providing flexible protection for personnel and devices and instrumentations against external mechanical impact and attack.

Effects of Interfacial Molecular Structures on Pressure-Driven Brine Flow in Silica Mesopores

In this work, we conduct molecular dynamics simulations to investigate pressure-driven brine flow in silica mesopores under typical reservoir conditions (323 K and 20 MPa). While surface counterions accumulate strongly in the vicinity of fully deprotonated silica surfaces, water forms multilayer structures due to hydrogen bonding, counterion hydration, and excluded-volume effect. Brine flow behaviors exhibit adsorption, transition, and bulk-like regions in fully deprotonated silica mesopores, while the transition region is negligible in fully protonated ones. In the adsorption region, strong surface hydrogen bonding and a high degree of counterion hydration collectively hinder water mobility. Even without surface hydrogen bonding, persistent ion hydration impedes water flow, leading to the transition region in fully deprotonated silica mesopores and higher viscosity of brine (with 10 wt % NaCl) in the bulk region. This work elucidates the collective effects of surface chemistry and interfacial water structures on brine flow behaviors in silica mesopores from molecular perspectives.

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.

Dynamics of Acrylamide Hydrogels, Polymers, and Monomers in Water Measured with Optical Heterodyne-Detected Optical Kerr Effect Spectroscopy

The ultrafast dynamics of acrylamide monomers (AAm), polyacrylamide (PAAm), and polyacrylamide hydrogels (PAAm-HG) in water were studied using optical heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy. Previous ultrafast infrared (IR) measurements of the water dynamics showed that at the same concentration of the acrylamide moiety, AAm, PAAm, and PAAm-HG exhibited identical water dynamics and that these dynamics slowed with increasing concentration. In contrast to the IR measurements, OHD-OKE experiments measure the dynamics of both the water and the acrylamide species, which occur on different time scales. In this study, the dynamics of all the acrylamide systems slowed with increasing concentration. We found that AAm exhibits tetraexponential decays, the longest component of which followed Debye-Stokes-Einstein behavior except for the highest concentration, 40% (w/v). Low concentrations of PAAm followed a single power law decay, while high concentrations of PAAm and all concentrations of PAAm-HG decayed with two power laws. The highest concentrations, 25% and 40%, of PAAm and PAAm-HG showed nearly identical dynamics. We interpreted this result as reflecting a similar extent of chain-chain interactions. At low concentrations, PAAm displays non-Markovian, single-chain dynamics (single power law), but PAAm displays entangled chain-chain interactions at high concentrations (two power laws). PAAm-HG has chain-chain interactions at all concentrations that arise from the cross-linking. At high concentrations, the dynamics of the entangled of PAAm become identical within error as those of the cross-linked PAAm-HG.

Distinct Hydrogen Bonding Dynamics Underlies the Microheterogeneity in DMF-Water Mixtures

The hydrogen bonding interaction between the amide functional group and water is fundamental to understanding the liquid-liquid heterogeneity in biological systems. Herein, the structure and dynamics of the N,N-dimethylformamide (DMF)-water mixtures have been investigated by linear and nonlinear IR spectroscopies, using the hydroxyl stretch and extrinsic probe of thiocyanate as local vibrational reporters. According to vibrational relaxation dynamics measurements, the orientational dynamics of water is not directly tied to those of DMF molecules. Wobbling-in-a-cone analysis demonstrates that the water molecules have varying degrees of angular restriction depending on their composition due to the formation of specific water-DMF networks. Because of the preferential solvation by DMF molecules, the rotational dynamics of the extrinsic probe is slowed significantly, and its rotational time constants are correlated to the change of solution viscosity. The unique structural dynamics observed in the DMF-water mixtures is expected to provide important insights into the underlying mechanism of microscopic heterogeneity in binary mixtures.

Water Dynamics in Aqueous Poly-N-Isopropylacrylamide Below and Through the Lower Critical Solution Temperature

Poly-N-isopropylacrylamide (PNIPAM) is a thermo-responsive polymer that exhibits a reversible structural change from extended chains to aggregates in aqueous solution above its lower critical solution temperature (LCST). Using polarization-selective IR pump-probe spectroscopy, the water orientational dynamics in PNIPAM from below to above the LCST were examined and compared to those of its monomer solution, N-isopropylacrylamide (NIPAM), polyacrylamide, and an acrylamide monomer solution, which are not thermo-responsive. The OD stretch of dilute HOD in H₂O is used as a vibrational probe of water orientational dynamics. Below the LCST of the polymer, NIPAM and PNIPAM solutions exhibited identical water dynamics that were significantly different from those of bulk water, containing both faster and slower components due to solute-water interactions. Therefore, there is no difference in the nature of water interactions with a single NIPAM moiety and a long polymer chain. For all systems, including PNIPAM below and above the LCST, the orientational dynamics were modeled with a bulk water component and a polymer/monomer-associated component based on previous experimental and computational findings. Above the LCST, PNIPAM showed fast water orientational relaxation but much slower long-time dynamics compared to those of NIPAM. The slow component in PNIPAM, which was too slow to be accurately measured due to the limited OD vibrational lifetime, is ascribed to water confined in small voids (<2 nm in diameter) of PNIPAM globules. These results highlight important details about thermo-responsive polymers and the dynamics of their solvation water as they undergo a significant structural change.