David Michael Eike

Liquid–liquid equilibria for soft-repulsive particles: Improved equation of state and methodology for representing molecules of different sizes and chemistry in dissipative particle dynamics

Three developments are presented that significantly expand the applicability of dissipative particle dynamics (DPD) simulations for symmetric and non-symmetric mixtures, where the former contain particles with equal repulsive parameter for self-interactions but a different repulsive parameter for cross-interactions, and the latter contain particles with different repulsive parameters also for the self-interactions. Monte Carlo and molecular dynamics simulations for unary phases covering a wide range of repulsive parameters and of densities for single-bead DPD particles point to deficiencies of the Groot and Warren equation of state (GW-EOS) [J. Chem. Phys. 107, 4423 (1997)]. A revised version, called rGW-EOS, is proposed here that is significantly more accurate over a wider range of parameters/densities. The second development is the generalization of the relationship between the Flory-Huggins χ parameter and the repulsive cross-interaction parameter when the two particles involved have different molecular volumes. The third aspect is an investigation of Gibbs ensemble Monte Carlo simulation protocols, which demonstrates the importance of volume fluctuations and excess volumes of mixing even for equimolar symmetric mixtures of DPD particles. As an illustrative example, the novel DPD methodology is applied to the prediction of the liquid-liquid equilibria for acetic anhydride/(n-hexane or n-octane) binary mixtures.

Accurate modeling of ionic surfactants at high concentration.

Molecular dynamics (MD) simulation is a useful tool for simulating formulations of surfactant mixtures from first-principles, which can be used to predict surfactant morphology and other industrially relevant thermodynamic properties. However, the surfactant structure is sensitive to the parameters used in MD simulations, and in the absence of extensive validation against experimental data, it is often not obvious a priori which range of parameter sets to choose. In this work, we compare the performance of ion parameters implemented in nonpolarizable classical MD simulations, and its effect on simulations of an idealized solution of sodium dodecyl sulfate (SDS). We find that previous artifacts reported in simulations of larger SDS constructs are a direct consequence of using parameters that poorly model ionic interactions at high concentration. Using osmotic pressure and/or other thermodynamic properties measured at finite concentration, such as Kirkwood-Buff integrals, is shown to be the most cost-effective means to validate and parametrize existing force fields. Our findings highlight the importance of optimizing intermolecular parameters for simulations of systems with a high local concentration, which may be applicable in other contexts, such as in molecular crowding, hotspot mapping, protein folding, and modeling pH effects.

Predicting proton titration in cationic micelle and bilayer environments.

Knowledge of the protonation behavior of pH-sensitive molecules in micelles and bilayers has significant implications in consumer product development and biomedical applications. However, the calculation of pKas in such environments proves challenging using traditional structure-based calculations. Here we apply all-atom constant pH molecular dynamics with explicit ions and titratable water to calculate the pKa of a fatty acid molecule in a micelle of dodecyl trimethylammonium chloride and liquid as well as gel-phase bilayers of diethyl ester dimethylammonium chloride. Interestingly, the pKa of the fatty acid in the gel bilayer is 5.4, 0.4 units lower than that in the analogous liquid bilayer or micelle, despite the fact that the protonated carboxylic group is significantly more desolvated in the gel bilayer. This work illustrates the capability of all-atom constant pH molecular dynamics in capturing the delicate balance in the free energies of desolvation and Coulombic interactions. It also shows the importance of the explicit treatment of ions in sampling the protonation states. The ability to model dynamics of pH-responsive substrates in a bilayer environment is useful for improving fabric care products as well as our understanding of the side effects of anti-inflammatory drugs.

Multiscale Modeling of the Effects of Salt and Perfume Raw Materials on the Rheological Properties of Commercial Threadlike Micellar Solutions

We link micellar structures to their rheological properties for two surfactant body-wash formulations at various concentrations of salts and perfume raw materials (PRMs) using molecular simulations and micellar-scale modeling, as well as traditional surfactant packing arguments. The two body washes, namely, BW-1EO and BW-3EO, are composed of sodium lauryl ethylene glycol ether sulfate (SLEnS, where n is the average number of ethylene glycol repeat units), cocamidopropyl betaine (CAPB), ACCORD (which is a mixture of six PRMs), and NaCl salt. BW-3EO is an SLE3S-based body wash, whereas BW-1EO is an SLE1S-based body wash. Additional PRMs are also added into the body washes. The effects of temperature, salt, and added PRMs on micellar lengths, breakage times, end-cap free energies, and other properties are obtained from fits of the rheological data to predictions of the Pointer Algorithm [ Zou , W. ; Larson , R.G. J. Rheol. 2014 , 58 , 1 - 41 ], which is a simulation method based on the Cates model of micellar dynamics. Changes in these micellar properties are interpreted using the Israelachvili surfactant packing argument. From coarse-grained molecular simulations, we infer how salt modifies the micellar properties by changing the packing between the surfactant head groups, with the micellar radius remaining nearly constant. PRMs do so by partitioning to different locations within the micelles according to their octanol/water partition coefficient POW and chemical structures, adjusting the packing of the head and/or tail groups, and by changing the micelle radius, in the case of a large hydrophobic PRM. We find that relatively hydrophilic PRMs with logu2009POW < 2 partition primarily to the head group region and shrink micellar length, decreasing viscosity substantially, whereas more hydrophobic PRMs, with logu2009POW between 2 and 4, mix with the hydrophobic surfactant tails within the micellar core and slightly enhance the viscosity and micelle length, which is consistent with the packing argument. Large and very hydrophobic PRMs, with logu2009POW > 4, are isolated deep inside the micelle, separating from the tails and swelling the radius of the micelle, leading to shorter micelles and much lower viscosities, leading eventually to swollen-droplet micelles.

Scission Free Energies for Wormlike Surfactant Micelles: Development of a Simulation Protocol, Application, and Validation for Personal Care Formulations

We present a scheme to calculate wormlike micelle scission free energies from a potential of mean force (PMF) derived from a weighted histogram analysis method (WHAM) applied to coarse grained dissipative particle dynamics (DPD) simulations. In contrast to previous related work, we use a specially chosen external potential based on a reaction coordinate that reversibly drives surfactants out of the nascent scission location. For the application to a model body wash formulation, we predict how addition of NaCl and small molecules such as perfume raw materials (PRMs) affect scission energies. The results show qualitative agreement and correct trends compared to recently determined scission energies for the same system; however, a more rigorous parametrization of the underlying DPD potential is required for quantitative agreement.

Probing Additive Loading in the Lamellar Phase of a Nonionic Surfactant: Gibbs Ensemble Monte Carlo Simulations Using the SDK Force Field

Understanding solute uptake into soft microstructured materials, such as bilayers and worm-like and spherical micelles, is of interest in the pharmaceutical, agricultural, and personal care industries. To obtain molecular-level insight on the effects of solutes loading into a lamellar phase, we utilize the Shinoda-Devane-Klein (SDK) coarse-grained force field in conjunction with configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The lamellar phase is comprised of a bilayer formed by triethylene glycol mono- n-decyl ether (C10E3) surfactants surrounded by water with a 50:50 surfactant/water weight ratio. We study both the unary adsorption isotherm and the effects on bilayer structure and stability caused by n-nonane, 1-hexanol, and ethyl butyrate at several different reduced reservoir pressures. The nonpolar n-nonane molecules load near the center of the bilayer. In contrast, the polar 1-hexanol and ethyl butyrate molecules both load with their polar bead close to the surfactant head groups. Near the center of the bilayer, none of the solute molecules exhibits a significant orientational preference. Solute molecules adsorbed near the polar groups of the surfactant chains show a preference for orientations perpendicular to the interface, and this alignment with the long axis of the surfactant molecules is most pronounced for 1-hexanol. Loading of n-nonane leads to an increase of the bilayer thickness, but does not affect the surface area per surfactant. Loading of polar additives leads to both lateral and transverse swelling. The reduced Henrys law constants of adsorption (expressed as a molar ratio of additive to surfactant per reduced pressure) are 0.23, 1.4, and 14 for n-nonane, 1-hexanol, and ethyl butyrate, respectively, and it appears that the SDK force field significantly overestimates the ethyl butyrate-surfactant interactions.