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

Molecularly informed field theory for estimating critical micelle concentrations of intrinsically disordered protein surfactants

The critical micelle concentration (CMC) is a crucial parameter in understanding the self-assembly behavior of surfactants. In this study, we combine simulation and experiment to demonstrate the predictive capability of molecularly informed field theories in estimating the CMC of biologically based protein surfactants. Our simulation approach combines the relative entropy coarse-graining of small-scale atomistic simulations with large-scale field-theoretic simulations, allowing us to efficiently compute the free energy of micelle formation necessary for the CMC calculation while preserving chemistry-specific information about the underlying surfactant building blocks. We apply this methodology to a unique intrinsically disordered protein platform capable of a wide variety of tailored sequences that enable tunable micelle self-assembly. The computational predictions of the CMC closely match experimental measurements, demonstrating the potential of molecularly informed field theories as a valuable tool to investigate self-assembly in bio-based macromolecules systematically.

A molecularly informed field-theoretic study of the complexation of polycation PDADMA with mixed micelles of sodium dodecyl sulfate and ethoxylated surfactants

The self-assembly and phase separation of mixtures of polyelectrolytes and surfactants are important to a range of applications, from formulating personal care products to drug encapsulation. In contrast to systems of oppositely charged polyelectrolytes, in polyelectrolyte-surfactant systems the surfactants micellize into structures that are highly responsive to solution conditions. In this work, we examine how the morphology of micelles and degree of polyelectrolyte adsorption dynamically change upon varying the mixing ratio of charged and neutral surfactants. Specifically, we consider a solution of the cationic polyelectrolyte polydiallyldimethylammonium, anionic surfactant sodium dodecyl sulfate, neutral ethoxylated surfactants (C[Formula: see text]EO[Formula: see text]), sodium chloride salt, and water. To capture the chemical specificity of these species, we leverage recent developments in constructing molecularly informed field theories via coarse-graining from all-atom simulations. Our results show how changing the surfactant mixing ratios and the identity of the nonionic surfactant modulates micelle size and surface charge, and as a result dictates the degree of polyelectrolyte adsorption. These results are in semi-quantitative agreement with experimental observations on the same system.

Acid Radical Tolerance of Silane Coatings on Calcium Silicate Hydrate Surfaces in Aggressive Environments: The Role of Nitrate/Sulfate Ratio

Silane is known as an effective coating for enhancing the resistance of concrete to harmful acids and radicals that are usually produced by the metabolism of microorganisms. However, the mechanism of silane protection is still unclear due to its nanoscale attributes. Here, the protective behavior of silane on the calcium silicate hydrate (C-S-H) surface is examined under the attack environment of nitrate/sulfate ions using molecular dynamics simulations. The findings revealed that silane coating improved the resistance of C-S-H to nitrate/sulfate ions. This resistance is considered the origin of silane protection against harmful ion attacks. Further research on the details of molecular structures suggests that the interaction between the oxygen in the silane molecule and the calcium in C-S-H, which can prevent the coordination of sulfate and nitrate to calcium on the C-S-H surface, is the cause of the silane molecules' strong adsorption. These results are also proved in terms of free energy, which found that the adsorption free energy on the C-S-H surface followed the order silane > sulfate > nitrate. This research confirms the excellent protection performance of silane on the nanoscale. The revealed mechanism can be further used to help the development of high-performance composite coatings.

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.

Predicting surfactant phase behavior with a molecularly informed field theory

HYPOTHESIS

The computational study of surfactants and self-assembly is challenging because 1) models need to reflect chemistry-specific interactions, and 2) self-assembled structures are difficult to equilibrate with conventional molecular dynamics. We propose to overcome these challenges with a multiscale simulation approach where relative entropy minimization transfers chemically-detailed information from all-atom (AA) simulations to coarse-grained (CG) models that can be simulated using field-theoretic methods. Field-theoretic simulations are not limited by intrinsic physical time scales like diffusion and allow for rigorous equilibration via free energy minimization. This approach should enable the study of properties that are difficult to obtain by particle-based simulations.

SIMULATION WORK

We apply this workflow to sodium dodecylsulfate. To ensure chemical fidelity we present an AA force field calibrated against interfacial tension experiments. We generate CG models from AA simulation trajectories and show that particle-based and field-theoretic simulations of the CG model reproduce AA simulations and experimental measurements.

FINDINGS

The workflow captures the complex balance of interactions in a multicomponent system ultimately described by an atomistic model. The resulting CG models can study complex 3D phases like double or alternating gyroids, and reproduce salt effects on properties like aggregation number and shape transitions.

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.

Free Energy Trajectory for Escape of a Single Chain from a Diblock Copolymer Micelle

We use umbrella sampling to compute the free energy trajectory of a single chain undergoing expulsion from an isolated diblock copolymer micelle. This approach elucidates the experimentally unobservable transition state, identifies the spatial position of the maximum free energy, and reveals the chain conformation of a single chain as it undergoes expulsion. Combining umbrella sampling with dissipative particle dynamics simulations of A₄B₈ micelles reveals that the core block (A) of the expelled chain remains partially stretched at the transition state, in contrast with the collapsed state assumed in some previous models. The free energy barrier increases linearly with the Flory-Huggins interaction parameter χ up to large interaction energies, where the structure of the otherwise spherical core apparently deforms near the transition state.