Study on the transformation from linear to branched wormlike micelles: An insight from molecular dynamics simulation

Ecofriendly Viscoelastic Solutions Formed from a Recyclable Rosin-Based Amine Oxide Surfactant

Innovations in molecular structures formed using bioresources are efficient means to prepare surfactant aggregates with unique properties. Here, a rosin-based amine oxide surfactant (R-11-3-AO) containing large hydrophobic groups was synthesized from rosin derivatives, namely, dehydroabietic acid and long-chain amino acids. Cryo-transmission electron microscopy showed that R-11-3-AO molecules formed extremely long wormlike micelles with a cross-sectional diameter of 4-5 nm at a concentration of approximately 7 mmol·L^-1. A gel-like system was obtained at approximately 30 mmol·L^-1 due to the dense entanglement of the wormlike micelles. The solutions also exhibited unique shear thickening behavior at a shear rate of approximately 10 s^-1 even at high concentrations. The large hydrophobic group contained in R-11-3-AO is the origin of the strong van der Waals interactions between the surfactant molecules, resulting in the rapid growth of wormlike micelles. This rosin-based surfactant is the first recoverable amine oxide surfactant from solutions through the salting-out effect with high recovery rates. This work demonstrates the unique capabilities of rosin-based surfactants for forming wormlike micelles and provides opportunities for the development of surfactant recovery technologies.

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 rise and fall of branching: A slowing down mechanism in relaxing wormlike micellar networks

A mean-field kinetic model suggests that the relaxation dynamics of wormlike micellar networks is a long and complex process due to the problem of reducing the number of free end-caps (or dangling ends) while also reaching an equilibrium level of branching after an earlier overgrowth. The model is validated against mesoscopic molecular dynamics simulations and is based on kinetic equations accounting for scission and synthesis processes of blobs of surfactants. A long relaxation time scale is reached with both thermal quenches and small perturbations of the system. The scaling of this relaxation time is exponential with the free energy of an end cap and with the branching free energy. We argue that the subtle end-recombination dynamics might yield effects that are difficult to detect in rheology experiments, with possible underestimates of the typical time scales of viscoelastic fluids.

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.

Aging of living polymer networks: a model with patchy particles

Microrheology experiments show that viscoelastic media composed by wormlike micellar networks display complex relaxations lasting seconds even at the scale of micrometers. By mapping a model of patchy colloids with suitable mesoscopic elementary motifs to a system of worm-like micelles, we are able to simulate its relaxation dynamics, upon a thermal quench, spanning many decades, from microseconds up to tens of seconds. After mapping the model to real units and to experimental scission energies, we show that the relaxation process develops through a sequence of non-local and energetically challenging arrangements. These adjustments remove undesired structures formed as a temporary energetic solution for stabilizing the thermodynamically unstable free caps of the network. We claim that the observed scale-free nature of this stagnant process may complicate the correct quantification of experimentally relevant time scales as the Weissenberg number.

Solvent-induced morphological transitions in methacrylate-based block-copolymer aggregates

Poly(ethylene oxide)-b-poly(butylmethacrylate) (PEO-b-PBMA) copolymers have recently been identified as excellent building blocks for the synthesis of hierarchical nanoporous materials. Nevertheless, while experiments have unveiled their potential to form bicontinuous phases and vesicles, a general picture of their phase and aggregation behavior is still missing. By performing Molecular Dynamics simulations, we here apply our recent coarse-grained model of PEO-b-PBMA to investigate its self-assembly in water and tetrahydrofuran (THF) and unveil the occurrence of a wide spectrum of mesophases. In particular, we find that the morphological phase diagram of this ternary system incorporates bicontinuous and lamellar phases at high copolymer concentrations, and finite-size aggregates, such as dispersed sheets or disk-like aggregates, spherical vesicles and rod-like vesicles, at low copolymer concentrations. The morphology of these mesophases can be controlled by tuning the THF/water relative content, which has a striking effect on the kinetics of self-assembly as well as on the resulting equilibrium structures. Our results disclose the fascinating potential of PEO-b-PBMA copolymers for the templated synthesis of nanostructured materials and offer a guideline to fine-tune their properties by accurately selecting the THF/water ratio.

Replica exchange dissipative particle dynamics method on threadlike micellar aqueous solutions

The self-assembly of surfactant molecules can spontaneously result in a variety of micelle morphologies, such as spherical micelles, threadlike micelles, and vesicles, and it is therefore crucial to predict and control the self-assembly to achieve a helpful process in the fields of materials chemistry and engineering. A dissipative particle dynamics (DPD) method used in a coarse-grained molecular simulation is applied to simulate various self-assembling soft matter systems because it can handle greater length and time scales than a typical molecular dynamics simulation (MD). It should be noted that the thorough sampling of a system is not assured at low temperatures because of large complex systems with coarse-grained representations. In this article, we demonstrate that the replica exchange method (REM) is very effective for even a DPD in which the energy barrier is comparatively lower than that of a MD. A replica exchange on DPD (REDPD) simulation for threadlike micellar aqueous solutions was conducted, and the values of the potential energy and the mean aggregation number were compared. As a result, the correct values and a self-assembled structure within a low-temperature range can only be obtained through the REDPD.

A Coarse-Grained Model for Microbial Lipopeptide Surfactin and Its Application in Self-Assembly

Biosurfactants exhibit outstanding interfacial properties and unique biological activities that fairly related to their self-assembly in solutions and at interfaces. Computational simulations provide structural details of biosurfactant aggregates at the molecular level relevant to thermodynamic properties, but the understanding of kinetics of self-assembly remains limited due to lower simulation efficiency. In this work, a coarse-grained model has been developed for microbial lipopeptide surfactin, and surfactin monolayer at the octane/water interface and micelle in aqueous solution were studied using molecular dynamics simulations. Interaction parameters were optimized and validated by comparing with results obtained from experiments and atomistic molecular dynamics simulations. In particular, self-assembly of surfactin in aqueous solution was studied using the optimized parameters. Results showed that coarse-grained simulations well reproduced structural properties of surfactin monolayer and micelle and the molecular behavior such as surfactin orientation and conformation. Self-assembly features of surfactin in different stages have been captured, and the aggregation numbers of dominant clusters were in accordance with experimental data. This report suggested that the present coarse-grained model and interaction parameters allowed surfactin simulations over longer timescales and larger systems, which provide insights into characterizing both the kinetics of surfactin self-assembly and the adsorption of surfactin onto varying interfaces.

Different strategies of foam stabilization in the use of foam as a fracturing fluid

An attractive alternative to mitigate the adverse effects of conventional water-based fluids on the efficiency of hydraulic fracturing is to inject foam-based fracking fluids into reservoirs. The efficiency of foaming fluids in subsurface applications largely depends on the stability and transportation of foam bubbles in harsh environments with high temperature, pressure and salinity, all of which inevitably lead to poor foam properties and thus limit fracturing efficiency. The aim of this paper is to elaborate popular strategies of foam stabilization under reservoir conditions. Specifically, this review first discusses three major mechanisms governing foam decay and summarizes recent progress in research on these phenomena. Since surfactants, polymers, nanoparticles and their composites are popular options for foam stabilization, their stabilizing effects, especially the synergies in composites, are also reviewed. In addition to reporting experimental results, the paper also reports recent advances in interfacial properties via molecular dynamical simulation, which provide new insights into gas/liquid interfacial properties under the influence of surfactants at molecular scale. The results of both experiments and simulations indicate that foam additives play an essential role in foam stability and the synergic effects of surfactants and nanoparticles exhibit more favorable performance.

Interfacial Interaction Enhanced Rheological Behavior in PAM/CTAC/Salt Aqueous Solution-A Coarse-Grained Molecular Dynamics Study

Interfacial interactions within a multi-phase polymer solution play critical roles in processing control and mass transportation in chemical engineering. However, the understandings of these roles remain unexplored due to the complexity of the system. In this study, we used an efficient analytical method-a nonequilibrium molecular dynamics (NEMD) simulation-to unveil the molecular interactions and rheology of a multiphase solution containing cetyltrimethyl ammonium chloride (CTAC), polyacrylamide (PAM), and sodium salicylate (NaSal). The associated macroscopic rheological characteristics and shear viscosity of the polymer/surfactant solution were investigated, where the computational results agreed well with the experimental data. The relation between the characteristic time and shear rate was consistent with the power law. By simulating the shear viscosity of the polymer/surfactant solution, we found that the phase transition of micelles within the mixture led to a non-monotonic increase in the viscosity of the mixed solution with the increase in concentration of CTAC or PAM. We expect this optimized molecular dynamic approach to advance the current understanding on chemical-physical interactions within polymer/surfactant mixtures at the molecular level and enable emerging engineering solutions.

Protobiotic Systems Chemistry Analyzed by Molecular Dynamics

Systems chemistry has been a key component of origin of life research, invoking models of life's inception based on evolving molecular networks. One such model is the graded autocatalysis replication domain (GARD) formalism embodied in a lipid world scenario, which offers rigorous computer simulation based on defined chemical kinetics equations. GARD suggests that the first pre-RNA life-like entities could have been homeostatically-growing assemblies of amphiphiles, undergoing compositional replication and mutations, as well as rudimentary selection and evolution. Recent progress in molecular dynamics has provided an experimental tool to study complex biological phenomena such as protein folding, ligand-receptor interactions, and micellar formation, growth, and fission. The detailed molecular definition of GARD and its inter-molecular catalytic interactions make it highly compatible with molecular dynamics analyses. We present a roadmap for simulating GARD's kinetic and thermodynamic behavior using various molecular dynamics methodologies. We review different approaches for testing the validity of the GARD model by following micellar accretion and fission events and examining compositional changes over time. Near-future computational advances could provide empirical delineation for further system complexification, from simple compositional non-covalent assemblies towards more life-like protocellular entities with covalent chemistry that underlies metabolism and genetic encoding.

Cooperative assembly of Janus particles and amphiphilic oligomers: the role of Janus balance

Cooperative assembly of nanoparticles and amphiphiles has emerged as a significant strategy for constructing hybrid nanocomposites with desired architectures and properties. It is of great significance to develop novel hybrid nanostructures with controlled spatial localization of nanoparticles within hybrid assemblies. Here, by adopting dissipative particle dynamics simulations, the cooperative assembly of Janus particles and amphiphilic oligomers is studied. We demonstrate that a variety of defined hybrid nanostructures such as balls, sticks, disks, lines, vesicles, and networks can be achieved by the cooperative assembly of Janus particles and amphiphilic oligomers. Furthermore, the investigation of the kinetic pathway illustrates that the formation of hybrid assemblies is an entropy-driven process. Our simulation results suggest that the Janus balance of nanoparticles can significantly affect the structure and size of hybrid aggregates and the spatial localization of Janus particles within hybrid assemblies. These findings not only enrich our understanding of the cooperative assembly of Janus nanoparticles and amphiphiles, but also offer a feasible strategy to prepare hybrid materials with controlled localization of nanoparticles.

Thermodynamic insights and molecular environments into catanionic surfactant systems: Influence of chain length and molar ratio

HYPOTHESIS

Imidazolium-based Ionic liquids as new generation cationic surfactants can provide designable alkyl chain length. In the catanionic surfactant systems, the alkyl chain lengths and molar ratios can greatly influence the interactions such as electrostatic and hydrophobic interaction. The variation in these interactions has a significant effect on the molecular environments of the self-assembly structure, and this process is always accompanied by the transition of aggregates and release or consumption of heat. Hence, it is of interest to study the relationship between intermolecular interactions, molecular environments, self-assembly structure and the change in energy of system in the catanionic surfactant mixed systems.

EXPERIMENTS

The enthalpy change ΔH of titrations the imidazolium-based into SDS micelle solution was studied to characterize the heat by using isothermal titration calorimetry (ITC) during the transitions of the aggregate structures. The corresponding self-assembly structure was monitored via cryogenic transmission electron microscopy (cryo-TEM). Employing proton magnetic resonance (¹H NMR), we focus on the interactions between imidazolium-based ILs and SDS based on the variations in the molecular environments of aggregates.

FINDINGS

Of these imidazolium-based ionic liquids/SDS system, the 1-octyl-3-methylimidazolium ([OMIM]Cl)/SDS system shows several features such as intense energy absorption and releasing processes, which indicate the formation of high entanglement wormlike micelles and vesicles. This is related to the formation of self-adjusting state between the SDS and [OMIM]Cl molecules due to the balance between the electrostatic interaction and hydrophobic interaction. Varying the alkyl chain length appears to cause significant differences to the molecular environments. From the molecular environments, three different models about the polarity of the catanionic surfactant molecules are used to explain the balance of the intermolecular interactions.

Enhanced stability of crosslinked and charged unimolecular micelles from multigeometry triblock copolymers with short hydrophilic segments: dissipative particle dynamics simulation

High micellar stability and well-performed drug loading and release are two conflicting factors for unimolecular micelles as an ideal drug delivery system. Achieving the formation of unimolecular micelles with short hydrophilic blocks is a challenging and promising approach to solve this bottleneck and limitation of current unimolecular micelle systems. In this work, dissipative particle dynamics (DPD) simulation is used to study the synergetic effect of crosslinking and electrostatic repulsion on stability of unimolecular micelles and to analyze the micro-mechanism and factors influencing this synergetic stabilization strategy. The strategy can generate unimolecular micelles with extremely high stability for various supramolecular polymers with short hydrophilic chains. Protonation of DEAEMA blocks leads to a large improvement in micellar hydrophilicity. The protonated middle layer further shrinks through crosslinking to produce the largest charge density, enlarging the electrostatic repulsion between colloidal particles. Additionally, the crosslinking and protonation treatment maximizes the extension degree of hydrophilic EO segments due to the increasing steric hindrance and poor compatibility between DEAHEMA and EO blocks. In this study, the relation between shrinkage degree of hydrophobic cores and stability of unimolecular micelles is first reported. The above-mentioned transition of micellar structures and properties results in the maximum degree of core shrinkage (R_g of MMA blocks) corresponding to the high stability of unimolecular micelles. Further study shows that the increasing cyclization degree, the mode of end cyclization, and the crosslinking and electrostatic repulsion of the middle layer all exert favorable effects on the stability of unimolecular micelles due to controlled shrinkage of hydrophobic cores.

Investigation of Active-Inactive Material Interdigitated Aggregates Formed by Wormlike Micelles and Cellulose Nanofiber

In this work, a novel active-inactive material interdigitated aggregates (AIMIAs) structure was constructed by self-assembled wormlike micelles (WLMs) and one-dimensional cellulose nanofiber (CNF). The rheological behaviors and microstructures of the AIMIA systems with different CNF concentrations were investigated by rheometer, cryogenic transmission electron microscopy, and environmental scanning electron microscope. Some key parameters, including zero-shear viscosity (η₀), relaxing time (τ_R), and contour length ( L), were calculated to analyze the changes in the properties of different systems. Meanwhile, a proper mechanism describing the interaction between CNF and WLMs was proposed. Through this work, we expect to deepen the understanding of the AIMIAs structure and widen its application.

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.

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.

Prediction of striped cylindrical micelles (SCMs) formed by dodecyl-β-d-maltoside (DDM) surfactants

Using fully atomistic and coarse-grained (CG) molecular dynamics (MD) simulations, we report, for the first time, the self-assembly of initially randomly dispersed dodecyl-β-d-maltoside (DDM) surfactants into a striped cylindrical micelle (SCM) with lamellae of surfactant heads and tails alternating along the cylindrical axis, with both heads and tails in contact with the water. By changing the interaction strength of the head group with water relative to itself, we find that such micelles are most likely for head groups with marginal solubility in the water solvent. Unlike the surfactants in a regular cylindrical micelle, whose tails are in the fluid micelle interior, the diffusion of DDM surfactants along the micelle body is blocked by the lamellar patterning. As a consequence, branches cannot slide along the micelle body and surfactant molecules cannot exchange between the micelle body and the branch, which should have a significant impact on the rheological properties of these micelles.