Salt- and pH-Dependent Viscosity of SDS/LAPB Solutions: Experiments and a Semiempirical Thermodynamic Model.

We present novel data on the composition-, pH-, and salt-dependent zero shear viscosity of the commercially important mixture of anionic sodium dodecyl sulfate (SDS) and zwitterionic lauramidopropyl betaine (LAPB). We show via proton NMR experiments that the notionally zwitterionic LAPB exhibits a large pKa shift in the presence of SDS and can become partially cationic at formulation-relevant pH ranges of 4.5-6.0-that is, the binary system is effectively a ternary system. This has a pronounced effect on the viscosity of the system at low pH, especially if the fraction of LAPB is high. We use theoretical arguments to motivate a semiempirical but practical approach to model the viscosity of the mixtures using thermodynamic parameters such as the excess chemical potentials or activity coefficients of the surfactants. We demonstrate this using an augmented regular solution theory-based mixed micelle thermodynamic model and develop robust regression models using Bayesian approaches. We also show how the pKa shift from NMR experiments can be used to parameterize the thermodynamic model. This framework should be extensible to other arbitrary surfactant mixtures in the future and hence will be of broad interest for the development of surfactant formulations for household, personal care, and other applications.

A new equation of state for homo-polymers in dissipative particle dynamics.

A chain-revised Groot-Warren equation of state (crGW-EOS) was developed and tested to describe systems of homo-oligomeric chains in the framework of dissipative particle dynamics (DPD). First, thermodynamic perturbation theory is applied to introduce correction terms that account for the reduction in pressure with an increasing number of bonds at constant bead number density. Then, this EOS is modified by introducing a set of switching functions that yields an accurate second virial coefficient in the low-density limit. The crGW-EOS offers several improvements over the revised Groot-Warren equation of state and Groot-Warren equation of state for chain molecules. We tested the crGW-EOS by using it to predict the pressure of oligomeric systems and the B2 virial coefficient of chain DPD particles for a range of bond lengths. Additionally, a method is developed for determining the strength of cross-interaction parameters between chains of different compositions and sizes and for thermal and athermal mixtures. We explored how different levels of coarse-graining affect the upper-critical solution temperature.

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 Henry's 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.

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.

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 log POW < 2 partition primarily to the head group region and shrink micellar length, decreasing viscosity substantially, whereas more hydrophobic PRMs, with log POW 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 log POW > 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.

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 pKa's 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.

Atomistic simulation of solid-liquid coexistence for molecular systems: application to triazole and benzene.

A method recently developed to rigorously determine solid-liquid equilibrium using a free-energy-based analysis has been extended to analyze multiatom molecular systems. This method is based on using a pseudosupercritical transformation path to reversibly transform between solid and liquid phases. Integration along this path yields the free energy difference at a single state point, which can then be used to determine the free energy difference as a function of temperature and therefore locate the coexistence temperature at a fixed pressure. The primary extension reported here is the introduction of an external potential field capable of inducing center of mass order along with secondary orientational order for molecules. The method is used to calculate the melting point of 1-H-1,2,4-triazole and benzene. Despite the fact that the triazole model gives accurate bulk densities for the liquid and crystal phases, it is found to do a poor job of reproducing the experimental crystal structure and heat of fusion. Consequently, it yields a melting point that is 100 K lower than the experimental value. On the other hand, the benzene model has been parametrized extensively to match a wide range of properties and yields a melting point that is only 20 K lower than the experimental value. Previous work in which a simple "direct heating" method was used actually found that the melting point of the benzene model was 50 K higher than the experimental value. This demonstrates the importance of using proper free energy methods to compute phase behavior. It also shows that the melting point is a very sensitive measure of force field quality that should be considered in parametrization efforts. The method described here provides a relatively simple approach for computing melting points of molecular systems.

Toward a robust and general molecular simulation method for computing solid-liquid coexistence.

A rigorous and generally applicable method for computing solid-liquid coexistence is presented. The method overcomes some of the technical difficulties associated with other solid-liquid simulation procedures and can be implemented within either a molecular dynamics or Monte Carlo framework. The method consists of three steps: First, relative Gibbs free energy curves are created for the solid and liquid phases using histogram reweighting. Next, the free energy difference between the solid and liquid phases is evaluated at a single state point by integrating along a pseudosupercritical transformation path that connects the two phases. Using this result, the solid and liquid free energy curves are referenced to a common point, allowing a single coexistence point to be determined. Finally, Gibbs-Duhem integration is used to determine the full coexistence curve. To evaluate its utility, this method is applied to the Lennard-Jones and NaCl systems. Results for solid-liquid coexistence agree with previous calculations for these systems. In addition, it is shown that the NaCl model does not correctly describe solid-liquid coexistence at high pressures. An analysis of the accuracy of the method indicates that the results are most sensitive to the transformation free energy calculation.