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

DPD simulations of anionic surfactant micelles: a critical role for polarisable water models

We investigate the effects of polarisable water models in dissipative particle dynamics (DPD) simulations, focussing on the influence these models have on the aggregation behaviour of sodium dodecyl sulfate solutions. Studies in the literature commonly report that DPD approaches underpredict the micellar aggregation number of ionic surfactants compared to experimental values. One of the proposed reasons for this discrepancy is that existing water models are insufficient to accurately model micellar solutions, as they fail to account for structural changes in water close to micellar surfaces. We show that polarisable DPD water models lead to more realistic counterion behaviour in micellar solutions, including the degree of counterion disassociation. These water models can also accurately reproduce changes in the dielectric constant of surfactant solutions as a function of concentration. We find evidence that polarisable water leads to the formation of more stable micelles at higher aggregation numbers. However, we also show that the choice of water model is not responsible for the underestimated aggregation numbers observed in DPD simulations. This finding addresses a key question in the literature surrounding the importance of water models in DPD simulations of ionic micellar solutions.

Microstructure-dependent CO2-responsive microemulsions for deep-cleaning of oil-contaminated soils

CO2-responsive microemulsion (ME) is considered a promising candidate for deep-cleaning and oil recovery from oil-contaminated soils. Understanding the responsive nature of different microstructures (i.e., oil-in-water (O/W), bicontinuous (B.C.) and water-in-oil (W/O)) is essential for unlocking the potential and mechanisms of CO2-responsive emulsions in complex multiphase systems and providing comprehensive guidance for remediation of oil-contaminated soils. Herein, the responsiveness of microstructures of ME to CO2 trigger was investigated using experimental designs and coarse-grained molecular dynamic simulations. MEs were formed for the first time by a weakly associated pseudo-Gemini surfactant of indigenous organic acids (naphthenic acids, NAs are a class of natural surface-active molecules in crude oil) and tetraethylenepentamine (TEPA) through fine tuning of co-solvent of dodecyl benzene sulfonic acid (DBSA) and butanol. The O/W ME exhibited an optimal CO2-responsive character due to easier proton migration in the continuous aqueous phase and more pronounced dependence of configuration on deprotonated NA ions. Conversely, the ME with W/O microstructure exhibited a weak to none responsive characteristic, most likely attributed to its high viscosity and strong oil-NA interactions. The O/W ME also showed superior cleaning efficiency and oil recovery from oil-contaminated soils. The results from this study provide insights for the design of CO2-responsive MEs with desired performance and guidance for choosing the favorable operating conditions in various industrial applications, such as oily solid waste treatment, enhanced oil recovery (EOR), and pipeline transportation. The insights from this work allow more efficient and tailored design of switchable MEs for manufacturing advanced responsive materials in various industrial sectors and formulation of household products.

Temperature-Derived Purification of Gold Nano-Bipyramids for Colorimetric Detection of Tannic Acid

Gold nanostructures have attracted broad attention. Among various nanostructures, gold nanobipyramids have shown great potential in sensing, biomedicine, environmental protection, chemical catalysis, and optics due to their unique physical and optical properties and ease of chemical functionalization. Compared with other plasmonic nanostructures, gold nanobipyramids possess narrow optical resonances, stronger plasmonic local field enhancement, and size- and shape-dependent surface plasmon resonance. However, the synthesis and purification of homogeneous gold nanobipyramids are very challenging. The gold nanobipyramids synthesized via the commonly used seed-mediated growth method have low yields and are often coproduced with spherical nanoparticles. In this study, we reported a temperature-derived purification method for the isolation of gold bipyramids. In the presence of salt, by altering the temperature of the solution, large gold bipyramids can be separated from small spherical nanoparticles. As a result, a yield of as high as 97% gold nanobipyramids can be achieved through a single round of purification, and correspondingly, the ratio between the longitudinal surface plasmon resonance (LSPR) and transverse SPR intensity significantly increases to as high as 6.7. The purified gold nanobipyramids can be used as a colorimetric probe in the detection of tannic acid with a detection limit of 0.86 μM and a linear detection range from 1.25 to 37.5 μM.

Molecular Insights on Morphology, Composition, and Stability of Mixed Micelles Formed by Ionic Surfactant and Nonionic Block Copolymer in Water Using Coarse-Grained Molecular Dynamics Simulations

The nanoscale association domains are the ultimate determinants of the macroscopic properties of complex fluids involving amphiphilic polymers and surfactants, and hence, it is foremost important to understand the role of polymer/surfactant concentration on these domains. We have used coarse-grained molecular dynamics simulations to investigate the effect of polymer/surfactant concentration on the morphology of poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) (PEO-PPO-PEO, i.e., pluronics or poloxamers) block copolymer, and ionic surfactants sodium dodecyl sulfate (SDS), mixed micelles in aqueous solution. The proclivity of the surfactant to form the mixed micelles is also probed using umbrella sampling simulations. In this study, we observed that the core of the pluronic + SDS formed mixed micelles consists of PPO, the alkyl tail of SDS, and some water molecules, whereas the PEO, water, and sulfate headgroups of SDS form a shell, consistent with experimental observations. The micelles are spherical at high-pluronic/low-SDS compositions, ellipsoidal at high-SDS/low-pluronic compositions, and wormlike-cylindrical at high-pluronic/high-SDS compositions. The transitions in micelle morphology are governed by the solvent accessible surface area of mixed aggregates, electrostatic repulsion between SDS-headgroups, and dehydration of PEO and PPO segments. The free energy barrier for SDS escape is much higher in mixed micelles than in pure SDS micelles, indicating a stronger tendency for SDS to form pluronic-SDS mixed micelles.

Probing self-assembled micellar topologies via micro-scale diffusive dynamics of surfactants

HYPOTHESIS

Surfactants spontaneously self-assemble in aqueous solutions and are critical in energy, biotechnology, and the environment. The self-assembled micelles may experience distinct topological transitions beyond a critical counter-ion concentration, yet the associated mechanical signatures are identical. By monitoring self-diffusion dynamics of individual surfactants in micelles via a non-invasive ¹H NMR diffusometry, we may distinguish various topological transitions overcoming challenges associated with traditional microstructural probing techniques.

EXPERIMENTS

Three micellar systems based on CTAB/5mS, OTAB/NaOA and CPCl/NaClO₃ are considered at various counter-ion concentrations, and their rheological properties are assessed. A systematic ¹H NMR diffusometry is conducted and the resulting signal attenuation is measured.

FINDINGS

With no counter-ion, surfactants self-diffuse freely with a mean squared displacement Z²∼T_diff in the micelles. As counter-ion concentration increases, self-diffusion becomes restricted with Z²∼T_diff^α, and α→0.5. Beyond the viscosity peak, for the OTAB/NaOA system that shows a linear-shorter linear micelle transition, Z²∼T_diff^0.5. Conversely, for the CTAB/5mS system that experiences a linear wormlike-vesicle transition above the viscosity peak, a free self-diffusion is recovered. The diffusion dynamics in CPCl/NaClO₃ are similar to those of OTAB/NaOA. Hence, a similar topological transition is surmised. These results highlight the unique sensitivity of the ¹H NMR diffusometry to micelles topological transitions.

Tween-80 on Water/Oil Interface: Structure and Interfacial Tension by Molecular Dynamics Simulations

Surfactants are used in many fields of the chemical industry in a wide range of applications. Generally, surfactants on water/oil interfaces reduce interfacial tension, enabling the formation of emulsions or providing greater stability to the emulsion formed. Although molecular dynamics has been extensively used and has achieved remarkable success in describing thermodynamic and molecular properties of systems with surface-active compounds, the traditional molecular simulation force fields considerably constrain the system size and the time scale of simulations. Here, we propose a coarse-grained model of polysorbate 80, a nonionic surfactant commercially known as Tween-80. Based on the proposed coarse-grained model, we evaluate the influence of the more internal ethylene oxide chain on the properties of the water/Tween-80/decane system. We verify with the simulation results that the model can reproduce the expected decrease in interfacial tension as the surfactant quantity increases in the simulations. Furthermore, we observe changes in the surfactant orientation as their quantity in the interface increases, indicating a preferential orientation for these molecules in the adsorption layer. We also assessed the partition coefficient of the surfactant between the two bulk phases by performing free energy calculations, which showed a higher affinity of Tween-80 surfactants with the water phase.

Design and simulation of a caprylic acid enzymatically modified phosphatidylcholine micelle using a coarse-grained molecular dynamics simulations approach

Computationally simulated micelle models provide useful structural information on the molecular and biological sciences. One strategy to study the self-aggregation process of surfactant molecules that make up a micelle is through molecular dynamics (MD) simulations. In this study, a theoretical approach with a coarse-grained MD simulation (CG-MD) was employed to evaluate the critical micellar concentration (CMC), the micellization process, building a tridimensional (3D) model system of a micelle using data from the experimentally enzymatically modified phospholipids (PL) by phospholipase A1 (PA1). This required enzymatic interesterification of soybean phosphatidylcholine (PC) with caprylic acid, along with purification and characterization by chromatographic techniques to measure the esterified fatty acids and the corresponding PL composition. The number of molecules used in the CG-MD simulation system was determined from the experimental CMC data which was 0.025%. The molecular composition of the system is: 1 C 18:2, 2 C 8:0/8:0, 3 C 8:0/18:3n-9, 4 C 8:0/18:0, 5 C8:0/18:2n-6, 6 C8:0/18:1n-9, and 7 C 8:0/16:0. According to our theoretical results, the micelle model is structurally stable with an average Rg of 3.64 ± 0.10 Å, and might have an elliptical form with a radius of 24.6 Å. Regarding CMC value there was a relationship between the experimental data of the modified PLs and the theoretical analysis by GC-MD, which suggest that the enzymatic modification of PLs does not affect their self-aggregation properties. Finally, the micellar system obtained in the current research can be used as a simple and useful model to design optimal biocompatible nanoemulsions as possible vehicles for bioactive small molecules.Communicated by Ramaswamy H. Sarma.

Rationalizing the Design of Pluronics-Surfactant Mixed Micelles through Molecular Simulations and Experiments

Aqueous systems comprising polymers and surfactants are technologically important complex fluids with tunable features dependent on the chemical nature of each constituent, overall composition in mixed systems, and solution conditions. The phase behavior and self-assembly of amphiphilic polymers can be changed drastically in the presence of conventional ionic surfactants and need to be clearly understood. Here, the self-aggregation dynamics of a triblock copolymer (Pluronics L81, EO₃PO₄₃EO₃) in the presence of three cationic surfactants (with a 12C long alkyl chain but with different structural features), viz., dodecyltrimethylammonium bromide (DTAB), didodecyldimethylammonium bromide (DDAB), and ethanediyl-1,2-bis(dimethyldodecylammonium bromide) (12-2-12), were investigated in an aqueous solution environment. The nanoscale micellar size expressed as hydrodynamic diameter (D_h) of copolymer-surfactant mixed aggregates was evaluated using dynamic light scattering, while the presence of a varied micellar geometry of L81-cationic surfactant mixed micelles were probed using small-angle neutron scattering. The obtained findings were further validated from molecular dynamics (MD) simulations, employing a simple and transferable coarse-grained molecular model based on the MARTINI force field. L81 remained molecularly dissolved up to ∼20 °C but phase separated, forming turbid/translucent dispersion, close to its cloud point (CP) and existed as unstable vesicles. However, it exhibited interesting solution behavior expressed in terms of the blue point (BP) and the double CP in the presence of different surfactants, leading to mixed micellar systems with a triggered morphology transition from unstable vesicles to polymer-rich micelles and cationic surfactant-rich micelles. Such an amendment in the morphology of copolymer nanoaggregates in the presence of cationic surfactants has been well observed from scattering data. This is further rationalized employing the MD approach, which validated the effective interactions between Pluronics-cationic surfactant mixed micelles. Thus, our experimental results integrated with MD yield a deep insight into the nanoscale interactions controlling the micellar aggregation (Pluronics-rich micelles and surfactant-rich micelles) in the investigated mixed system.

Prediction of Micellar Thermodynamics of Nonionic Surfactants Based on the Square Gradient Theory

A geometry-dependent contribution based on the square gradient theory of van der Waals is proposed as a predictive modification of the interfacial energy contribution for the micellar thermodynamic theory. The model has an analytic prediction for the spherical and cylindrical geometries. For ellipsoidal geometry, a simple yet physically meaningful approximation is proposed. The critical micelle concentration (CMC) and the surface tension isotherm under the new contribution are compared with the classical theory. The modified model describes qualitatively the available experimental data and the surface isotherm, showing an improvement in the predictions of the CMC.

Coarse-Grained Model Incorporating Short- and Long-Range Effective Potentials for the Fast Simulation of Micelle Formation in Solutions of Ionic Surfactants

A coarse-grained model comprising short- and long-range effective potentials, parametrized with the iterative Boltzmann inversion (IBI) method, is presented for capturing micelle formation in aqueous solutions of ionic surfactants using as a model system sodium dodecyl sulfate (SDS). In the coarse-grained (CG) model, each SDS molecule is represented as a sequence of four beads while each water molecule is modeled as a single bead. The proposed CG scheme involves ten potential energy functions: four of them describe bonded interactions and control the distribution functions of intramolecular degrees of freedom (bond lengths, valence angles, and dihedrals) along an SDS molecule while the other six account for intermolecular interactions between pairs of SDS and water beads and control the radial distribution functions. The nonbonded effective potentials between coarse-grained SDS molecules extend up to about 12 nm and capture structural and morphological features of the micellar solution both at short and long distances. The long-range component of these potentials, in particular, captures correlations between surfactant molecules belonging to different micelles and is essential to describe ordering associated with micelle formation. A new strategy is introduced for determining the effective potentials through IBI by using information (target distribution functions) extracted from independent atomistic simulations of a micellar reference system (a salt-free SDS solution at total surfactant concentration c_T equal to 103 mM, temperature T equal to 300 K, and pressure P equal to 1 atm) obtained through a multiscale approach described in an earlier study. It employs several optimization steps for bonded and nonbonded interactions and a gradual parametrization of the short- and long-range components of the latter, followed by reparametrization of the bonded ones. The proposed CG model can reproduce remarkably accurately the microstructure and morphology of the reference system within only a few hours of computational time. It is therefore very promising for future studies of structural and morphological behavior of various liquid surfactant formulations.

Scission energy and topology of micelles controlled by the molecular structure of additives

We employ coarse-grained (CG) molecular dynamics simulations (MD) to investigate the effects of the molecular structure of additives on the scission energy and morphology of charged micelles. Considering sodium dodecyl sulfate (SDS) as a representative charged surfactant and taking trimethylphenylammonium chloride (TMPAC) and octyltrimethylammonium bromide (OTAB) as oppositely charged additives, we show that the scission energy and topology of micelles vary significantly depending on the molecular structure of the hydrophobic part of the additives. The cyclic aromatic tail of the TMPAC disrupts the core structure of the SDS micelle and hence decreases the micelle scission energy, whereas the linear alkyl tail of the OTAB packs very well with the micelle core and increases the scission energy. Although both the additives have similar head structures, they lead to very different micelle morphologies because of the difference in the shape of their tail structures; ring-like or toroidal shaped micelles are formed in SDS/TMPAC solution whereas bicelle-like structures are formed in SDS/OTAB solution when the additive to surfactant ratio is higher than a certain value.

The phase behaviour of cetyltrimethylammonium chloride surfactant aqueous solutions at high concentrations: an all-atom molecular dynamics simulation study

We explore the phase behaviour of aqueous solutions of the cetyltrimethyl ammonium chloride (CTAC) surfactant and in particular the transition from the micellar phase (L₁) to the hexagonal columnar phase (H₁) by employing all-atom (AA) molecular dynamics (MD) simulations for six CTAC concentrations in the range of 34.1 wt% to 70.5 wt%, at the temperature of 318 K and pressure of 1 atm. For the concentrations considered, we examine the spontaneous occurrence of the H₁ phase by testing a number of plausible values for the linear density (molecules per unit length) along the cylindrical columns. Using large simulation cells and starting from random initial configurations, the MD simulations demonstrate that the micellar phase occurs for concentrations up to 50.0 wt%, with CTAC molecules self-assembling into a mixture of spherical and rod-like micelles. At even higher concentrations, the system self-organizes into the H₁ phase in accordance with the available experimental data. For the analysis of the MD trajectories, we devise a clustering algorithm based on Voronoi tesselation which enables (a) the thorough characterization of the shape and structure of both molecules and assemblies, and (b) the investigation of the positional and orientational order in the system that are further scrutinised using radial pair correlation functions and X-ray diffraction patterns. Our work paves the way for the investigation of the phase behaviour at high concentrations of other surfactants.

Ion Atmosphere of Wormlike Micelles Profiled by Contrast Variation Small-Angle Neutron Scattering

Structural studies of wormlike micelles have so far mostly focused on the conformational properties of surfactant aggregates. The diffuse ionic atmosphere, which has a profound influence on various micellization phenomena such as thermodynamic stability and structural polymorphism, remains largely unexplored experimentally. In this report a strategy of contrast variation small-angle neutron scattering for this crucial structural study is outlined. Underlined by a general criterion established for unbiasedly identifying the length scale relevant to charge association from the spectral evolution, our analytical framework can provide a quantitative description of counterion distribution in a mathematically tractable manner. Our method can be conveniently extended to facilitate structural studies of complex multicomponent systems using contrast variation neutron scattering.

Temperature-Transferable Coarse-Grained Model for Poly(propylene oxide) to Study Thermo-Responsive Behavior of Triblock Copolymers

Thermo-responsive behavior of ethylene oxide (EO)-propylene oxide (PO) copolymers makes them suitable for many potential applications. Reproducing the origins of the tunable properties of EO-PO copolymers using coarse-grained (CG) models such as the MARTINI force field is critically important for building a better understanding of their behavior. In the present work, we have investigated the effects of coarse-graining on the water-polymer interaction across a temperature range. We compared the performance of different all-atom force fields to find the most appropriate one for the purpose of PO block parameterization in the MARTINI platform. We parameterized a CG temperature-dependent PO model based on the reproduction of the atomistic free energy of transfer of propylene oxide trimer from octane to water over a range of temperatures (20-60 °C) and compared the atomistic bond and angle distributions. Then, we used the model to study the effects of EO/PO ratio, molecular weight, and concentration on the thermo-responsive behavior of EO-PO copolymers in water. The results show an excellent agreement with experiments in different areas. Our temperature-dependent model reproduces (1) micellar phase above critical micelle temperature (CMT) and unimer phase below CMT for different Pluronics (a class of EO-PO triblock copolymers) spanning many EO/PO ratios and molecular weights; (2) spherical-to-rodlike micellar shape transition for Pluronics with 60 wt % of PO content or more; (3) diffusion coefficients for Pluronics with high PO content (P104 Pluronic with a PO mass of 3500 g mol^-1) across a broad range of temperatures; and (4) micelle core size and micelle diameter similar to experimental results. Overall, our model improves the temperature sensitivity of EO-PO copolymers of existing models significantly, particularly for copolymers that are dominated by PO agents.

High-viscosity Pickering emulsion stabilized by amphiphilic alginate/SiO2 via multiscale methodology for crude oil-spill remediation

The separation of crude oil from oily water and collection of the emulsion constituents has attracted significant attention. We demonstrate that the relationships between inherent dynamic factors and the performance of a Pickering emulsion stabilized by SiO₂ particles with adsorbed hydrophobically modified sodium alginate derivatives (HMSA), a natural pH-sensitive polysaccharide, can be clarified via a multi-scale methodology. Functionalization of the silica surface with HMSA controls particle dispersibility, as verified by turbidity and stability analyses, the zeta potential, and transmission electron microscopy measurements. The interaction mechanism between HMSA and SiO₂ nanoparticles was elucidated by both experimental adsorption measurements and computer simulations, which showed qualitative consistency. The aggregation/disaggregation of HMSA/SiO₂ particles achieved by tuning the pH of the solution facilitated reversible dispersibility/collectability behavior. Overall, a high-viscosity Pickering emulsion system based on particle-particle and droplet-droplet interactions, which can be filtered for the recovery of spilled crude oil, was demonstrated.

Self-assembly of the cationic surfactant n-hexadecyl-trimethylammonium chloride in methyltrimethoxysilane aqueous solution: classical and reactive molecular dynamics simulations

A flexible aerogel polymerized from methyltrimethoxysilane (MTMS) shows great promise as a high-performance insulator owing to its substantially low thermal conductivity and mechanical flexibility, attributed to its porous microstructure and organic-inorganic hybridization, respectively, which promote its industrial applications. Conventionally, the cationic surfactant n-hexadecyltrimethylammonium chloride (CTAC) is utilized to experimentally control the nanoscale microstructure and, consequently, the flexibility of the MTMS aerogel; however, the mechanism through which CTAC prevents MTMS aggregation in the solution is not yet fully understood. This study unravels the role of CTAC in preventing MTMS aggregation in aqueous solution using both classical and reactive molecular dynamics simulations. We found that CTAC molecules can form self-aggregates even when the polymerization of MTMS progresses and then the MTMS-derived oligomer turns to be hydrophobic in aqueous solution. In summary, the self-assemblies of CTAC disperse among the MTMS associations and effectively prevent MTMS clustering, and this is considered as the key mechanism underlying the formation of a flexible microstructure of the hybrid aerogel.

Surfactant desorption and scission free energies for cylindrical and spherical micelles from umbrella-sampling molecular dynamics simulations

HYPOTHESIS

The free energies associated with adsorption/desorption of individual surfactants from micelles and the fusion/scission of long micelles can be used to estimate the rate constants for micellar kinetics as functions of surfactant and salt concentration.

EXPERIMENTS

We compute the escape free energies △G_esc of surfactant from micelles and the scission free energies △G_sciss of long micelles from coarse-grained molecular dynamics simulations coupled with umbrella sampling, for micelles of both sodium dodecylsulfate (SDS) in sodium chloride (NaCl) and cetyltrimethylammonium chloride (CTAC) in sodium salicylate (NaSal).

FINDINGS

For spherical micelles, △G_esc values have maxima at certain aggregation numbers, and at salt-to-surfactant molar concentration ratios R near unity, consistent with experiments. For cylindrical micelles, SDS/NaCl shows a minimum, and CTAC/NaSal a maximum in △G_esc, both at R ~ 0.7, while △G_sciss of CTAC micelles also peaks at around R ~ 0.7 and that of SDS micelles increases monotonically with R. We explain the non-monotonic dependence of escape and scission free energies on R by a combination of electrostatic screening and the decrease of micelle radius with increasing R. Transitions from predominantly spherical to cylindrical micelles, and between adsorption/desorption and fusion/scission kinetics with changing salt concentration can be inferred from the free energies for CTAC/NaSal.

Tuning Cationic Micelle Properties with an Antioxidant Additive: A Molecular Perspective

In this work, we characterize the micellization and morphology transition induced in aqueous cetyltrimethylammonium bromide (CTAB) solution by the addition of the antioxidant propyl gallate (PG) using tensiometry, rheology, and small-angle neutron scattering (SANS) techniques combined with the molecular dynamics (MD) simulation approach. The adsorption of CTAB at the air-water interface in the presence of varying [PG] revealed a progressive decrease in the critical micelle concentration (CMC), while the changes in different interfacial parameters indicated enhancement of the hydrophobicity induced by PG in the CTAB micellar system. The dynamic rheology behavior indicated an increase in the flow viscosity (η) as a function of [PG]. Moreover, the rheological components (storage modulus, G', and loss modulus, G″) depicted the viscoelastic features. SANS measurements depicted the existence of ellipsoidal micelles with varying sizes and aggregation number (Nagg) as a function of [PG] and temperature. Computational simulation performed using density functional theory (DFT) calculations and molecular dynamics (MD) provided an insight into the atomic composition of the examined system. The molecular electrostatic potential (MEP) analysis depicted a close proximity of CTAB, i.e., emphasized favorable interactions between the quaternary nitrogen of CTAB and the hydroxyl group of the PG monomer, further validated by the two-dimensional nuclear Overhauser enhancement spectroscopy (2D-NOESY), which showed the penetration of PG inside the CTAB micelles. In addition, various dynamic properties, viz., the radial distribution function (RDF), the radius of gyration (Rg), and solvent-accessible surface area (SASA), showed a significant microstructural evolution of the ellipsoidal micelles in the examined CTAB-PG system, where the changes in the micellar morphology with a more elongated hydrophobic chain and the increased Rg and SASA values indicated the notable intercalation of PG in the CTAB micelles.

Molecular Basis for the Morphological Transitions of Surfactant Wormlike Micelles Triggered by Encapsulated Nonpolar Molecules

Surfactant wormlike micelles are prone to experience morphological changes, including the transition to spherical micelles, upon the addition of nonpolar additives. These morphological transitions have profound implications in diverse technological areas, such as the oil and personal-care industries. In this work, additive-induced morphological transitions in wormlike micelles were studied using a molecular theory that predicts the equilibrium morphology and internal molecular organization of the micelles as a function of their composition and the molecular properties of their components. The model successfully captures the transition from wormlike to spherical micelles upon the addition of a nonpolar molecule. Moreover, the predicted effects of the concentration, molecular structure, and degree of hydrophobicity of the nonpolar additive on the wormlike-to-sphere transition are shown to be in good agreement with experimental trends in the literature. The theory predicts that the location of the additive in the micelle (core or hydrophobic-hydrophilic interface) depends on the additive hydrophobicity and content, and the morphology of the micelles. Based on the results of our model, simple molecular mechanisms were proposed to explain the morphological transitions of wormlike micelles upon the addition of nonpolar molecules of different polarities.

Self-assembly and viscosity changes of binary surfactant solutions: A molecular dynamics study

HYPOTHESIS

Structure and self-assembly of surfactants in the solution shows a fundamental influence on its viscosity. Through molecular simulations using Martini force field, synergistic effects in aggregation as well as the viscosity changes of a binary ionic surfactant system can be modelled. Simulations: Coarse-grained molecular dynamics simulations are performed to model the SDS/CAPB binary surfactant solution, and both equilibrium and non-equilibrium methods are used to calculate the viscosity of the equilibrated micellar systems.

FINDINGS

Our simulation results indicate that the new version of the Martini force field can provide more reasonable self-assembly of surfactant, both single and binary system. Synergistic effects in micelle formation for SDS/CAPB have been successfully reproduced, that is, the formation of cylindrical micelles or even wormlike micelles at a lower concentration when compared with the pure system. Meanwhile, both equilibrium and non-equilibrium methods provide quantitatively comparable viscosity for each system. For pure micellar system, the viscosity linearly increases with the total concentration. Nevertheless, our simulation fails to capture the viscosity enhancement of the solution in corresponding with the formation of rodlike or wormlike micelles, and a full parameter optimization of force field is still necessary.

Supramolecular Packing Drives Morphological Transitions of Charged Surfactant Micelles

The shape and size of self-assembled structures upon local organization of their molecular building blocks are hard to predict in the presence of long-range interactions. Combining small-angle X-ray/neutron scattering data, theoretical modelling, and computer simulations, sodium dodecyl sulfate (SDS), over a broad range of concentrations and ionic strengths, was investigated. Computer simulations indicate that micellar shape changes are associated with different binding of the counterions. By employing a toy model based on point charges on a surface, and comparing it to experiments and simulations, it is demonstrated that the observed morphological changes are caused by symmetry breaking of the irreducible building blocks, with the formation of transient surfactant dimers mediated by the counterions that promote the stabilization of cylindrical instead of spherical micelles. The present model is of general applicability and can be extended to all systems controlled by the presence of mobile charges.

A combined experimental and computational approach reveals how aromatic peptide amphiphiles self-assemble to form ion-conducting nanohelices

Reported here is a combined experimental-computational strategy to determine structure-property-function relationships in persistent nanohelices formed by a set of aromatic peptide amphiphile (APA) tetramers with the general structure K S XEK S , where KS= S-aroylthiooxime modified lysine, X = glutamic acid or citrulline, and E = glutamic acid. In low phosphate buffer concentrations, the APAs self-assembled into flat nanoribbons, but in high phosphate buffer concentrations they formed nanohelices with regular twisting pitches ranging from 9-31 nm. Coarse-grained molecular dynamics simulations mimicking low and high salt concentrations matched experimental observations, and analysis of simulations revealed that increasing strength of hydrophobic interactions under high salt conditions compared with low salt conditions drove intramolecular collapse of the APAs, leading to nanohelix formation. Analysis of the radial distribution functions in the final self-assembled structures led to several insights. For example, comparing distances between water beads and beads representing hydrolysable KS units in the APAs indicated that the KS units in the nanohelices should undergo hydrolysis faster than those in the nanoribbons; experimental results verified this hypothesis. Simulation results also suggested that these nanohelices might display high ionic conductivity due to closer packing of carboxylate beads in the nanohelices than in the nanoribbons. Experimental results showed no conductivity increase over baseline buffer values for unassembled APAs, a slight increase (0.4 × 102 μS/cm) for self-assembled APAs under low salt conditions in their nanoribbon form, and a dramatic increase (8.6 × 102 μS/cm) under high salt conditions in their nanohelix form. Remarkably, under the same salt conditions, these self-assembled nanohelices conducted ions 5-10-fold more efficiently than several charged polymers, including alginate and DNA. These results highlight how experiments and simulations can be combined to provide insight into how molecular design affects self-assembly pathways; additionally, this work highlights how this approach can lead to discovery of unexpected properties of self-assembled nanostructures.

Molecular Dynamics Simulations and Density Functional Theory on Unraveling Photoresponsive Behavior of Wormlike Micelles Constructed by 12-2-12·2Br- and trans-ortho-Methoxy Cinnamate

Photoresponsive systems with controllable self-assembly morphologies and adjustable rheological properties have attracted widespread interest by researchers in the past few years. Among them, the photoresponsive systems consisting of ortho-methoxycinnamic (OMCA) and Gemini surfactants are endowed with rich self-assemblies with different states and in different scales including spherical micelles, wormlike micelles, vesicles, aqueous two-phase system (ATPS), etc. All these self-assemblies display excellent photoresponsive behavior. However, the mechanism of these photoresponsive behaviors has not been unraveled systematically so far. In this study, molecular dynamics (MD) simulations, density functional theory (DFT) calculations, transmission electron microscopy, and rheology are employed to investigate the photoresponsive behaviors of wormlike micelles caused by photoisomerization of trans-OMCA in 12-2-12·2Br^-/trans-OMCA solutions and to unravel the underlying mechanisms of these photoresponsive behaviors. The experimental results show that 12-2-12·2Br^-/trans-OMCA micelles display photoresponsiveness after UV-light irradiation, with the transformation of micellar morphologies from wormlike micelle to spherical micelles. In MD simulations, certain micelle morphologies in experiments and the specific packing between 12-2-12·2Br^-/OMCA were successfully captured. The larger three-dimensional structure and steric hindrance of cis-OMCA disturb the interior structure of micelles. The stronger hydrophilicity of cis-OMCA induces the escape of cis-OMCA from the interval of micelles to the solution. The energy results prove that trans-OMCA associates more strongly with 12-2-12·2Br^- than cis-OMCA. These causes lead to the fission and repacking of wormlike micelles.

Simulation of diblock copolymer surfactants. II. Micelle kinetics

Molecular dynamics (MD) simulations are used to measure dynamical properties of a simple bead-spring model of A-B diblock copolymer molecules, and to characterize rates and mechanisms of several dynamical processes. Dynamical properties are analyzed within the context of a kinetic population model that allows for both stepwise insertion and expulsion of individual free molecules and occasional fission and fusion of micelles. Kinetic coefficients for stepwise processes and micelle fission have been extracted from MD simulations of individual micelles. Insertion of a free surfactant molecule into a preexisting micelle is shown to be a completely diffusion-controlled process for the model studied here. Estimates are given for rates of rare events that create and destroy entire micelles by competing mechanisms involving stepwise association and dissociation or fission and fusion. Both mechanisms are shown to be relevant over the range of parameters studied here, with association and dissociation dominating in systems with more soluble surfactants and fission and fusion dominating in systems with less soluble surfactants.

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.

Stretch and Breakage of Wormlike Micelles under Uniaxial Strain: A Simulation Study and Comparison with Experimental Results

We use coarse-grained (CG) molecular dynamics simulations to determine the effect of uniaxial strain on the stress, scission stress, and scission energy of solutions of wormlike micelles of cetyltrimethylammonium chloride/sodium salicylate (NaSal). We find that the breaking stress, stretch modulus, and scission energy of the charged micelles are nonmonotonic functions of oppositely charged hydrotrope (NaSal) concentration. While the stretch modulus shows a peak at a value of surfactant-to-hydrotrope concentration ratio ( R) close to unity as expected due to neutralization of head-group charge at R = 1, the breaking stress and scission energy produce a peak at R < 1.0 because of thinning of the micelle diameter with increased R. The breaking stress from the simulations depends on the rate of deformation and roughly agrees with the experimental values of Rothstein ( J. Rheol. 2003 , 47 , 1227 ) after extrapolation to the much lower experimental rates. The method and results can be used to predict the effects of flow and mechanical stress on rates of micellar breakage, which is important in the rheology of wormlike micellar solutions.

Self-Assembled Micellar Structures of Lipopeptides with Variable Number of Attached Lipid Chains Revealed by Atomistic Molecular Dynamics Simulations

We present atomistic molecular dynamics simulation study of the self-assembly behavior of toll-like agonist lipopeptides (Pam _nCSK4) in aqueous solutions. The variable number of hexadecyl lipid chains ( n = 1, 2, 3) per molecule has been experimentally suggested to have remarkable influence on their self-assembled nanostructures. Starting from preassembled spherical or bilayer configurations, the aggregates of lipopeptides, PamCSK4 and Pam₂CSK4, which contain peptide sequences CSK4 linked to either mono- or dilipid chains (Pam), evolve into spherical-like micelles within 30 ns, whereas the self-assembled structure of trilipidated lipopeptides, Pam₃CSK4, relaxes much slower and reaches an equilibrium state of flattened wormlike micelle with a bilayer packing structure. The geometric shapes and sizes, namely the gyration radii of spherical micelles and thickness of the flattened wormlike micelle, are found to be in good agreement with experimental measurements, which effectively validates the simulation models and employed force fields. Detailed analyses of molecular packing reveal that these self-assembled nanostructures all consist of a hydrophobic core constructed of lipid chains, a transitional layer, and a hydrophilic interfacial layer composed of peptide sequences. The average area per peptide head at the interfaces is found to be nearly constant for all micellar structures studied. The packing parameter of the lipopeptide molecules thus increases with the increase of the number of linked lipid chains, giving rise to the distinct micellar shape transition from spherical-like to flattened wormlike geometry with bilayer stacking, which is qualitatively different from the shape transitions of surfactant micelles induced by variation of concentration or salt type. To facilitate the close-packing of the lipid chains in the hydrophobic core, the lipopeptide molecules typically take the bent conformation with average tilt angles between the peptide sequences and the lipid chains ranging from 110° to 140°. This consequently affects the orientation angles of the lipid chains with respect to the radial or normal direction of the spherical-like or flattened wormlike micelles. In addition, the secondary structures of the peptides may also be altered by the number of lipid chains to which they are linked and the resultant micellar structures. Our simulation results on the microscopic structural features of the lipopeptide nanostructures may provide potential insights into their bioactivities and contribute to the design of bioactive medicines or drug carriers. The force fields built for these lipopeptides and the geometric packing discussions could also be adopted for simulating and understanding the self-assembly behavior of other bioactive amiphiphiles with similar chemical compositions.

Nonmonotonic Scission and Branching Free Energies as Functions of Hydrotrope Concentration for Charged Micelles

Using coarse-grained molecular dynamics simulations and an umbrella sampling method that uses local surfactant density as a reaction coordinate, we directly calculate, for the first time, both the scission and branching free energies of a model charged micelle [cationic cetyltrimethylammonium chloride (CTAC)] in the presence of inorganic and organic salts (hydrotropes). We find that while inorganic salt only weakly affects the micelle scission energy, organic hydrotropes produce a strong, nonmonotonic dependence of both scission energy and branching on salt concentration. The nonmonotonicity in scission energy is traced to a competition between electrostatic screening of the repulsions among the surfactant head groups and thinning of the micellar core, which result from attachment of the hydrotropes to the micelle surface. We are able to correlate the nonmonotonicity in the scission energy of CTAC micelles with the peak observed experimentally in viscosity versus hydrotrope concentration and the location of this peak in CTAC solutions.

Growth of wormlike micelles in nonionic surfactant solutions: Quantitative theory vs. experiment

Despite the considerable advances of molecular-thermodynamic theory of micelle growth, agreement between theory and experiment has been achieved only in isolated cases. A general theory that can provide self-consistent quantitative description of the growth of wormlike micelles in mixed surfactant solutions, including the experimentally observed high peaks in viscosity and aggregation number, is still missing. As a step toward the creation of such theory, here we consider the simplest system - nonionic wormlike surfactant micelles from polyoxyethylene alkyl ethers, C_iE_j. Our goal is to construct a molecular-thermodynamic model that is in agreement with the available experimental data. For this goal, we systematized data for the micelle mean mass aggregation number, from which the micelle growth parameter was determined at various temperatures. None of the available models can give a quantitative description of these data. We constructed a new model, which is based on theoretical expressions for the interfacial-tension, headgroup-steric and chain-conformation components of micelle free energy, along with appropriate expressions for the parameters of the model, including their temperature and curvature dependencies. Special attention was paid to the surfactant chain-conformation free energy, for which a new more general formula was derived. As a result, relatively simple theoretical expressions are obtained. All parameters that enter these expressions are known, which facilitates the theoretical modeling of micelle growth for various nonionic surfactants in excellent agreement with the experiment. The constructed model can serve as a basis that can be further upgraded to obtain quantitative description of micelle growth in more complicated systems, including binary and ternary mixtures of nonionic, ionic and zwitterionic surfactants, which determines the viscosity and stability of various formulations in personal-care and house-hold detergency.

Simulations Study of Single-Component and Mixed n-Alkyl-PEG Micelles

Here, we study one-component and mixed n-alkyl-poly(ethylene glycol) (C mE n) micelles with varying poly(ethylene glycol) (PEG) chain lengths n using coarse-grained molecular simulations. These nonionic alkyl-PEG surfactants and their aggregates are widely used in bio and chemical technology. As expected, the simulations show that increasing the PEG chain length decreases the alkyl-PEG micelle core diameter and the aggregation number but also enhances PEG chain penetration to the core region and spreads the micelle corona. Both the core and corona density are heavily dependent on the PEG chain length and decrease with increasing PEG length. Furthermore, we find that the alkyl-PEG surfactants exhibit two distinct micellization modes: surfactants with short PEG chains as their hydrophilic heads aggregate with the PEG heads relatively extended. Their aggregation number and the PEG corona density are dictated by the core carbon density. For longer PEG chains, the PEG sterics, that is, the volume occupied by the PEG head group, becomes the critical factor limiting the aggregation. Finally, simulations of binary mixtures of alkyl-PEGs of two different PEG chain lengths show that even in the absence of core-freezing, the surfactants prefer the aggregate size of their single-component solutions with the segregation propelled via enthalpic contributions. The findings, especially as they provide a handle on the density and the density profile of the aggregates, raise attention to effective packing shape as a design factor of micellar systems, for example, drug transport, solubilization, or partitioning.