Free energy and size distributions of micelles in solution

Aerosol-OT Surfactant Forms Stable Reverse Micelles in Apolar Solvent in the Absence of Water

Normal micelle aggregates of amphiphilic surfactant in aqueous solvents are formed by a process of entropically driven self-assembly. The self-assembly of reverse micelles from amphiphilic surfactant in a nonpolar solvent in the presence of water is considered to be an enthalpically driven process. Although the formation of normal and reverse surfactant micelles has been well characterized in theory and experiment, the nature of dry micelle formation, from amphiphilic surfactant in a nonpolar solvent in the absence of water, is poorly understood. In this study, a theory of dry reverse micelle formation is developed. Variation in free energy during micelle assembly is derived for the specific case of aerosol-OT surfactant in isooctane solvent using atomistic molecular dynamics simulation analyzed using the energy representation method. The existence and thermodynamic stability of dry reverse micelles of limited size are confirmed. The abrupt occurrence of monodisperse aggregates is a clear signature of a critical micelle concentration, commonly observed in the formation of normal surfactant micelles. The morphology of large dry micelles provides insight into the nature of the thermodynamic driving forces stabilizing the formation of the surfactant aggregates. Overall, this study provides detailed insight into the structure and stability of dry reverse micelles assembly in a nonpolar solvent.

Multiscale Molecular Dynamics Simulations of Model Hydrophobically Modified Ethylene Oxide Urethane Micelles

The flower-like micelles of various aggregation numbers of a model hydrophobically modified ethylene oxide urethane (HEUR) molecule, C16E45C16, and their corresponding starlike micelles, containing the surfactants C16E22 and C16E23, were studied by atomistic and coarse-grained molecular dynamic (MD) simulations. We used free energies from umbrella sampling to calculate the size distribution of micelle sizes and the average time for escape of a hydrophobic group from the micelle. Using the coarse-grained MARTINI force field, the most probable size of the model HEUR molecule was thereby determined to be about 80 hydrophobes per micelle and the average hydrophobe escape time to be about 0.1 s, both of which are consistent with previous experimental studies. Atomistic simulations reveal that hydrogen bond formation and the mean lifetime of hydration waters of the poly(ethylene oxide) (or PEO) groups are location-dependent in the HEUR micelle, with PEO groups immediately adjacent to the C16 groups forming the fewest hydrogen bonds with water and having hydration waters with longer lifetimes than those of the PEO groups located further away from the C16 groups.

Stability of biaxial nematic phase for systems with variable molecular shape anisotropy

We study the influence of fluctuations in molecular shape on the stability of the biaxial nematic phase by generalizing the mean-field model of Mulder and Ruijgrok [Physica A 113, 145 (1982)]. We limit ourselves to the case when the molecular shape anisotropy, represented by the alignment tensor, is a random variable of an annealed type. A prototype of such behavior can be found in lyotropic systems--a mixture of potassium laurate, 1-decanol, and D2O , where distribution of the micellar shape adjusts to actual equilibrium conditions. Further examples of materials with the biaxial nematic phase, where molecular shape is subject to fluctuations, are thermotropic materials composed of flexible trimericlike or tetrapodlike molecular units. Our calculations show that the Gaussian equilibrium distribution of the variables describing molecular shape (dispersion force) anisotropy gives rise to new classes of the phase diagrams, absent in the original model. Depending on properties of the shape fluctuations, the stability of the biaxial nematic phase can be either enhanced or depressed, relative to the uniaxial nematic phases. In the former case the splitting of the Landau point into two triple points with a direct phase transition line from isotropic to biaxial phase is observed.

A molecular dynamics study of free energy of micelle formation for sodium dodecyl sulfate in water and its size distribution

Free energy of micelle formation has been evaluated for spherical sodium dodecyl sulfate (SDS) in water by a thermodynamic integration method combined with a series of large-scale molecular dynamics calculations following the chemical species model. In particular, free energy change delta mu(n+1)0 with respect to the addition of one surfactant molecule to the spherical micelle of size n was obtained as a function of n. The free energy profile showed a minimum followed by a maximum, which corresponds to a peak in the size distribution. The calculated peak size n = 57 near its critical micelle concentration is in good agreement with the experimental averaged aggregation number n = 55-75 of the SDS micelle. The distribution showed a rather sharp peak, indicating that the size is almost a monodisperse one. The size is likely to be insensitive to the total concentration of the surfactant.

Classical density functional study of mixed amphiphile mesostructures

This paper presents a density functional study of mixed amphiphile solutions. The amphiphiles are modeled using four to seven fused hard spheres where the connectivity is such that the amphiphiles have a proper head group. We present calculations of lamellar and vesicle density distributions, paying particular attention to the morphology of the bilayer structures.

Density functional theory for semiflexible and cyclic polyatomic fluids

The effects of bond angle and chain stiffness on the structures of semiflexible polyatomic fluids are investigated by incorporating the bending potential into a density functional theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 2368 (2002)] that combines a modified fundamental measure theory for the excluded-volume effects and the first-order thermodynamics perturbation theory for the chain connectivity. The refined density functional theory faithfully reproduces the density profiles and conformational properties of a variety of triatomic fluids near a hard wall in which extensive Monte Carlo simulation data are available. In particular, the theory is able to capture the structures of rigid cyclic trimers where all segments are identical. The variation of local density profiles with respect to the chain length of confined polyatomic fluids is also explored. For quadratomic fluids confined in slit pores, the density profile of the middle segments exhibits novel double peaks that are absent in a fully flexible chain model. In addition, the density functional theory is applied to predicting the conformational properties and adsorption behavior of heterogeneous triatomic fluids of type "ABB" mimicking surfactant molecules. The competition between surface adsorption and self-association of trimers consisting of surface active and self-binding "A" segments and neutral "B" segment is explored.