Sphere-to-rod transitions of nonionic surfactant micelles in aqueous solution modeled by molecular dynamics simulations

Self-assembly and structures of nanoscale double emulsion droplets through coarse-grained molecular dynamics simulations

Examples of self-assembled multiple emulsion droplets on the nanometre scale are very rare. In this work, we use coarse-grained (CG) molecular dynamics simulations to study the self-assembly of ternary mixtures consisting of water, n-heptane, and nonionic surfactant tetraethylene glycol monododecyl ether (C12E4). The water volume fractions studied are 1%, 3%, and 5%, respectively. Various nanoscale emulsions are obtained in a spontaneous process. When the water/surfactant volume ratio vm/s = 1.0/1.0, the obtained emulsion droplets are identified as oil-in-water-in-oil (O/W/O) double types, consisting of an oil core, an inner surfactant layer, a water layer, and an outer surfactant layer. The water molecules are distributed around the hydrophilic ends of the surfactants, while the hydrophobic ends of the surfactants wrap the oil cores and penetrate into the oil bulk. Hydrogen-bond interactions among water and the hydrophilic ends of the surfactants form cross-links that stabilize the double emulsion droplets. The sizes of all the oil cores inside the droplets are <6 nm in diameter, even with the highest water volume fraction of 5%. Both the concentration of free water molecules on the order of 10-6 mol/cm3 and the favourable energy change during emulsion formation indicate that the emulsion droplets are thermodynamically stable. In contrast, for vm/s = 1.0/5.5, no double emulsion but a simple water-in-oil emulsion was observed, with morphologies evolving from oblate to bicontinuous phases with an increase in the water volume fraction from 1% to 5%. Our coarse-grained molecular dynamics simulations provide valuable insight for the preparation of nanoscale double emulsions and the characterization of their structures.

Micelle kinetics of photoswitchable surfactants: Self-assembly pathways and relaxation mechanisms

HYPOTHESIS

A key question in the kinetics of surfactant self-assembly is whether exchange of unimers or fusion/fission of entire micelles is the dominant pathway. In this study, an isomerizable surfactant is used to explore fundamental out-of-equilibrium kinetics and mechanisms for growth and dissolution of micelles.

EXPERIMENTS

The kinetics of cationic surfactant 4-butyl-4'-(3-trimethylammoniumpropoxy)-phenylazobenzene was studied using molecular dynamics simulations. The fusion and exchange processes were investigated using umbrella sampling. Equilibrium states were validated by comparison with small-angle X-ray scattering data. The photo-isomerization event was simulated by modifying the torsion potential of the photo-responsive group to emulate the trans-to-cis transition.

FINDINGS

Micelle growth is dominated by unimer exchange processes, whereas, depending on the conditions, dissolution can occur both through fission and unimer expulsion. Fusion barriers increase steeply with the aggregation number making this an unlikely pathway to equilibrium for micelles of sizes that fit with the experimental data. The barriers for unimer expulsion remain constant and are much lower for unimer insertion, making exchange more likely at high aggregation. When simulating photo-conversion events, both fission and a large degree of unimer expulsion can occur depending on the extent of the out-of-equilibrium stress that is put on the system.

Phase Behavior of Alkyl Ethoxylate Surfactants in a Dissipative Particle Dynamics Model

We present a dissipative particle dynamics (DPD) model capable of capturing the liquid state phase behavior of nonionic surfactants from the alkyl ethoxylate (C_nE_m) family. The model is based upon our recent work [Anderson et al. J. Chem. Phys. 2017, 147, 094503] but adopts tighter control of the molecular structure by setting the bond angles with guidance from molecular dynamics simulations. Changes to the geometry of the surfactants were shown to have little effect on the predicted micelle properties of sampled surfactants, or the water-octanol partition coefficients of small molecules, when compared to the original work. With these modifications the model is capable of reproducing the binary water-surfactant phase behavior of nine surfactants (C₈E₄, C₈E₅, C₈E₆, C₁₀E₄, C₁₀E₆, C₁₀E₈, C₁₂E₆, C₁₂E₈, and C₁₂E₁₂) with a good degree of accuracy.

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.

Coarse-Grained Molecular Dynamics Simulations of Paclitaxel-Loaded Polymeric Micelles

A continuous manufacturing technology based on coaxial turbulent jet in coflow was previously developed to produce paclitaxel-loaded polymeric micelles. Herein, coarse-grained molecular dynamics (CG-MD) simulations were implemented to better understand the effect of the material attributes (i.e., the drug-polymer ratio and the ethanol concentration) and process parameters (i.e., temperature) on the self-assembly process of polymeric micelles as well as to provide molecular details on micelle instability. An all-atom (AA) poly (ethylene glycol)-poly (lactic acid) (PEG-PLA) polymer model was developed as the reference for parameterizing a coarse-grained (CG) model, and the AA polymer model was further validated with experimental glass transition temperature (T_g). The model transferability was verified by comparing structural properties between the AA and CG models. The CG model was further validated with experimental data, including micelle particle size measurements and drug encapsulation efficiency. Furthermore, the encapsulation of paclitaxel into the polymeric micelles was included in the simulations, taking into consideration the interactions between the paclitaxel and the polymers. The results from various points of view demonstrated a strong dependence of the shape of the micelles on the drug encapsulation, with micelles transitioning from spherical to ellipsoidal structures with an increasing paclitaxel amount. Simulation data were also used to identify the critical aggregation number (i.e., the number of polymer and drug molecules required for transition from one shape to another). Improved micellar structural stability was found with a larger micellar size and less solvent accessibility. Lastly, an evaluation was performed on the micellar dissociation free energy using a steered molecular dynamics simulation over a range of temperatures and ethanol concentrations. These simulations revealed that at higher ethanol and temperature conditions, micelles become destabilized, resulting in greater paclitaxel release. The increased drug release was determined to originate from the solvation of the hydrophobic core, which promoted micellar swelling and an associated reduction in hydrophobic interactions, leading to a loosely packed micellar structure.

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.

Structural Insights into Self-Assembled Aerosol-OT Aggregates in Aqueous Media Using Atomistic Molecular Dynamics

In water, the surfactant dioctyl sulfosuccinate (Aerosol-OT or AOT) exhibits diverse aggregate structures, ranging from micelles to lamella. An atomic-level understanding, however, of the formation and structure of these aggregates is lacking. Herein, using atomistic molecular dynamics (MD) with microsecond-long simulations, self-assembly of AOT in water is studied for concentrations of 1, 7.2, and 20 wt % at 293 K and for 7.2 wt % at 353 K. Assembly proceeds through stepwise association and dissociation of single AOT molecules, and the fusion and fission of AOT clusters. At 293 K, AOT self-assembles into either (i) spherical micelles (1 wt %), (ii) biphasic systems consisting of rod-like and prolate spheroidal micelles (7.2 wt %), or (iii) bilayers (20 wt %). We hypothesize that the observed rod-like structure is a precursor to lamellar microdomains found experimentally in biphasic dispersions. Increasing temperature to 353 K at 7.2 wt % results in a system consisting of prolate micelles but no rod-like micelles. Simulated phase behavior agrees with previously published experimental observations. Individual aggregates formed during self-assembly are identified using graph theory. Structural metrics of these aggregates like the radius of gyration, shape anisotropy, and prolateness are presented. Trends in structural metrics quantitatively reflect how shapes and sizes of AOT aggregates vary with surfactant concentration and temperature. These simulations provide deeper insight into open questions in the scientific community and demonstrate a method to generate physics-based micelle structures that can be used to rationalize experimental observations.

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.

The Martini Model in Materials Science

The Martini model, a coarse-grained force field initially developed with biomolecular simulations in mind, has found an increasing number of applications in the field of soft materials science. The model's underlying building block principle does not pose restrictions on its application beyond biomolecular systems. Here, the main applications to date of the Martini model in materials science are highlighted, and a perspective for the future developments in this field is given, particularly in light of recent developments such as the new version of the model, Martini 3.

Unveiling the phase behavior of CiEj non-ionic surfactants in water through coarse-grained molecular dynamics simulations

Poly(oxyethylene) alkyl ethers, usually denoted by CiEj surfactants, exhibit a rich phase behavior in water, self-assembling to form a variety of 3-D structures with a controllable morphology that find multiple applications across different industrial segments. Hence, being able to describe and understand the effect of molecular structure on the phase behavior of these systems is highly relevant for the efficient design of new materials and their applications. Considering the promising results obtained over the last decade using the MARTINI model to describe ethylene-oxide containing compounds, an extensive assessment of the ability of such a model to describe the phase behavior of CiEj in water was carried out and results are presented here. Given the overall poor temperature transferability of the MARTINI model, mostly due to the lack of an accurate representation of hydrogen bonding, simulations were carried out at a single temperature of 333 K, where most phases are expected to occur according to experiments. Different chain lengths of both the hydrophobic and hydrophilic moieties, spanning a wide range of hydrophilic-lipophilic balance values, were investigated and the phase diagrams of various CiEj surfactants explored over a wide concentration range. The model was able to satisfactorily describe the effect of surfactant structure and concentration on mesophase formation. The stability and dimensions of the obtained phases, and the prediction of some unique features such as the characterization of a singular lamellar phase are presented. The results obtained in this work highlight both the predictive ability and the transferability of the MARTINI forcefield in the description of such systems. Moreover, the model was shown to provide adequate descriptions of the micellar phase in terms of micelle dimensions, critical micelle concentration, and average aggregation number.

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 (N_agg) 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 (R_g), 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 R_g and SASA values indicated the notable intercalation of PG in the CTAB micelles.

Aqueous phase behavior of the PEO-containing non-ionic surfactant C12E6: A molecular dynamics simulation study

HYPOTHESIS

Non-ionic surfactants containing polyethylene oxide (PEO) chains are widely used in drug formulations, cosmetics, paints, textiles and detergents. High quality molecular dynamics models for PEO surfactants can give us detailed, atomic-scale information about the behavior of surfactant/water mixtures.

SIMULATIONS

We used two molecular dynamics force fields (FFs), 2016H66 and 53A6_DBW, to model the simple non-ionic PEO surfactant, hexaoxyethylene dodecyl ether (C₁₂E₆). We investigated surfactant/water mixtures that span the phase diagram of starting from randomly distributed arrangements. In some cases, we also started with prebuilt, approximate models. The simulations results were compared with the experimentally observed phase behavior.

FINDINGS

Overall, this study shows that the spontaneous self-assembly of PEO non-ionic surfactants into different colloidal structures can be accurately modeled with MD simulations using the 2016H66 FF although transitions to well-formed hexagonal phase are slow. Of the two FFs investigated, the 2016H66 FF better reproduces the experimental phase behavior across all regions of the C₁₂E_6/water phase diagram.

Hydrophilic dangling chain interfacial segregation in polyurethane networks at aqueous interfaces and its underlying mechanisms: molecular dynamics simulations

Polymer networks with hydrophilic dangling chains are ideal candidates for many submerged applications, e.g., protein non-adhesive coatings with non-fouling behavior. The dangling chains segregate from the polymer network towards the water and form a brush-like structure at the interface. Several factors such as the polymer network structure, dangling chain length, and water/dangling chain interaction may all affect the interfacial performance of the polymer. Therefore, we employed a Martini based coarse-grained (CG) molecular dynamics (MD) simulation to elucidate the influences of the abovementioned parameters on dangling chain interfacial segregation. We built up several polyurethane (PU) networks based on poly(tetra methylene glycol) (PTMG), as a macrodiol, and methoxy poly(ethylene glycol) (mPEG), as a dangling chain, with varying molecular weights. We found out that the macrodiol/dangling chain length ratio considerably smaller than one impedes the migration of dangling chains towards the water interface, while the dangling chain hydrophilicity and length determine the polymer interfacial layer density/thickness. Then, we artificially changed the dangling chain affinity to water from an intermediate to a very attractive water/dangling chain interaction. We justified that a brush-like structure forms in two consecutive steps: first, a longitudinal, and then a lateral migration of dangling chains in water. The latter step results in a uniform interfacial layer over the polymer interface that mainly occurs in the case of the attractive water/dangling chain interaction.

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.

Molecular Dynamics Simulations of Aqueous Nonionic Surfactants on a Carbonate Surface

The interactions and structure of secondary alcohol ethoxylates with 15 and 40 ethoxylate units in water near a calcite surface are studied. It is found that water binds preferentially to the calcite surface. Prediction of the free-energy landscape for surfactant molecules shows that single-surfactant molecules do not adsorb because they cannot get close enough to the surface because of the water layer for attractive ethoxylate-calcite or dispersion interactions to be significant. Micelles can adsorb onto the surface even with the intervening water layer because of the integrative effect of the attractive interactions of all the surfactant molecules. Adsorption is found to increase because of the closer proximity of the micelles to the surface due to a weakened water layer at higher temperatures. The free-energy well and barrier values are used to estimate surface to bulk partition coefficients for different surfactants and temperatures, and qualitative agreement is found with experimental observations. The combined effect of surfactant-water and surfactant-solid interactions is found to be responsible for an increased adsorption for nonionic surfactants as the system approaches the cloud point.

Can Polarity-Inverted Surfactants Self-Assemble in Nonpolar Solvents?

We investigate the self-assembly process of a surfactant with inverted polarity in water and cyclohexane using both all-atom and coarse-grained hybrid particle-field molecular dynamics simulations. Unlike conventional surfactants, the molecule under study, proposed in a recent experiment, is formed by a rigid and compact hydrophobic adamantane moiety, and a long and floppy triethylene glycol tail. In water, we report the formation of stable inverted micelles with the adamantane heads grouping together into a hydrophobic core and the tails forming hydrogen bonds with water. By contrast, microsecond simulations do not provide evidence of stable micelle formation in cyclohexane. Validating the computational results by comparison with experimental diffusion constant and small-angle X-ray scattering intensity, we show that at laboratory thermodynamic conditions the mixture resides in the supercritical region of the phase diagram, where aggregated and free surfactant states coexist in solution. Our simulations also provide indications as to how to escape this region to produce thermodynamically stable micellar aggregates.

A Surface Site Interaction Point Method for Dissipative Particle Dynamics Parametrization: Application to Alkyl Ethoxylate Surfactant Self-Assembly

Dissipative particle dynamics (DPD) is a coarse-grained approach to the simulation of large supramolecular systems, but one limitation has been that the parameters required to describe the noncovalent interactions between beads are not readily accessible. A first-principles computational method has been developed so that bead interaction parameters can be calculated directly from ab initio gas-phase molecular electrostatic potential surfaces of the molecular fragments that represent the beads. A footprinting algorithm converts the molecular electrostatic potential surfaces into a discrete set of surface site interaction points (SSIPs), and these SSIPs are used in the SSIMPLE (surface site interaction model for the properties of liquids at equilibrium) algorithm to calculate the free energies of transfer of one bead into a solution of any other bead. The bead transfer free energies are then converted into the required DPD interaction parameters for all pairwise combinations of different beads. The reliability of the parameters was demonstrated using DPD simulations of a range of alkyl ethoxylate surfactants. The simulations reproduce the experimentally determined values of the critical micelle concentration and mean aggregation number well for all 22 surfactants studied.

Applications of Micelles in Catalyzing Organic Reactions

:

Micellar chemistry is gaining considerable interest among organic chemists because these reactions are carried out in environmentally benign solvents like water. Owing to the exhaustive use of toxic solvents in carrying out the different chemical reactions, there is a pressing need for alternative approaches either environmental friendly or having minimum impact on the environment. In this article, we aim to discuss the various aspects of micellar chemistry viz-a-viz its role in guiding the chemical reactions. Micelles help to drive various kinds of organic reactions including oxidations, reductions, carbon-carbon bond formation, carbon-heteroatom bond formation, multi-component reactions, Pd-coupling reaction, olefin metathesis reaction, Morita-Baylis-Hillman reaction, etc. in water.

Self-Assembly of Phosphocholine Derivatives Using the ELBA Coarse-Grained Model: Micelles, Bicelles, and Reverse Micelles

The ELBA coarse-grained force field was originally developed for lipids, and its water model is described as a single-site Lennard-Jones particle with electrostatics modeled by an embedded point-dipole, while other molecules in this force field have a three (or four)-to-one mapping scheme. Here, ELBA was applied to investigate the self-assembly processes of dodecyl-phosphocholine (DPC) micelle, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1,2-dihexaoyl-sn-glycero-3-phosphocholine (DPPC/DHPC) bicelles, and DPPC/cyclohexane/water reverse micelles through coarse-grained molecular dynamics (MD) simulations. New parameters were obtained using a simplex algorithm-based calibration procedure to determine the Lennard-Jones parameters for cyclohexane, dodecane, and cyclohexane-dodecane cross-interactions. Density, self-diffusion coefficient, surface tension, and mixture excess volume were found to be in fair agreement with experimental data. These new parameters were used in the simulations, and the obtained structures were analyzed for shape, size, volume, and surface area. Except for the shape of DPC micelles, all other properties match well with available experimental data and all-atom simulations. Remarkably, in agreement with experiments the rodlike shape of the DPPC reverse micelle is well described by ELBA, while all-atom data in the literature predicts a disclike shape. To further check the consistency of the force field in reproducing the correct shapes of reverse micelles, additional simulations were performed doubling the system size. Two distinct reverse micelles were obtained both presenting the rodlike shape and correct aggregation number.

Effect of the Anomeric Configuration on the Micellization of Hexadecylmaltoside Surfactants

The self-assembly of the two anomeric forms of n-hexadecyl-d-maltopyranoside (denoted α-C₁₆G₂ and β-C₁₆G₂) has been studied in dilute aqueous solution by means of surface tension measurements, scattering methods (dynamic light scattering, static light scattering, and small-angle X-ray and neutron scattering), and cryo-transmission electron microscopy at different surfactant concentrations and temperatures. Surface tension measurements demonstrate differences in the surfactant adsorption at the air-water interface, where α-C₁₆G₂ shows a lower CMC than β-C₁₆G₂. Similarly, micelle morphology was found to profoundly depend on anomerism. β-C₁₆G₂ preferentially forms very elongated micelles with large persistence lengths, whereas α-C₁₆G₂ assembles into smaller micelles for which the structure varies with concentration and temperature. The differences between the two surfactant anomers in terms of self-assembly can be attributed to the interaction between neighboring headgroups. Specifically, β-C₁₆G₂ allows for a closer packing in the palisade layer, hence reducing the micelle curvature and promoting the formation of more elongated micelles. Strong intermolecular headgroup interactions may also account for the observed rigidity of the micelles.

Structure and Dynamics of Spherical and Rodlike Alkyl Ethoxylate Surfactant Micelles Investigated Using NMR Relaxation and Atomistic Molecular Dynamics Simulations

Predicting and controlling the properties of amphiphile aggregate mixtures require understanding the arrangements and dynamics of the constituent molecules. To explore these topics, we study molecular arrangements and dynamics in alkyl ethoxylate nonionic surfactant micelles by combining NMR relaxation measurements with large-scale atomistic molecular dynamics simulations. We calculate parameters that determine relaxation rates directly from simulated trajectories, without introducing specific functional forms to describe the dynamics. NMR relaxation rates, which depend on relative motions of interacting atom pairs, are influenced by wide distributions of dynamic time scales. We find that relative motions of neighboring atom pairs are rapid and liquidlike but are subject to structural constraints imposed by micelle morphology. Relative motions of distant atom pairs are slower than nearby atom pairs because changes in distances and angles are smaller when the moving atoms are further apart. Large numbers of atom pairs undergoing these slow relative motions contribute to predominantly negative cross-relaxation rates. For spherical micelles, but not for cylindrical micelles, cross-relaxation rates are positive only for surfactant tail atoms connected to the hydrophilic headgroup. This effect is related to the lower packing density of these atoms at the hydrophilic-hydrophobic boundary in spherical vs cylindrical arrangements, with correspondingly rapid and less constrained motion of atoms at the boundary.

Simulation of diblock copolymer surfactants. I. Micelle free energies

Semigrand hybrid Monte Carlo simulations are used to measure equilibrium properties of micelles formed in a simple bead-spring model of asymmetric A-B diblock copolymer surfactant molecules in an A homopolymer solvent, over a range of values of surfactant solubility. Simulations are used to accurately measure the free energy of formation of micellar clusters as a function of aggregation number over a wide range of values, and to characterize the crossover from spherical to rodlike micelle shape with increasing aggregation number. Dynamical properties of the same model are considered in an accompanying paper [Phys. Rev. E 100, 012603 (2019)10.1103/PhysRevE.100.012603].

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.

Polymer Amphiphiles for Photoregulated Anticancer Drug Delivery

We report the synthesis of amphiphilic polymers featuring lipophilic stearyl chains and hydrophilic poly(ethylene glycol) polymers that are connected through singlet oxygen-cleavable alkoxyanthracene linkers. These amphiphilic polymers assembled in water to form micelles with diameters of ∼20 nm. Reaction of the alkoxyanthracene linkers with light and O₂ cleaved the ether C-O bonds, resulting in formation of the corresponding 9,10-anthraquinone derivatives and concomitant disruption of the micelles. These micelles were loaded with the chemotherapeutic agent doxorubicin, which was efficiently released upon photo-oxidation. The drug-loaded reactive micelles were effective at killing cancer cells in vitro upon irradiation at 365 nm, functioning through both doxorubicin release and photodynamic mechanisms.

The key to controlling the morphologies of quantum nanocrystals: spherical carborane ligands

Carboranyl ligands with different functional moieties have been used to obtain different morphologies in the quantum regime using a colloidal route.

How Hydrophobic Hydration Destabilizes Surfactant Micelles at Low Temperature: A Coarse-Grained Simulation Study

Micelles are self-assembled aggregates of amphiphilic surfactant molecules that are important in a variety of applications, including drug delivery, detergency, and catalysis. It is known that the micellization process is driven by the same physiochemical forces that drive protein folding, aggregation, and biological membrane self-assembly. Nevertheless, the molecular details of how micelle stability changes in water at low temperature are not fully clear. We develop and use a coarse-grained model to investigate how the interplay between nonionic surfactants and the surrounding water at the nanoscale affects the stability of micelles at high and low temperatures. Simulations of preformed C₁₂E₅ micelles in explicit water at a range of temperatures reveal the existence of two distinct surfactant conformations within the micelle, a bent structure and an extended structure, the latter being more prevalent at low temperature. Favorable interactions of the surfactant with more ordered solvation water stabilizes the extended configuration, allowing nanoscale wetting of the dry, hydrophobic core of the micelle, leading to the micelle breaking. Taken together, our coarse-grained simulations unravel how energetic and structural changes of the surfactant and the surrounding water destabilize micelles at low temperature, which is a direct consequence of the weakened hydrophobicity. Our approach thus provides an effective mean for extracting the molecular-level changes during hydrophobicity-driven destabilization of molecular self-assembly, which is important in a wide range of fields, including biology, polymer science, and nanotechnology.

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.

Transferable MARTINI Model of Poly(ethylene Oxide)

Motivated by the deficiencies of the previous MARTINI models of poly(ethylene oxide) (PEO), we present a new model featuring a high degree of transferability. The model is parametrized on (a) a set of 8 free energies of transfer of dimethoxyethane (PEO dimer) from water to solvents of varying polarity; (b) the radius of gyration in water at high dilution; and (c) matching angle and dihedral distributions from atomistic simulations. We demonstrate that our model behaves well in five different areas of application: (1) it produces accurate densities and phase behavior or small PEO oligomers and water mixtures; (2) it yields chain dimensions in good agreement with the experiment in three different solvents (water, diglyme, and benzene) over a broad range of molecular weights (∼1.2 kg/mol to 21 kg/mol); (3) it reproduces qualitatively the structural features of lipid bilayers containing PEGylated lipids in the brush and mushroom regime; (4) it is able to reproduce the phase behavior of several PEO-based nonionic surfactants in water; and (5) it can be combined with the existing MARTINI PS to model PS-PEO block copolymers. Overall, the new PEO model outperforms previous models and features a high degree of transferability.

Modeling micelle formation and interfacial properties with iSAFT classical density functional theory

Surfactants reduce the interfacial tension between phases, making them an important additive in a number of industrial and commercial applications from enhanced oil recovery to personal care products (e.g., shampoo and detergents). To help obtain a better understanding of the dependence of surfactantproperties on molecular structure, a classical density functional theory, also known as interfacial statistical associating fluid theory, has been applied to study the effects of surfactant architecture on micelle formation and interfacial properties for model nonionic surfactant/water/oil systems. In this approach, hydrogen bonding is explicitly included. To minimize the free energy, the system minimizes interactions between hydrophobic components and hydrophilic components with water molecules hydrating the surfactant head group. The theory predicts micellar structure, effects of surfactant architecture on critical micelle concentration, aggregation number, and interfacial tension isotherm of surfactant/water systems in qualitative agreement with experimental data. Furthermore, this model is applied to study swollen micelles and reverse swollen micelles that are necessary to understand the formation of a middle-phase microemulsion.

Engineering single-polymer micelle shape using nonuniform spontaneous surface curvature

Conventional micelles, composed of simple amphiphiles, exhibit only a few standard morphologies, each characterized by its mean surface curvature set by the amphiphiles. Here we demonstrate a rational design scheme to construct micelles of more general shape from polymeric amphiphiles. We replace the many amphiphiles of a conventional micelle by a single flexible, linear, block copolymer chain containing two incompatible species arranged in multiple alternating segments. With suitable segment lengths, the chain exhibits a condensed spherical configuration in solution, similar to conventional micelles. Our design scheme posits that further shapes are attained by altering the segment lengths. As a first study of the power of this scheme, we demonstrate the capacity to produce long-lived micelles of horseshoe form using conventional bead-spring simulations in two dimensions. Modest changes in the segment lengths produce smooth changes in the micelle's shape and stability.

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.