Hydration structure in concentrated aqueous lithium chloride solutions: A reverse Monte Carlo based combination of molecular dynamics simulations and diffraction data

On the Elusive Links between Solution Microstructure, Dynamics, and Solvation Thermodynamics: Demystifying the Path through a Bridge over Troubled Conjectures and Misinterpretations

We build a fundamentally based bridge between the solute-induced microstructural perturbation of the species environment and the dynamic as well as thermodynamic responses of the fluid system, regardless of the state conditions, composition, nature of the solvent, and either the magnitude or the type of solute-solvent intermolecular-interaction asymmetries. For that purpose, we advance a fluctuation-based solvation formalism of fluid mixtures to provide meaningful descriptors of solvation phenomena, the microstructural signatures of their solute-solvent intermolecular interaction asymmetry, and the thermodynamic manifestations linked to the solution nonideality. The rigorous foundations afford us to address some crucial issues frequently invoked in the literature including the microstructural perturbation domain, its proper identification and molecular-based meaning toward the interpretation of the solvation process, and the potential impact of the local differential behavior between anions and cations on the actual salt-induced perturbation of the solvent microstructure. Indeed, we link the precisely characterized species solvation behavior to fundamental thermodynamic residual-property relations, and the dynamics associated with either the viscous flow or diffusive behavior of the solvent, to finally illustrate their outcome with experimental data of aqueous electrolyte solutions from the available literature. Ultimately, this effort provides a highly desirable unambiguous identification of the cause-effect connections between the microstructurally perturbed domains and the experimentally measured macroscopic solvation properties, including their effect on the dynamics of the solvent environment. More importantly, it lends a well-established solvation framework to bridge rigorously the microstructural details of the mixture, its dynamics, and its solvation thermodynamics to enhance our understanding of well-defined ranked Hofmeister series, i.e., by avoiding ad hoc conjectures and unsupported microscopic interpretations of solvation phenomena.

Thermodynamic calculations using reverse Monte Carlo: Simultaneously tuning multiple short-range order parameters for 2D lattice adsorption problem

Lattice simulations are an important class of problems in crystalline solids, surface science, alloys, adsorption, absorption, separation, catalysis, to name a few. We describe a fast computational method for performing lattice thermodynamic calculations that is based on the use of the reverse Monte Carlo (RMC) technique and multiple short-range order (SRO) parameters. The approach is comparable in accuracy to the Metropolis Monte Carlo (MC) method. The equilibrium configuration is determined in 5-10 Newton-Raphson iterations by solving a system of coupled nonlinear algebraic SRO growth rate equations. This makes the RMC-based method computationally more efficient than MC, given that MC typically requires sampling of millions of configurations. The technique is applied to the interacting 2D adsorption problem. Unlike grand canonical MC, RMC is found to be adept at tackling geometric frustration, as it is able to quickly and correctly provide the ordered c(2 × 2) adlayer configuration for Cl adsorbed on a Cu (100) surface.

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.

A probabilistic microkinetic modeling framework for catalytic surface reactions

We present a probabilistic microkinetic modeling (MKM) framework that incorporates the short-ranged order (SRO) evolution for adsorbed species (adspecies) on a catalyst surface. The resulting model consists of a system of ordinary differential equations. Adsorbate-adsorbate interactions, surface diffusion, adsorption, desorption, and catalytic reaction processes are included. Assuming that the adspecies ordering/arrangement is accurately described by the SRO parameters, we employ the reverse Monte Carlo (RMC) method to extract the relevant local environment probability distributions and pass them to the MKM. The reaction kinetics is faithfully captured as accurately as the kinetic Monte Carlo (KMC) method but with a computational time requirement of few seconds on a standard desktop computer. KMC, on the other hand, can require several days for the examples discussed. The framework presented here is expected to provide the basis for wider application of the RMC-MKM approach to problems in computational catalysis, electrocatalysis, and material science.

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.

Structures of ions accommodated in salty ice Ih crystals

Frozen aqueous electrolytes are ubiquitous and involved in various phenomena occurring in the natural environment. Although salts are expelled from ice during freezing of aqueous solutions, minor amounts of the constituent ions are accommodated in the crystal lattice of ice. This phenomenon was associated with the generation of the Workman-Reynolds freezing potential. Molecular simulations also confirmed the ion incorporation in the crystal lattice of ice Ih upon freezing of aqueous electrolytes and identified possible local structures of the ions. However, no experimental information is available on the structure of ions accommodated in the crystal lattice of ice Ih. In this work, we use X-ray absorption fine structure (XAFS) to study the local structures of K⁺ and Cl^- accommodated in ice Ih single crystals. Previous molecular simulations predicted that ions are trapped in the hexagonal cavities of the ice structure or replace two water molecules in the crystal lattice. Four possible configurations are considered and optimized by the calculations using ONIOM (QM/QM/QM). The results are evaluated in terms of the agreement between the experimental XAFS spectra and those simulated from the optimized structures. The spectra are most reasonably interpreted by assuming that K⁺ replaces one water molecule in the ice crystal lattice and is accommodated in a tetrahedral coordination cage. Similarly, Cl^- probably adopts the same configuration, because it explains the coordination number better than other structures, such as that assuming the replacement of two water molecules belonging to the same hexagonal planes.

Chloride ions as integral parts of hydrogen bonded networks in aqueous salt solutions: the appearance of solvent separated anion pairs

Characteristic hydrogen bonded motifs, including solvent separated anion pairs, in concentrated aqueous LiCl solutions.

Sensitivity of Electron Transfer Mediated Decay to Ion Pairing

Ion pairing in electrolyte solutions remains a topic of discussion despite a long history of research. Very recently, nearest-neighbor mediated electronic de-excitation processes of core hole vacancies (electron transfer mediated decay, ETMD) were proposed to carry a spectral fingerprint of local solvation structure and in particular of contact ion pairs. Here, for the first time, we apply electron-electron coincidence detection to a liquid microjet, and record ETMD spectra of Li 1s vacancies in aqueous solutions of lithium chloride (LiCl) in direct comparison to lithium acetate (LiOAc). A change in the ETMD spectrum dependent on the electrolyte anion identity is observed for 4.5 M salt concentration. We discuss these findings within the framework of the formation and presence of contact ion pairs and the unique sensitivity of ETMD spectroscopy to ion pairing.

Acidic properties of aqueous phosphoric acid solutions: a microscopic view

We report on new neutron and x-ray diffraction data on D2O:D3PO4 solutions at two concentrations, 1:1 and 3:1. The experimental datasets were modelled simultaneously by the reverse Monte Carlo (RMC) method. From the resulting models, partial radial distribution functions (prdf) and coordination numbers were obtained. The acidity was found to decrease with increasing D3PO4 concentration. The ratio of dissociated acidic protons was estimated by dedicated simulation runs using average coordination number constraints. It was found that in the saturated solution the ratio of dissociated protons cannot exceed 20%.