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.

Anionic surfactant solutions under shear using dissipative particle dynamics.

We present a dissipative particle dynamics study of surfactant solutions under shear, which allows us to investigate their rheological properties. We consider a variety of concentrations and phase structures, including micellar solutions and liquid crystal phases. It is shown that the viscosity of micellar solutions increases as a function of concentration, in agreement with what is expected from experimental data. We also show that micelles can exhibit shear-thinning behavior when a shear force is applied, which is a result of micelles breaking down into smaller aggregates. Lamellar and hexagonal phases are found to orientate under the application of shear, in agreement with experimental observations. It is normally suggested that lamellar phases under shear can exhibit a transition between orientations as the shear rate is increased, usually as a result of lower viscosity. We calculate the viscosity for different lamellar phase orientations, showing that, although the viscosity of perpendicular orientations is lower than that of parallel orientations, we do not observe a transition to the perpendicular phase at high shear rates. Finally, we show that the choice of Schmidt number has a significant impact on the results, which is important for determining the correct behavior via simulations.

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.

Investigating anionic surfactant phase diagrams using dissipative particle dynamics: development of a transferable model.

Dissipative particle dynamics (DPD) provides a powerful coarse-grained simulation technique for the study of a wide range of soft matter systems. Here, we investigate the transferability of DPD models to the prediction of anionic surfactant phase diagrams, taking advantage of fast parameter sweeps to optimise the choice of DPD parameters for these systems. Parameters are developed which provide a good representation of the phase diagrams of SDS (sodium dodecyl sulfate) and three different isomeric forms of LAS (linear alkylbenzene sulfonates) across an extensive concentration range. A high degree of transferability is seen, with parameters readily transferable to other systems, such as AES (alkyl ether sulfates). Excellent agreement is obtained with experimentally measured quantities, such as the lamellar layer spacing. Isosurfaces are produced from the surfactant head group, from which the second moment M of the isosurface normal distribution is calculated for different phase structures. Lyotropic liquid crystalline phases are characterised by a combination of the eigenvalues of M, radial distribution functions, and visual inspections.

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.