Interface Adsorption versus Bulk Micellization of Surfactants: Insights from Molecular Simulations

Surfactants play essential roles in many commonplace applications and industrial processes. Although significant progress has been made over the past decades with regard to model-based predictions of the behavior of surfactants, important challenges have remained. Notably, the characteristic time scales of surfactant exchange among micelles, interfaces, and the bulk solution typically exceed the time scales currently accessible with atomistic molecular dynamics (MD) simulations. Here, we circumvent this problem by introducing a framework that combines the general thermodynamic principles of self-assembly and interfacial adsorption with atomistic MD simulations. This approach provides a full thermodynamic description based on equal chemical potentials and connects the surfactant bulk concentration, the experimental control parameter, with the surfactant surface density, the suitable control parameter in MD simulations. Self-consistency is demonstrated for the nonionic surfactant C₁₂EO₆ (hexaethylene glycol monododecyl ether) at an alkane/water interface, for which the adsorption and pressure isotherms are computed. The agreement between the simulation results and experiments is semiquantitative. A detailed analysis reveals that the used atomistic model captures well the interactions between surfactants at the interface but less so their adsorption affinities to the interface and incorporation into micelles. Based on a comparison with other recent studies that pursued similar modeling challenges, we conclude that the current atomistic models systematically overestimate the surfactant affinities to aggregates, which calls for improved models in the future.