No abstract

No abstract

We present a systematic, top-down, thermodynamic parametrization scheme for dissipative particle dynamics (DPD) using water-octanol partition coefficients, supplemented by water-octanol phase equilibria and pure liquid phase density data. We demonstrate the feasibility of computing the required partition coefficients in DPD using brute-force simulation, within an adaptive semi-automatic staged optimization scheme. We test the methodology by fitting to experimental partition coefficient data for twenty one small molecules in five classes comprising alcohols and poly-alcohols, amines, ethers and simple aromatics, and alkanes (i.e., hexane). Finally, we illustrate the transferability of a subset of the determined parameters by calculating the critical micelle concentrations and mean aggregation numbers of selected alkyl ethoxylate surfactants, in good agreement with reported experimental values.

No abstract

We use dissipative particle dynamics (DPD) to study micelle formation in alkyl sulfate surfactants, with alkyl chain lengths ranging from 6 to 12 carbon atoms. We extend our recent DPD force field [ J. Chem. Phys. 2017 , 147 , 094503 ] to include a charged sulfate chemical group and aqueous sodium ions. With this model, we achieve good agreement with the experimentally reported critical micelle concentrations (CMCs) and can match the trend in mean aggregation numbers versus alkyl chain length. We determine the CMC by fitting a charged pseudophase model to the dependence of the free surfactant on the total surfactant concentration above the CMC and compare it with a direct operational definition of the CMC as the point at which half of the surfactant is classed as micellar and half as monomers and submicellar aggregates. We find that the latter provides the best agreement with experimental results. Finally, with the same model, we are able to observe the sphere-to-rod morphological transition for sodium dodecyl sulfate (SDS) micelles and determine that it corresponds to SDS concentrations in the region of 300-500 mM.