Dissipative particle dynamics (DPD) study of the interfacial tension for alkane/water systems by using COSMO-RS to calculate interaction parameters

Dissipative particle dynamics parametrisation using infinite dilution activity coefficients: the impact of bonding

Dissipative particle dynamics (DPD) simulations have proven to be a valuable coarse-grained simulation technique for studying complex systems such as surfactant and polymer solutions. However, the best method to use in parametrising DPD systems is not universally agreed. One common approach is to map infinite dilution activity coefficients to the DPD simulation 'beads' that represent molecular fragments. However, we show that here that this approach can lead to serious errors when bonding beads together to create molecules. We show errors arise from the verlaps between bonded beads, which alters their solubility. In this article, we demonstrate how these bonding errors can be accounted for when defining DPD force fields using simple theoretical methods to account for the overlapping volumes, and we demonstrate the validity of our approach by calculating the partition coefficients for a series of solutes into two immiscible solvents.

Dissipative Particle Dynamics Simulation of Protein Folding in Explicit and Implicit Solvents: Coarse-Grained Model for Atomic Resolution

Advancements have been made to dissipative particle dynamics (DPD), a robust coarse-grained (CG) simulation method, to study the folded structures of four miniproteins (1L2Y, 1WN8, 1YRF, and 2I9M) in explicit and implicit solvents. In this endeavor, we aim to establish model parametrization and enhance computational efficiency. Unlike traditional CG models that use empirical force parameters, ex-force parameters (r0(ex), a ~ , δd, δp) of DPD particles constructed for specific research purposes can be obtained from atomistic molecular dynamics simulations. On the other hand, im-force parameters (r0(im), c, σ) can be derived from ex-DPD simulations, according to the underlying thermodynamic theory. Based on a mapping scheme proposed for the modeling of amino acids, all-atom proteins can be converted into a CG model. Both ex-/im-DPDs are then carried out to investigate the folding pathways of the four mini-proteins. Structural analysis of the RMSDs shows that the im-simulated proteins have greater structural similarity to native proteins than the ex-simulated ones. The constructed CG models achieve a resolution of Angstrom (Å), a level normally associated with atomic models. Additionally, speed tests reveal that im-DPD accelerates the simulation process and significantly improves simulation efficiency.

Many-body dissipative particle dynamics simulations of micellization of sodium alkyl sulfates

We present a study of micelle formation in alkyl sulfate surfactants using the simulation method of many-body dissipative particle dynamics (MDPD). We parametrise our model by tuning the intermolecular interactions in order to reproduce experimental values for the chemical potential and density at room temperature. Using this approach, we find that our model shows good agreement with experimental values for the critical micelle concentration (CMC). Furthermore, we show that our model can accurately predict CMC trends, which result from varying properties such as surfactant tail length and the salt concentration. We apply our model to investigate the effect of aggregation number on various micellar properties, such as the shape of individual micelles and the fraction of bound counterions. We show that micelles become aspherical at large aggregation numbers, in line with experimental predictions, and that longer tail surfactants are generally more spherical at all aggregation numbers compared to those which are shorter. We find excellent agreement between our simulations and experimental values for the degree of counterion binding, a factor that is crucial to accurately studying micellar shape, but one that is typically overlooked in the existing literature.

Instability of membranes containing ionizable cationic lipids: Effects of the repulsive range of headgroups and tail structures

The stability of membranes formed by ionizable cationic lipids, which constitute the primary components in lipid nanoparticles capable of endosomal escape, is explored using coarse-grained dissipative particle dynamics. Three types of ionizable model lipids with different tail structures are considered. Endosome acidification causes the ionization of lipids, leading to an increased repulsive range between their headgroups. When electrostatic repulsion is modeled as a conservative force with a long-range cutoff distance (rc,HH), the membrane and vesicle experience a loss of structural integrity and develop holes as rc,HH is beyond a critical value, which varies with the tail structure. When Coulombic repulsion is explicitly incorporated and intensified, a fully ionized lipid membrane undergoes a loss of structural integrity, displaying a qualitative similarity to the effect observed with the increase in rc,HH on the membrane stability. Qualitatively similar results are obtained for partially ionized membranes as the fraction of charged lipids increases. The stability of a mixed lipid membrane containing both ionizable and conventional lipids is also investigated. The disruption of the bilayer structure occurs for a sufficiently high charged fraction. The membrane instability can be attributed to the decrease in the packing parameter, which significantly deviates from unity as the interaction range increases.

Dissipative Particle Dynamics-Based Simulation of the Effect of Asphaltene Structure on Oil-Water Interface Properties

Asphaltenes are the main substances that stabilize emulsions during the production, processing, and transport of crude oil. The purpose of this research is to investigate the process of asphaltenes forming interfacial films at the oil-water interface by means of dissipative particle dynamics (DPD) and the effect of asphaltenes of different structures on the oil-water interface during the formation of interfacial film. It is demonstrated that the thickness of the interfacial film formed at the oil-water interface gradually increases as the asphaltene concentration rises and the amount of asphaltene adsorbed at the oil-water interface gradually multiplies. Both the number and type of heteroatoms in asphaltenes affect the interfacial behavior of asphaltenes. The interface activity of asphaltenes can be enhanced by increasing the number of heteroatoms in the asphaltene, and the type of heteroatom affects as well the interfacial activity of the asphaltene as it affects the aggregation behavior of the asphaltene in the system. As the number of asphaltene aromatic rings increases, the oil-water interfacial tension (IFT) trends down gradually, while the effect of alkyl side chains on the reduction of IFT of asphaltenes is different, and asphaltenes with medium length alkyl side chains can reduce IFT more efficiently.

A many-body dissipative particle dynamics parametrisation scheme to study behaviour at air-water interfaces

In this article, we present a general parametrisation scheme for many-body dissipative particle dynamics (MDPD). The scheme is based on matching model components to experimental surface tensions and chemical potentials. This allows us to obtain the correct surface and mixing behaviours of complex, multicomponent systems. The methodology is tested by modelling the behaviour of nonionic polyoxyethylene alkyl ether surfactants at an air/water interface. In particular, the influence of the number of ethylene oxide units in the surfactant head group is investigated. We find good agreement with many experimentally obtained parameters, such as minimum surface area per molecule; and a decrease in the surface tension with increasing surfactant surface density. Moreover, we observe an orientational transition, from surfactants lying directly on the water surface at low surface coverage, to surfactants lying parallel or tilted with respect to the surface normal at high surface coverage. The parametrisation scheme is also extended to cover the zwitterionic surfactant lauryldimethylamine oxide (LDAO), where we provide good predictions for the surface tension at maximum surface coverage. Here, if we exceed this coverage, we are able to demonstrate the spontaneous production of micelles from the surface surfactant layer.

Comprehensive review of the interfacial behavior of water/oil/surfactant systems using dissipative particle dynamics simulation

A comprehensive understanding of interfacial behavior in water/oil/surfactant systems is critical to evaluating the performance of emulsions in various industries, specifically in the oil and gas industry. To gain fundamental knowledge regarding this interfacial behavior, atomistic methods, e.g., molecular dynamics (MD) simulation, can be employed; however, MD simulation cannot handle phenomena that require more than a million atoms. The coarse-grained mesoscale methods were introduced to resolve this issue. One of the most effective mesoscale coarse-grained approaches for simulating colloidal systems is dissipative particle dynamics (DPD), which bridges the gap between macroscopic time and length scales and molecular-scale simulation. This work reviews the fundamentals of DPD simulation and its progress on colloids and interface systems, especially surfactant/water/oil mixtures. The effects of temperature, salt content, a water/oil ratio, a shear rate, and a type of surfactant on the interfacial behavior in water/oil/surfactant systems using DPD simulation are evaluated. In addition, the obtained results are also investigated through the lens of the chemistry of surfactants and emulsions. The outcome of this comprehensive review demonstrates the importance of DPD simulation in various processes with a focus on the colloidal and interfacial behavior of surfactants at water-oil interfaces.

Phase Diagram Study of Catanionic Surfactants Using Dissipative Particle Dynamics

Dissipative particle dynamics (DPD) simulations has been performed to study the phase transition of a mixture of cationic and anionic surfactants in an aqueous solution as a function of the total concentration in water and the relative ratio of surfactants. The impact of the relative difference between the tail lengths of the cationic and anionic surfactants on the phase diagram has been simulated by tuning the number of DPD beads in the simulation model. This research also discusses the impact of the frequently used values of the parameters associated with the harmonic bonds among the bonded DPD beads on the obtained self-assemblies. We find remarkable differences in the resultant self-assemblies based on different choices of harmonic bond parameters. The performed simulations show an enhanced spectrum of self-assemblies with augmented tail lengths and disparate harmonic bond parameters. The obtained self-assemblies are quite unique and can potentially be used in the future for various applications. We also compare the simulation results of the vesicle structures obtained by modeling the electrostatic interaction in the simulation among the charged beads by explicitly introducing charges with a long-range interaction with those obtained by tuning the implicit electrostatic interaction without the long-range interaction. The effects of the chain length of the model and the harmonic bond parameters on the internal density of DPD beads and stress profiles within the vesicles are examined closely. These results are a significant contribution to understanding the stability of the phases and tailoring of the desired vesicles.

Interfacial Assembly of Anisotropic Core-Shell and Hollow Microgels

Microgels, cross-linked polymers with submicrometer size, are ideal soft model systems. While spherical microgels have been studied extensively, anisotropic microgels have hardly been investigated. In this study, we compare the interfacial deformation and assembly of anisotropic core-shell and hollow microgels. The core-shell microgel consists of an elliptical core of hematite covered with a thin silica layer and a thin shell made of poly(N-isopropylacrylamide). The hollow microgels were obtained after a two-step etching procedure of the inorganic core. The behavior of these microgels at the oil-water interface was investigated in a Langmuir-Blodgett trough combined with ex situ atomic force microscopy. First, the influence of the architecture of anisotropic microgels on their spreading at the interface was investigated experimentally and by dissipative particle dynamic simulations. Hereby, the importance of the local shell thickness on the lateral and longitudinal interfacial deformation was highlighted as well as the differences between the core-shell and hollow architectures. The shape of the compression isotherms as well as the dimensions, ordering, and orientation of the microgels at the different compressions were analyzed. Due to their anisotropic shape and stiffness, both anisotropic microgels were found to exhibit significant capillary interactions with a preferential side-to-side assembly leading to stable microgel clusters at low interfacial coverage. Such capillary interactions were found to decrease in the case of the more deformable hollow anisotropic microgels. Consequently, anisotropic hollow microgels were found to distribute more evenly at high surface pressure compared to stiffer core-shell microgels. Our findings emphasize the complex interplay between the colloid design, anisotropy, and softness on the interfacial assembly and the opportunities it therefore offers to create more complex ordered interfaces.

Dissipative particle dynamics simulations in colloid and Interface science: a review

Dissipative particle dynamics (DPD) is one of the most efficient mesoscale coarse-grained methodologies for modeling soft matter systems. Here, we comprehensively review the progress in theoretical formulations, parametrization strategies, and applications of DPD over the last two decades. DPD bridges the gap between the microscopic atomistic and macroscopic continuum length and time scales. Numerous efforts have been performed to improve the computational efficiency and to develop advanced versions and modifications of the original DPD framework. The progress in the parametrization techniques that can reproduce the engineering properties of experimental systems attracted a lot of interest from the industrial community longing to use DPD to characterize, help design and optimize the practical products. While there are still areas for improvements, DPD has been efficiently applied to numerous colloidal and interfacial phenomena involving phase separations, self-assembly, and transport in polymeric, surfactant, nanoparticle, and biomolecules systems.

Micellization of Rhamnolipid Biosurfactants and Their Applications in Oil Recovery: Insights from Mesoscale Simulations

The dissipative particle dynamics (DPD) mesoscopic method is used to investigate the self-assembly of rhamnolipid congeners and their aggregation behaviors with paraffins including nonane and pentadecane. The coarse-grained force field is parameterized by combining molecular dynamics (MD) simulations, COSMOtherm calculations, and available experimental data. This model reproduces the vesicular formation of α-l-rhamnopyranosyl-β-hydroxydecanoyl-β-hydroxydecanoate (Rha-C10-C10) reported by all-atom MD simulations. The vesicle composed of Rha-C10-C10 is found to be most stable at a surfactant concentration of 100-146 mM based on asphericity analysis. The architecture of rhamnolipid congeners affects the morphology of their aggregates. Di-rhamno-di-lipidic dRha-C16-C16 forms vesicles with a thicker unilamellar layer of 3.2 nm. Rha-C16-C16 forms vesicles at a lower concentration of 70 mM, but the enclosed water space collapses when the surfactant concentration increases. dRha-C10-C10 forms wormlike micelles, which agglomerate into a torus and interconnected network at higher concentrations. In the presence of alkane molecules, dRha-C10-C10 maintains its wormlike micellar morphology with alkane molecules wrapped inside the aggregates. For Rha-C10-C10, Rha-C16-C16, and dRha-C16-C16, nonane molecules are distributed in the hydrophobic subdomain formed by rhamnolipid molecules. Spherical vesicles are formed at a surfactant concentration of 50 mM and then develop into ellipsoidal vesicles when the concentration increases to 125 mM. When mixed with pentadecane, the alkane molecules are aggregated and surrounded by surfactants forming a core-shell structure at a low surfactant concentration of 20 mM. At higher alkane and surfactant concentrations, the morphologies develop into disk micelles, wormlike micelles, and vesicles, with pentadecane molecules being distributed and packed with rhamnolipids. The obtained simulation results suggest that these biosurfactants have potential as environmental remediation agents.

A Surface Site Interaction Point Method for Dissipative Particle Dynamics Parametrization: Application to Alkyl Ethoxylate Surfactant Self-Assembly

Dissipative particle dynamics (DPD) is a coarse-grained approach to the simulation of large supramolecular systems, but one limitation has been that the parameters required to describe the noncovalent interactions between beads are not readily accessible. A first-principles computational method has been developed so that bead interaction parameters can be calculated directly from ab initio gas-phase molecular electrostatic potential surfaces of the molecular fragments that represent the beads. A footprinting algorithm converts the molecular electrostatic potential surfaces into a discrete set of surface site interaction points (SSIPs), and these SSIPs are used in the SSIMPLE (surface site interaction model for the properties of liquids at equilibrium) algorithm to calculate the free energies of transfer of one bead into a solution of any other bead. The bead transfer free energies are then converted into the required DPD interaction parameters for all pairwise combinations of different beads. The reliability of the parameters was demonstrated using DPD simulations of a range of alkyl ethoxylate surfactants. The simulations reproduce the experimentally determined values of the critical micelle concentration and mean aggregation number well for all 22 surfactants studied.

COSMOplex: self-consistent simulation of self-organizing inhomogeneous systems based on COSMO-RS

During the past 20 years, the efficient combination of quantum chemical calculations with statistical thermodynamics by the COSMO-RS method has become an important alternative to force-field based simulations for the accurate prediction of free energies of molecules in liquid systems. While it was originally restricted to homogeneous liquids, it later has been extended to the prediction of the free energy of molecules in inhomogeneous systems such as micelles, biomembranes, or liquid interfaces, but these calculations were based on external input about the structure of the inhomogeneous system. Here we report the rigorous extension of COSMO-RS to a self-consistent prediction of the structure and the free energies of molecules in self-organizing inhomogeneous systems. This extends the application range to many new areas, such as the prediction of micellar structures and critical micelle concentrations, finite loading effects in micelles and biomembranes, the free energies and structure of liquid interfaces, microemulsions, and many more related topics, which often are of great practical importance.

Simulations of Interfacial Tension of Liquid-Liquid Ternary Mixtures Using Optimized Parametrization for Coarse-Grained Models

In this work, liquid-liquid systems are studied by means of coarse-grained Monte Carlo simulations (CG-MC) and Dissipative Particle Dynamics (DPD). A methodology is proposed to reproduce liquid-liquid equilibrium (LLE) and to provide variation of interfacial tension (IFT), as a function of the solute concentration. A key step is the parametrization method based on the use of the Flory-Huggins parameter between DPD beads to calculate solute/solvent interactions. Parameters are determined using a set of experimental compositional data of LLE, following four different approaches. These approaches are evaluated, and the results obtained are compared to analyze advantages/disadvantages of each one. These methodologies have been compared through their application on six systems: water/benzene/1,4-dioxane,water/chloroform/acetone, water/benzene/acetic acid, water/benzene/2-propanol, water/hexane/acetone, and water/hexane/2-propanol. CG-MC simulations in the Gibbs (NVT) ensemble have been used to check the validity of parametrization approaches for LLE reproduction. Then, CG-MC simulations in the osmotic (μ_soluteN_solventP _zzT) ensemble were carried out considering the two liquid phases with an explicit interface. This step allows one to work at the same bulk concentrations as the experimental data by imposing the precise bulk phase compositions and predicting the interface composition. Finally, DPD simulations were used to predict IFT values for different solute concentrations. Our results on variation of IFT with solute concentration in bulk phases are in good agreement with experimental data, but some deviations can be observed for systems containing hexane molecules.