Implicit Solvent Model Simulations of Surfactant Self-Assembly in Aqueous Solutions

A Surfactant Concentration Model for the Systematic Determination of the Critical Micellar Concentration and the Transition Width

The critical micellar concentration (cmc) is a fundamental property of surfactant solutions. Many proposed methods for the definition and determination of the cmc from property-concentration plots yield values, which depend on the studied property, on the specific technique used for its analysis and in many cases on the subjective choice of the chosen type of plot and concentration interval. In this focus review, we revise the application of a surfactant concentration model we proposed earlier that defines the cmc directly based on the surfactant concentration. Known equations for the concentration-dependence of different surfactant properties can then be combined with this concentration model and fitted to experimental data. This modular concept makes it possible to determine the cmc and the transition width in a systematic and unambiguous way. We revise its use in the literature in different contexts: the determination of the cmc of surfactants and their mixtures from different properties (electrical conductivity, NMR chemical shift, self-diffusion, surface tension, UV-Vis absorption, fluorescence intensity and fluorescence correlation). We also revise the dependence of the width of the transition region on composition, detailed studies of the properties of fluorescent probes and the aggregation of non-surfactant systems, namely amyloid peptides.

Molecular Modeling of Surfactant Micellization Using Solvent-Accessible Surface Area

We report a new implicit solvent simulation model for studying the self-assembly of surfactants, where the hydrophobic interactions were captured by calculating the relative changes of the solvent-accessible surface area (SASA) of the hydrophobic domains. Using histogram-reweighting grand canonical Monte Carlo simulations, we demonstrate that this approach allows us to match both the experimental critical micelle concentrations (cmc) and micellar aggregation numbers simultaneously with a single phenomenological surface tension γ_SASA for the poly(oxyethylene) monoalkyl ether (C _mE _n) surfactants in aqueous solutions. Excellent transferability is observed: the same model can accurately predict the experimental cmc and aggregation numbers for the C _mE _n surfactants with the alkyl lengths m between 6 and 12 and the poly(oxyethylene) lengths n between 1 and 9. The SASA-based implicit solvent model put forward in this work is general and may be applied to study more complex amphiphilic systems such as surfactants with branched alkyl chains or surfactant-hydrocarbon mixtures.

Explicit- and implicit-solvent simulations of micellization in surfactant solutions

In this article, we focus on simulation methodologies to obtain the critical micelle concentration (cmc) and equilibrium distribution of aggregate sizes in dilute surfactant solutions. Even though it is now relatively easy to obtain micellar aggregates in simulations starting from a fully dispersed state, several major challenges remain. In particular, the characteristic times of micelle reorganization and transfer of monomers from micelles to free solution for most systems of practical interest exceed currently accessible molecular dynamics time scales for atomistic surfactant models in explicit solvent. In addition, it is impractical to simulate highly dilute systems near the cmc. We have demonstrated a strong dependence of the free surfactant concentration (frequently, but incorrectly, taken to represent the cmc in simulations) on the overall concentration for ionic surfactants. We have presented a theoretical framework for making the necessary extrapolations to the cmc. We find that currently available atomistic force fields systematically underpredict experimental cmc's, pointing to the need for the development of improved models. For strongly micellizing systems that exhibit strong hysteresis, implicit-solvent grand canonical Monte Carlo simulations represent an appealing alternative to atomistic or coarse-grained, explicit-solvent simulations. We summarize an approach that can be used to obtain quantitative, transferrable effective interactions and illustrate how this grand canonical approach can be used to interpret experimental scattering results.

13C NMR chemical shift calculations of charged surfactants in water — A combined density functional theory (DFT) and molecular dynamics (MD) methodological study

Structural and magnetic properties of one anionic and one cationic amphiphile molecule (sodium octanoate and hexadecyltrimethylammonium chloride, respectively) in water are studied comparing different methods to account for the presence of the solvent. Calculated 13C NMR chemical shifts are used as the probe for accuracy of the theoretical electronic structures obtained with different descriptions of the surfactants in water solution. The best agreement with the experimental data are obtained by averaging 13C NMR isotropic chemical shifts over a large number of conformational structures of sodium octanoate while considering the electronic structure of the solvent molecules. The 13C chemical shift values of the hexadecyltrimethylammonium alkane chain are systematically overestimated by 10–15 ppm even if an extensive conformational sampling and water as the polarized continuum medium have been taken into consideration. The role of the basis set quality has been studied and discussed as well.

A model for monomer and micellar concentrations in surfactant solutions: application to conductivity, NMR, diffusion, and surface tension data

An empirical model for the concentrations of monomeric and micellized surfactants in solution is presented as a consistent approach for the quantitative analysis of data obtained with different experimental techniques from surfactant solutions. The concentration model provides an objective definition of the critical micelle concentration (cmc) and yields precise and well defined values of derived physical parameters. The use of a general concentration model eliminates subjective graphical procedures, reduces methodological differences, and thus allows one to compare directly the results of different techniques or to perform global fits. The application and validity of the model are demonstrated with electrical conductivity, surface tension, NMR chemical shift, and self-diffusion coefficient data for the surfactants SDS, CTAB, DTAB, and LAS. In all cases, the derived models yield excellent fits of the data. It is also shown that there is no need to assume the existence of different premicellar species in order to explain the chemical shifts and self-diffusion coefficients of SDS as claimed recently by some authors.

Coarse-grained molecular dynamics simulations of the sphere to rod transition in surfactant micelles

Surfactant molecules self-assemble in aqueous solutions to form various micellar structures such as spheres, rods, or lamellae. Although phase transitions in surfactant solutions have been studied experimentally, their molecular mechanisms are still not well understood. In this work, we show that molecular dynamics (MD) simulations using the coarse-grained (CG) MARTINI force field and explicit CG solvent, validated against atomistic MD studies, can accurately represent micellar assemblies of cetyltrimethylammonium chloride (CTAC). The effect of salt on micellar structures is studied for aromatic anionic salts, e.g., sodium salicylate, and simple inorganic salts, e.g., sodium chloride. Above a threshold concentration, sodium salicylate induces a sphere to rod transition in the micelle. CG MD simulations are shown to capture the dynamics of this shape transition and support a mechanism based on the reduction in the micelle-water interfacial tension induced by the adsorption of the amphiphilic salicylate ions. At the threshold salt concentration, the interface is nearly saturated with adsorbed salicylate ions. Predictions of the effect of salt on the micelle structure in different CG solvent models, namely, single-site standard water and three-site polarizable water, show qualitative agreement. This suggests that phase transitions in aqueous micelle solutions could be investigated by using standard CG water models which allow for 3 orders of magnitude reduction in the computational time as compared to that required for atomistic MD simulations.