Theoretical and simulations-based modeling of micellization in linear and branched surfactant systems

Abstract

Surfactants are chemically-heterogeneous molecules possessing hydrophilic (head) and hydrophobic (tail) moieties. This dual nature of surfactants leads to interesting phase behavior in aqueous solution as a function of surfactant concentration, including: (i) formation of surfactant monolayers at surfaces and interfaces, and (ii) self-assembly into finite aggregates (micelles) in the bulk solution beyond the critical micelle concentration (cmc). This concentration-dependent phase behavior induces changes in solution properties. For example, the surface activity of surfactants can decrease the surface tension, and self-assembly in bulk solution can lead to changes in viscosity, equivalent conductivity, solubilization capacity, and other bulk properties. These effects make surfactants quite attractive and unique for use in product formulations, where they are utilized as detergents, dispersants, emulsifiers, solubilizers, surface and interfacial tension modifiers, and in other contexts. The specific chemical structure of the surfactant head and tail is essential in determining the overall performance properties of a surfactant in aqueous media. The surfactant tail drives the self-assembly process through the hydrophobic effect, while the surfactant head imparts a certain extent of solubility to the surfactant in aqueous solution through preferential interactions with the hydrogen-bonding network of water. The interplay between these two effects gives rise to the particular phase diagram of a surfactant, including the specific cmc at which micelles begin to form. In addition to serving as a quantitative indicator of micelle formation, the cmc represents a limit to surface monolayer formation, and hence to surface and interfacial tension reduction, because surfactant adsorption at interfaces remains approximately constant beyond the cmc. In addition, the cmc represents the onset of changes in bulk solution properties. This Thesis is concerned with the prediction of cmc's and other micellization properties for a variety of linear and branched surfactant chemical architectures which are commonly encountered in practice. Single-component surfactant solutions are investigated, in order to clarify the specific contributions of the surfactant head and tail to the free energy of micellization, a quantity which determines the cmc and all other aspects of micellization. First, a molecular-thermodynamic (MT) theory is presented which makes use of bulk-phase thermodynamics and a phenomenological thought process to describe the energetics related to the formation of a micelle from its constituent surfactant monomers. Second, a combined computer-simulation/molecular-thermodynamic (CSMT) framework is discussed which provides a more detailed quantification of the hydrophobic effect using molecular dynamics simulations. A novel computational strategy to identify surfactant head and tail using an iterative dividing surface approach, along with simulated micelle results, is proposed. Force-field development for novel surfactant structures is also discussed. Third, a statistical-thermodynamic, single-chain, mean-field theory for linear and branched tail packing is formulated, which enables quantification of the specific energetic penalties related to confinement and constraint of surfactant tails within micelles. Finally, these theoretical and simulations-based strategies are used to predict the micellization behavior of 55 linear surfactants and 28 branched surfactants. Critical micelle concentration and optimal micelle properties are reported and compared with experiment, demonstrating good agreement across a range of surfactant head and tail types. In particular, the CSMT framework is found to provide improved agreement with experimental cmc's for the branched surfactants considered.