DPD Parameters Estimation for Simultaneously Simulating Water-Oil Interfaces and Aqueous Nonionic Surfactants

The outcome of a coarse-grained simulation within the dissipative particle dynamics framework strongly depends on the choice of the repulsive parameter between different species. Different methodologies have been used in the literature to determine these parameters toward reproducing selected experimental system properties. In this work, a systematic investigation on possible procedures for estimating the simulation parameters is conducted. We compare methods based on the Hildebrand and the Hansen solubility parameter theories, mapped into the Flory-Huggins model. We find that using the Hansen solubility parameters it is possible to achieve a high degree of coarse graining, with parameters that yield realistic values for the interfacial tension. The procedure was first applied to the water/benzene system and then validated for water/n-octane, water/1,1-dichloroethane, water/methyl cyclohexane, and water/isobutyl acetate. In all these cases, the experimental interfacial tension could be reproduced by adjusting a single correction factor. In the case of the water benzene system, the dissipative particle dynamics parameters derived using our approach were able to simultaneously describe both the interfacial tension and micellar properties of aqueous nonionic surfactants representative of the octyl poly(ethylene oxide) C8H17O(C2H4O)(m)H family. We show how the parameters can be used, within the dissipative particle dynamics framework, to simulate the water/oil interface in the presence of surfactants at varying concentrations. The results show, as expected, that as the surfactant concentration increases, the interfacial tension decreases, and micelles form in bulk water.