Many-body dissipative particle dynamics simulations of micellization of sodium alkyl sulfates.

We present a study of micelle formation in alkyl sulfate surfactants using the simulation method of many-body dissipative particle dynamics (MDPD). We parametrise our model by tuning the intermolecular interactions in order to reproduce experimental values for the chemical potential and density at room temperature. Using this approach, we find that our model shows good agreement with experimental values for the critical micelle concentration (CMC). Furthermore, we show that our model can accurately predict CMC trends, which result from varying properties such as surfactant tail length and the salt concentration. We apply our model to investigate the effect of aggregation number on various micellar properties, such as the shape of individual micelles and the fraction of bound counterions. We show that micelles become aspherical at large aggregation numbers, in line with experimental predictions, and that longer tail surfactants are generally more spherical at all aggregation numbers compared to those which are shorter. We find excellent agreement between our simulations and experimental values for the degree of counterion binding, a factor that is crucial to accurately studying micellar shape, but one that is typically overlooked in the existing literature.

Using Raman Spectroscopy and Molecular Dynamics to Study Conformation Changes of Sodium Lauryl Ether Sulfate Molecules.

A study using both Raman spectroscopy and molecular dynamics (MD) simulations was carried out for alkyl ethoxysulfate (AES) surfactants at various concentrations in solution. Direct comparison between experiment and simulation shows that the conformational changes observed in MD are in good agreement with those obtained via Raman spectroscopy. We show that there is an increase in the relative number of trans conformations with increasing concentration and illustrate the relationship between phase structure and molecular conformation, which is often speculated but difficult to confirm. Our results open up the possibility of applying MD to other surfactants, with the aim of analyzing conformational behavior, which can typically be difficult to study experimentally using spectroscopy methods, due to large numbers of vibrational modes present in large complex molecules.

A many-body dissipative particle dynamics parametrisation scheme to study behaviour at air-water interfaces.

In this article, we present a general parametrisation scheme for many-body dissipative particle dynamics (MDPD). The scheme is based on matching model components to experimental surface tensions and chemical potentials. This allows us to obtain the correct surface and mixing behaviours of complex, multicomponent systems. The methodology is tested by modelling the behaviour of nonionic polyoxyethylene alkyl ether surfactants at an air/water interface. In particular, the influence of the number of ethylene oxide units in the surfactant head group is investigated. We find good agreement with many experimentally obtained parameters, such as minimum surface area per molecule; and a decrease in the surface tension with increasing surfactant surface density. Moreover, we observe an orientational transition, from surfactants lying directly on the water surface at low surface coverage, to surfactants lying parallel or tilted with respect to the surface normal at high surface coverage. The parametrisation scheme is also extended to cover the zwitterionic surfactant lauryldimethylamine oxide (LDAO), where we provide good predictions for the surface tension at maximum surface coverage. Here, if we exceed this coverage, we are able to demonstrate the spontaneous production of micelles from the surface surfactant layer.

Studying the Structure of Sodium Lauryl Ether Sulfate Solutions Using Dissipative Particle Dynamics.

Sodium lauryl ether sulfate (SLES) is a common anionic surfactant used in a large number of personal care products. Commercial products typically contain a distribution in the number of ethoxy groups; despite this, there is limited existing work studying the effect of the ethoxy groups on the phase formation and structure. This is particularly important for the effect the structure has on the viscosity, an important consideration for commercial products. Dissipative particle dynamics is used to simulate the full phase diagram of SLES in water, including both micellar and lyotropic liquid crystal phases. Phase transitions occur at locations which are in good agreement with experimental data, and we find that these boundaries can shift as a result of varying the number of ethoxy groups. Varying the ethoxy groups has a significant effect on the micellar shape and crystalline spacing, with a reduction leading to more nonspherical micelles and decreased periodic spacing of the hexagonal and lamellar phases. Finally, while typical commercial products contain a distribution of ethoxy groups, computational work tends to focus on simulations containing a single chain length. We show that it is valid to use monodisperse simulations to infer behavior about solutions with a polydisperse chain length, based on its mean molecular length.