Dependence of thermodynamic properties of model systems on some dissipative particle dynamics parameters

Compression and Ordering of Hollow Microgels in Monolayers Formed at Liquid-Liquid Interfaces

Monolayers of polymer microgels with a spherical cavity adsorbed at the liquid-liquid interface were studied using mesoscopic computer simulations. One liquid, named water, was always considered as a good solvent, while the microgel solubility in the second liquid, named oil, was varied. The symmetric and asymmetric cases of vanishing and the strong differences in solubility between the network particles and the liquids were considered. The simulations provided us with an insight into the shape and volume changes of the microgels upon compression, making it possible to relate the response of the individual network with the collective order and structure of the monolayer. Similar to regular microgels, the compression of the monolayer of hollow particles led to a decrease in lateral sizes accompanied by shape transformation from a flattened to a nearly spherical shape. However, the presence of a cavity filled with solvent caused some unique differences in the behavior of the system. The adsorption pathway of hollow microgels at the liquid interface predefines: (a) the position of the particles with respect to the interface and (b) the structure of the monolayer. A striking discovery is that in the symmetric case of similar solubility of the microgel in both liquids, it is possible to produce a monolayer in which one part of the network faces the aqueous phase and the other part faces the oil phase. The polymer concentration profiles plotted along the normal to the interface reveal a redistribution of polymeric mass of the microgels relative to the interface, distinguishing between the microgels whose cavities are filled with water and oil, respectively. Moreover, the ratio between the microgels faced in water and oil does not change upon compression and predetermines the response and order of the monolayer.

Investigating anionic surfactant phase diagrams using dissipative particle dynamics: development of a transferable model

Dissipative particle dynamics (DPD) provides a powerful coarse-grained simulation technique for the study of a wide range of soft matter systems. Here, we investigate the transferability of DPD models to the prediction of anionic surfactant phase diagrams, taking advantage of fast parameter sweeps to optimise the choice of DPD parameters for these systems. Parameters are developed which provide a good representation of the phase diagrams of SDS (sodium dodecyl sulfate) and three different isomeric forms of LAS (linear alkylbenzene sulfonates) across an extensive concentration range. A high degree of transferability is seen, with parameters readily transferable to other systems, such as AES (alkyl ether sulfates). Excellent agreement is obtained with experimentally measured quantities, such as the lamellar layer spacing. Isosurfaces are produced from the surfactant head group, from which the second moment M of the isosurface normal distribution is calculated for different phase structures. Lyotropic liquid crystalline phases are characterised by a combination of the eigenvalues of M, radial distribution functions, and visual inspections.

Swelling, collapse and ordering of rod-like microgels in solution: Computer simulation studies

Polymer microgels have proven to be highly promising macromolecular objects for a wide variety of applications. In particular, the soft particles of an anisotropic (rod-like) shape are of special interest because of their potential use in tissue engineering or materials design. However, a little is known about the physical behavior of such microgels in solution, which inspired us to study them using mesoscopic computer simulations. For single networks, depending on the solvent quality, the dimensional characteristics were obtained for microgels of different molecular weight, crosslinking density and aspect ratio. In particular, the conditions for the rod-to-rod (preserving the nonspherical shape) and rod-to-sphere collapse were found. In addition, the effect of the liquid-crystalline (LC) ordering was demonstrated for the ensemble of rod-like microgels at different swelling ratios, and the influence of microgel aspect ratio on the volume fraction of the LC transition was shown.

The viscosity of polyelectrolyte solutions and its dependence on their persistence length, concentration and solvent quality

In this work, a comprehensive study of the influence on shear viscosity of polyelectrolyte concentration, persistence length, salt concentration and solvent quality is reported, using numerical simulations of confined solutions under stationary Poiseuille flow. Various scaling regimes for the viscosity are reproduced, both under good solvent and theta solvent conditions. The key role played by the electrostatic interactions in the viscosity is borne out when the ionic strength is varied. It is argued that these results are helpful for the understanding of viscosity scaling in entangled polyelectrolyte solutions for both rigid and flexible polyelectrolytes in different solvents, which is needed to perform intelligent design of new polyelectrolytes capable of fine tuning the viscosity in complex fluids.

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.

Modeling of Branched Thickening Polymers under Poiseuille Flow Gives Clues as to How to Increase a Solvent's Viscosity

The viscosity enhancement of a solvent produced by the addition of thickening branched polymers is predicted as a function of polymer concentration, branch length and persistence length, and strength of the covalent bonding interactions. Nonequilibrium, stationary-state Poiseuille numerical simulations are performed using the dissipative particle dynamics model to obtain the viscosity of the fluid. It is found that the clustering of the polymers into aggregates increases the viscosity and that it is more strongly affected by the strength of the bonding interactions. General scaling relationships are found for the viscosity as a function of the variables studied, which are expected to be useful for the design and synthesis of new viscosifying polymers. It is argued that our results can be applied to aqueous thickeners, of importance for colloidal fluids such as paints and coatings and also for nonpolar fluids such as supercritical CO₂, which is a promising nonhydraulic fracking fluid also useful in enhanced oil recovery.

Compression and Ordering of Microgels in Monolayers Formed at Liquid-Liquid Interfaces: Computer Simulation Studies

Monolayers of polymer microgels adsorbed at the liquid interfaces were studied by dissipative particle dynamics simulations. The results demonstrated that the compressibility of the monolayers can be widely tuned by varying the cross-linking density of the microgels and their (in)compatibility with the immiscible liquids. In particular, the compression of the monolayers (increase of 2D concentration of the microgels) leads to the decrease of their lateral size. Herewith, the shape of the individual soft particles gradually changes from oblate (diluted 2D system) to nearly spherical (compressed monolayer). The polymer concentration profiles plotted along the normal to the interface reveal a nonmonotonous shape with a sharp maximum at the interface. This is a consequence of the shielding effect: saturation of the interface by monomer units of the subchains is driven by minimization of unfavorable contacts between the immiscible liquids and is opposed by elasticity of the network. The decrease of the interfacial tension upon concentration (compression) of the monolayer is quantified. It has been demonstrated that the interfacial tension significantly differs if the solubility of the polymer chains of the microgel network in the liquids changes. These results correlate well with experimental data. The examination of the microgels' crystalline ordering in monolayers demonstrated a nonmonotonous dependency on the compression degree (microgel concentration). Finally, the worsening of the solvent quality leads to the collapse of the microgels in monolayer and nonmonotonous behavior of the interfacial tension.

Quantification of Ostwald Ripening in Emulsions via Coarse-Grained Simulations

Ostwald ripening is a diffusional mass transfer process that occurs in polydisperse emulsions, often with the result of threatening the emulsion stability. In this work, we design a simulation protocol that is capable of quantifying the process of Ostwald ripening at the molecular level. To achieve experimentally relevant time scales, the dissipative particle dynamics (DPD) simulation protocol is implemented. The simulation parameters are tuned to represent two benzene droplets dispersed in water. The coalescence between the two droplets is prevented via the introduction of membranes, which allow diffusion of benzene from one droplet to the other. The simulation results are quantified in terms of the changes in the droplet volume as a function of time. The results are in qualitative agreement with experiments. The agreement with the Lifshitz-Slyozov-Wagner theory becomes quantitative when the simulated solubility and diffusion coefficient of benzene-in-water are considered. The effect of two different surfactants was also investigated. In agreement with both experimental observations and theory, the addition of surfactants at moderate concentrations decreased the Ostwald ripening rate because of the reduction in the interfacial tension between benzene and water; as the surfactant film becomes dense, other phenomena are likely to further delay the Ostwald ripening. In fact, the results suggest that the surfactant that yields higher density at the benzene-water interface delayed more effectively Ostwald ripening. The formation of micelles can also affect the ripening rate, in qualitative agreement with experiments, although our simulations are not conclusive on such effects. Our simulations show that the coarse-grained DPD formalism is able to capture the molecular phenomena related to Ostwald ripening and reveal molecular level features that could help to understand experimental observations. The results could be useful for predicting and eventually controlling the long-term stability of emulsions.

Assessing the transferability of common top-down and bottom-up coarse-grained molecular models for molecular mixtures

Assessing the performance of top-down and bottom-up coarse-graining approaches.

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.

Surfactant chain length and concentration influence on the interfacial tension of two immiscible model liquids: a coarse-grained approach

The interfacial tension between immiscible liquids is studied as a function of a model linear surfactant length and concentration using coarse-grained, dissipative particle dynamics numerical simulations. The adsorption isotherms obtained from the simulations are found to be in agreement with Langmuir's model. The reduction of the interfacial tension with increasing surfactant concentration is found to display some common characteristics for all the values of chain length modeled, with our predictions being in agreement with Szyszkowski's equation. Lastly, the critical micelle concentration is predicted for all surfactant lengths, finding exponentially decaying behavior, in agreement with Kleven's model. It is argued that these findings can be helpful guiding tools in the interpretation of available experiments and in the design of new ones with new surfactants and polymers.

Electrostatics in dissipative particle dynamics using Ewald sums with point charges

A proper treatment of electrostatic interactions is crucial for the accurate calculation of forces in computer simulations. Electrostatic interactions are typically modeled using Ewald-based methods, which have become some of the cornerstones upon which many other methods for the numerical computation of electrostatic interactions are based. However, their use with charge distributions rather than point charges requires the inclusion of ansatz for the solutions of the Poisson equation, since there is no exact solution known for smeared out charges. The interest in incorporating electrostatic interactions at the scales of length and time that are relevant for the study the physics of soft condensed matter has increased considerably. Using mesoscale simulation techniques, such as dissipative particle dynamics (DPD), allows us to reach longer time scales in numerical simulations, without abandoning the particulate description of the problem. The main problem with incorporating electrostatics into DPD simulations is that DPD particles are soft and those particles with opposite charge can form artificial clusters of ions. Here we show that one can incorporate the electrostatic interactions through Ewald sums with point charges in DPD if larger values of coarse-graining degree are used, where DPD is truly mesoscopic. Using point charges with larger excluded volume interactions, the artificial formation of ionic pairs with point charges can be avoided and one obtains correct predictions. We establish ranges of parameters useful for detecting boundaries where artificial formation of ionic pairs occurs. Lastly, using point charges we predict the scaling properties of polyelectrolytes in solvents of varying quality, and obtain predictions that are in agreement with calculations that use other methods and with recent experimental results.

Induced polar order in sedimentation equilibrium of rod-like nanoswimmers

The active diffusion and sedimentation equilibrium of rod-like nanoswimmers with length L are investigated by dissipative particle dynamics. In the absence of propulsion, the nanorod has the rotational correlation time τθ and mobility μ0, which vary with L. On the basis of the mean squared displacement, the diffusive behavior of rod-like nanoswimmers subject to active force Fa is found to follow (D - D0) ∝ (μ0Fa)(2)τθ, where D0 depicts the Brownian diffusivity. When the nanoswimmer suspension is under the external force Fe, the balance between the downward migration and upward active diffusion yields the sedimentation length δ = D/(μ0Fe), which no longer obeys δ0 = kBT/Fe obtained from the Einstein-Smoluchowski relationship, D0 = μ0kBT. Different from the suspension of passive rods, the polar order is clearly seen for active rods. The local polar order is essentially constant within the distance of about 2δ from the bottom wall but decays as the distance is further increased. In this work, the active Peclet number is small compared to unity and the maximum polar order grows linearly with Fe/Fa. The polar order arises because it is easier for the nanoswimmer with the swimming direction opposite to the external force to escape from the bottom wall.

Nanotribology of biopolymer brushes in aqueous solution using dissipative particle dynamics simulations: an application to PEG covered liposomes in a theta solvent

We undertake the investigation of sheared polymer chains grafted onto flat surfaces to model liposomes covered with polyethylene glycol brushes as a case study for the mechanisms of efficient drug delivery in biologically relevant situations, for example, as carriers for topical treatments of illnesses in the human vasculature. For these applications, specific rheological properties are required, such as low viscosity at high shear rates, to improve the transport of the liposomes. Therefore, extensive non-equilibrium, coarse-grained dissipative particle dynamics simulations of polymer brushes of various lengths and shear rates are performed to obtain the average viscosity and the friction coefficient of the system as functions of the shear rate and polymerization degree under theta-solvent conditions, and we find that the brushes experience considerable shear thinning at large shear rates. The viscosity (η) and the friction coefficient (μ) are shown to obey the scaling laws η ∼ γ dot above (-0.31) and μ ∼ γ dot above (0.69) at high shear rates (γ dot above) in a theta solvent, irrespective of the degree of polymerization of brushes. These results confirm recent scaling predictions and reproduce very well trends in measurements of the viscosity at a high shear rate (γ dot above) of red blood cells in a liposome containing medium.