Electrostatic potential and counterion partition between flat and spherical interfaces

J Chem Phys. 2019 Feb 21;150(7):074704. doi: 10.1063/1.5078686.

Abstract

While the electrostatic potential and the counterion distribution produced by interfaces with idealized geometries can be well-described by analytical models, the same does not hold true for the interaction between surfaces with different and arbitrary geometries. Besides, the geometry of a charged interface may also affect the counterion adsorption, potentially modulating the electrostatic potential and the solvent organization close to the interfaces, demanding molecular details to be taken into account. The complex electrostatics of a sodium dodecyl sulfate micelle in the presence of monolayers of the same surfactant at the water-vapor interface was assessed by a set of molecular dynamics simulations. The electrostatic potential was evaluated numerically, and its total magnitude was decomposed into contributions arising from each species comprising the system. The counterion adsorption was stronger at the flat interfaces due to the more favorable formation of sodium bridges, where the same counterion is bounded to two or more anionic heads, while water reorientation was more pronounced near the micelle. These opposing effects counteracted each other so that the overall electrostatic potential changes were similar for both interfaces. The increase in the counterion concentration between the micelle and the interface originates a double layer mediated repulsion amounting to a free energy barrier of at least 14 kJ/mol, preventing the micelle to get closer to the monolayers. It is noteworthy that the hydrophobic regions had electrostatic potential contributions as large as those arising from the hydrophilic regions, mostly due to the orderly orientation of the terminal methyl groups.