Nonlinear dynamics in micellar surfactant solutions. I. Kinetics

Predicting surfactant phase behavior with a molecularly informed field theory

HYPOTHESIS

The computational study of surfactants and self-assembly is challenging because 1) models need to reflect chemistry-specific interactions, and 2) self-assembled structures are difficult to equilibrate with conventional molecular dynamics. We propose to overcome these challenges with a multiscale simulation approach where relative entropy minimization transfers chemically-detailed information from all-atom (AA) simulations to coarse-grained (CG) models that can be simulated using field-theoretic methods. Field-theoretic simulations are not limited by intrinsic physical time scales like diffusion and allow for rigorous equilibration via free energy minimization. This approach should enable the study of properties that are difficult to obtain by particle-based simulations.

SIMULATION WORK

We apply this workflow to sodium dodecylsulfate. To ensure chemical fidelity we present an AA force field calibrated against interfacial tension experiments. We generate CG models from AA simulation trajectories and show that particle-based and field-theoretic simulations of the CG model reproduce AA simulations and experimental measurements.

FINDINGS

The workflow captures the complex balance of interactions in a multicomponent system ultimately described by an atomistic model. The resulting CG models can study complex 3D phases like double or alternating gyroids, and reproduce salt effects on properties like aggregation number and shape transitions.

Nonlinear dynamics in micellar surfactant solutions. II. Diffusion

We discuss diffusion in micellar surfactant solutions in a form appropriate for analyzing experiments that involve large deviations from equilibrium. A general nonlinear dynamical model for inhomogeneous systems is developed that describes the effects of diffusion and micelle kinetics as a set of coupled partial differential equations for unimer concentration, micelle number concentration, average micelle aggregation number, and, optionally, the variance of the micelle aggregation number. More specialized models are developed to describe slow dynamics in situations in which the system stays in a state of partial local equilibrium or full local equilibrium. As an illustrative example of a nonlinear transport phenomenon, we discuss a simple model of diffusion from an initially homogeneous micellar solution to a rapidly created absorbing interface with fast unimer adsorption.