Role of the Counterions in the Surface Tension of Aqueous Surfactant Solutions. A Computer Simulation Study of Alkali Dodecyl Sulfate Systems

Contribution of Different Molecules and Moieties to the Surface Tension in Aqueous Surfactant Solutions. II: Role of the Size and Charge Sign of the Counterions

Understanding the role of the counterion species in surfactant solutions is a complicated task, made harder by the fact that, experimentally, it is not possible to vary independently bulk and surface quantities. Here, we perform molecular dynamics simulations at constant surface coverage of the liquid/vapor interface of lithium, sodium, potassium, rubidium, and cesium dodecyl sulfate aqueous solutions. We investigate the effect of counterion type and charge sign on the surface tension of the solution, analyzing the contribution of different species and moieties to the lateral pressure profile. The observed trends are qualitatively compatible with the Hofmeister series, with the notable exception of sodium. We point out a possible shortcoming of what is at the moment, in our experience, the most realistic nonpolarizable force field (CHARMM36) that includes the parametrization for the whole series of alkali counterions. In the artificial system where the counterion and surfactant charges are inverted in sign, the counterions become considerably harder. This charge inversion changes considerably the surface tension contributions of the counterions, surfactant headgroups, and water molecules, stressing the key role of the hardness of the counterions in this respect. However, the hydration free energy gain of the counterions, occurring upon charge inversion, is compensated by the concomitant free energy loss of the headgroups and water molecules, leading to a negligible change in the surface tension of the entire system.

Surface Affinity of Alkali and Halide Ions in Their Aqueous Solution: Insight from Intrinsic Density Analysis

The surface tension of all aqueous alkali halide solutions is higher than that of pure water. According to the Gibbs adsorption equation, this indicates a net depletion of these ions in the interfacial region. However, simulations and experiments show that large, soft ions, such as I^-, can accumulate at the liquid/vapor interface. The presence of a loose hydration shell is usually considered to be the reason for this behavior. In this work, we perform computer simulations to characterize the liquid-vapor interface of aqueous alkali chloride and sodium halide solutions systematically, considering all ions from Li⁺ to Cs⁺ and from F^- to I^-. Using computational methods for the removal of surface fluctuations, we analyze the structure of the interface at a dramatically enhanced resolution, showing that the positive excess originates in the very first molecular layer and that the next 3-4 layers account for the net negative excess. With the help of a fictitious system with charge-inverted ion pairs, we also show that it is not possible to rationalize the surface affinity of ions in solutions in terms of the properties of anions and cations separately. Moreover, the surface excess is generally dominated by the smaller of the two ions.