Recent progress in simulating microscopic ion transport mechanisms at liquid-liquid interfaces

Reparameterization of Polarizable Force Fields for Studying Ion Transfer across Liquid-Liquid Interfaces

We have developed a general scheme for refining classical polarizable molecular dynamics (MD) force fields that can accurately describe the molecular interactions in systems with liquid-liquid interfaces. While ab initio MD (AIMD) simulations can naturally describe molecular interactions, they are often so computationally expensive that simulating large system sizes and/or long time scales is usually infeasible. To resolve this, we parameterized efficient and accurate classical polarizable force fields that use AIMD reference data by minimizing both the relative entropy and the root mean squared deviation in atomic forces. We utilized our new multiscale models to study chloride ion transfer across the water-dichloromethane (DCM) interface with and without the tetraethylammonium cation as the phase-transfer catalyst. Our calculated free-energy barrier for the water-DCM interface is consistent with the other reported simulation results. We further analyzed the ion-transfer process by studying the hydration shell structures around the chloride ion and the ion-pair formation to better understand the mechanism. We observed that electronic polarizability is an important factor for the studied phase-transfer catalyst to lower the free-energy barrier of the ion transfer. Using the water-benzene interface system as an additional example, we show that our parameterization scheme provides a general route for modeling different liquid-liquid interface systems even when the experimental data or force field parameters are not readily available.

Aqueous Interfaces in Chemical Separations

Chemical separations play a vital role in refinery and reprocessing of critical materials, such as platinum group metals, rare earths, and actinides. The choice of separation system─whether it is liquid-liquid extraction (LLE), sorbents, or membranes─depends on specific needs and applications. In almost all separation processes, the desired metal ions adsorb or transfer across an aqueous interface, such as the solid/liquid interface in sorbents or oil/water interfaces in LLE. Despite these separation technologies being extensively used for decades, our understanding of the molecular-scale mechanisms governing ion adsorption and transport at interfaces remains limited. This knowledge gap presents a significant challenge in meeting the increasing demands for these critical materials due to their growing use in advanced technologies. Fortunately, recent advancements in surface-specific experimental and computational techniques offer promising avenues to bridge this gap and facilitate the development of next-generation separation systems. Interestingly, unanswered questions regarding interfacial phenomena in chemical separations hold great relevance to various fields, including energy storage, geochemistry, and atmospheric chemistry. Therefore, the model interfacial systems developed for studying chemical separations, such as amphiphilic molecules assembled at a solid/water, air/water, or oil/water interface, may have far-reaching implications, extending beyond separations and opening doors to addressing a wide range of scientific inquiries. This perspective discusses recent interfacial studies elucidating amphiphile-ion interactions in chemical separations of metal ions. These studies provide direct, molecular-scale information about solute and solvent behavior at aqueous interfaces, including multivalent and complex ions in highly concentrated solutions, which play key roles in LLE of critical materials.

Acidic Conditions Impact Hydrophobe Transfer across the Oil-Water Interface in Unusual Ways

Molecular dynamics simulation and enhanced free energy sampling are used to study hydrophobic solute transfer across the water-oil interface with explicit consideration of the effect of different electrolytes: hydronium cation (hydrated excess proton) and sodium cation, both with chloride counterions (i.e., dissociated acid and salt, HCl and NaCl). With the Multistate Empirical Valence Bond (MS-EVB) methodology, we find that, surprisingly, hydronium can to a certain degree stabilize the hydrophobic solute, neopentane, in the aqueous phase and including at the oil-water interface. At the same time, the sodium cation tends to "salt out" the hydrophobic solute in the expected fashion. When it comes to the solvation structure of the hydrophobic solute in the acidic conditions, hydronium shows an affinity to the hydrophobic solute, as suggested by the radial distribution functions (RDFs). Upon consideration of this interfacial effect, we find that the solvation structure of the hydrophobic solute varies at different distances from the oil-liquid interface due to a competition between the bulk oil phase and the hydrophobic solute phase. Together with an observed orientational preference of the hydroniums and the lifetime of water molecules in the first solvation shell of neopentane, we conclude that hydronium stabilizes to a certain degree the dispersal of neopentane in the aqueous phase and eliminates any salting out effect in the acid solution; i.e., the hydronium acts like a surfactant. The present molecular dynamics study provides new insight into the hydrophobic solute transfer across the water-oil interface process, including for acid and salt solutions.

On the mechanisms of ion adsorption to aqueous interfaces: air-water vs. oil-water

The adsorption of ions to water-hydrophobe interfaces influences a wide range of phenomena, including chemical reaction rates, ion transport across biological membranes, and electrochemical and many catalytic processes; hence, developing a detailed understanding of the behavior of ions at water-hydrophobe interfaces is of central interest. Here, we characterize the adsorption of the chaotropic thiocyanate anion (SCN^-) to two prototypical liquid hydrophobic surfaces, water-toluene and water-decane, by surface-sensitive nonlinear spectroscopy and compare the results against our previous studies of SCN^- adsorption to the air-water interface. For these systems, we observe no spectral shift in the charge transfer to solvent spectrum of SCN^-, and the Gibb's free energies of adsorption for these three different interfaces all agree within error. We employed molecular dynamics simulations to develop a molecular-level understanding of the adsorption mechanism and found that the adsorption for SCN^- to both water-toluene and water-decane interfaces is driven by an increase in entropy, with very little enthalpic contribution. This is a qualitatively different mechanism than reported for SCN^- adsorption to the air-water and graphene-water interfaces, wherein a favorable enthalpy change was the main driving force, against an unfavorable entropy change.

Structure, Thermodynamics, and Dynamics of Thiocyanate Ion Adsorption and Transfer across the Water/Toluene Interface

Molecular dynamics simulations are used to examine in detail the structure, thermodynamics, and dynamics involve in the adsorption and transfer of the thiocyanate ion (SCN^-) across the water/toluene interface. Free energy, hydration structure, and several dynamical properties as a function of the ion location along the interface normal are calculated and contrasted with recent experiments. The free energy profile exhibits a local minimum near the interface corresponding to adsorption free energy relative to bulk water of -6 kJ/mol, in reasonable agreement with experiments. The simulations provide insight into the water surface fluctuations that are coupled to the ion transfer, demonstrating formation of water finger-like structures assisting the transfer process.

A hydrogen bond-modulated soft nanoscale water channel for ion transport through liquid-liquid interfaces

Ion transport through interfaces is of ubiquitous importance in many fields such as electrochemistry, emulsion stabilization, phase transfer catalysis, liquid-liquid extraction and enhanced oil recovery. However, the knowledge of interfacial structures that significantly affect ion transport through liquid-liquid interfaces is still lacking due to the difficulty of observing nanoscale interfaces. We studied here the evolution of interfacial structures during ion transport through the decane-water interface under different ionic concentrations and external forces using molecular dynamics simulations. The roles of hydrogen bonds in ion transport through interfaces are revealed. We identified a soft nanoscale channel during ion transport through liquid-liquid interfaces and the decane phase under specific external force. The stability of the water channel and the ion transport velocity both increase with ionic concentration due to the layered ordering structures of the water near the channel surface. We observed that the stability and connectivity of the water channel in the decane phase are remarkably improved both by the high increase of the number of hydrogen bonds in the water channel with increasing ionic concentration, and by the conformational change in water molecules near the water channel surface. Our discovery of a soft nanoscale water channel by molecular simulations implies that there is a potential stable passage for ion transport through liquid-liquid interfaces.