Thermodynamically Consistent Force Field for Coarse-Grained Modeling of Aqueous Electrolyte Solution

A Coarse-Grained Model Describing the Critical Micelle Concentration of Perfluoroalkyl Surfactants in Ionic Aqueous Phase

In this study, dissipative particle dynamics (DPD) simulations were employed to determine the critical micelle concentration (CMC) of perfluoroalkyl and polyfluoroalkyl substances (PFAS) in ionic aqueous solutions. This approach provides precise CMC data for PFAS surfactants in the presence of various ionic species, thereby addressing a gap in the current literature. Additionally, this study contributes to the development of open-source molecular force fields for charged perfluorinated compounds, which are currently limited. These models incorporate hydration free energy values obtained from density functional theory (DFT) and account for ionic interactions through a well-established linear relationship. Hydrophobic interactions between the surfactant tail and water were fine-tuned to match the CMC of chosen surfactants. Then, the DPD models successfully predicted CMC values for a diverse range of surfactants, including those based on hydrocarbons and PFAS, demonstrating the ability to represent realistic salinities encountered in natural waters. Experimental validation of the methodology was conducted using sodium n-nonyl sulfate (SNS) and sodium n-dodecyl sulfate (SDS) via interfacial tension measurements, confirming the accurate representation of the CMC changes with salinity. This study enhances our understanding of the behavior of PFAS surfactants in ionic aqueous solutions and provides a valuable tool for predicting CMC values in complex environmental systems.

Partial osmotic pressures of ions in electrolyte solutions and the Gibbs-Guggenheim uncertainty principle

The concept of the partial osmotic pressure of ions in an electrolyte solution is critically examined. In principle these can be defined by introducing a solvent-permeable wall and measuring the force per unit area which can certainly be attributed to individual ions. Here I demonstrate that although the total wall force balances the bulk osmotic pressure as required by mechanical equilibrium, the individual partial osmotic pressures are extrathermodynamic quantities dependent on the electrical structure at the wall, and as such they resemble attempts to define individual ion activity coefficients. The limiting case where the wall is a barrier to only one species of ion is also considered, and with ions on both sides the classic Gibbs-Donnan membrane equilibrium is recovered, thus providing a unifying treatment. The analysis can be extended to illustrate how the electrical state of the bulk is affected by the nature of the walls and the container handling history, thus supporting the "Gibbs-Guggenheim uncertainty principle," namely the notion that the electrical state is unmeasurable and usually accidentally determined. Since this uncertainty is conferred also onto individual ion activities, it has implications for the current (2002) IUPAC definition of pH.

Impacts of targeting different hydration free energy references on the development of ion potentials

Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol^-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H⁺). The choice of ΔG°(H⁺) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K⁺, Mg²⁺, Al³⁺, and Cl^- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H⁺) for the cation and an increase in ΔG°(H⁺) for the anion. The increase in ΔG°(H⁺) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H⁺) value of -259.8 kcal mol^-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H⁺) for Mg²⁺ model is better for a reasonable description of the metal binding network.

Systematic Parameterization of Ion-Surfactant Interactions in Dissipative Particle Dynamics Using Setschenow Coefficients

Dissipative particle dynamics (DPD) simulations of nonionic surfactants with an added salt show that the Setschenow relationship is reproduced; that is, the critical micelle concentration is log-linearly dependent on the added salt concentration. The simulated Setschenow coefficients depend on the DPD bead–bead repulsion amplitudes, and matching to the experimentally determined values provides a systematic method to parameterize the interactions between salt ion beads and surfactant beads. The optimized ion-specific interaction parameters appear to be transferrable and follow the same trends as the empirical Hofmeister series.

Electroosmotic Flow Grows with Electrostatic Coupling in Confining Charged Dielectric Surfaces

In this work, the effects of polarization of confining charged planar dielectric surfaces on induced electroosmotic flow are investigated. To this end, we use dissipative particle dynamics to model solvent and ionic particles together with a modified Ewald sum method to model electrostatic interactions and surfaces polarization. A relevant difference between counterions number density profiles, velocity profiles, and volumetric flow rates obtained with and without surface polarization for moderate and high electrostatic coupling parameters is observed. For low coupling parameters, the effect is negligible. An increase of almost 500% in volumetric flow rate for moderate/high electrostatic coupling and surface separation is found when polarizable surfaces are considered. The most important result is that the increase in electrostatic coupling substantially increases the electroosmotic flow in all studied range of separations when the dielectric constant of electrolytes is much higher than the dielectric constant of confining walls. For the higher separation simulated, an increase of around 340% in volumetric flow rate when the electrostatic coupling is increased by a factor of two orders of magnitude is obtained.