A. Gama Goicochea

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

The influence of systematic perturbation of input interaction parameters on thermodynamic equilibrium properties is studied employing dissipative particle dynamics (DPD) simulations. The values of both the excess pressure and the surface tension are found to be very sensitive to the values of the soft repulsion parameter between unlike DPD particles for high values of the coarse-graining level (number of water molecules per DPD particle). For the case in which a molecular surfactant is present at the interface we have determined the dependence of these properties on the values of the parameters that characterize the bonding force between polymer beads. No significant differences were found between linear and branched surfactants.

Modeling the temperature dependent interfacial tension between organic solvents and water using dissipative particle dynamics.

The interfacial tension between organic solvents and water at different temperatures is predicted using coarse-grained, mesoscopic Dissipative Particle Dynamics (DPD) simulations. The temperature effect of the DPD repulsive interaction parameters, aij, for the different components is calculated from the dependence of the Flory-Huggins χ parameter on temperature, by means of the solubility parameters. Atomistic simulations were carried out for the calculation of the solubility parameters for different organic compounds at different temperatures in order to estimate χ and then the aij coefficients. We validate this parametrization through the study of the interfacial tension in a mixture of benzene and water, and cyclohexane and water, varying the temperature. The predictions of our simulations are found to be in good agreement with experimental data taken from the literature, and show that the use of the solubility parameter at different temperatures to obtain the repulsive DPD parameters is a good alternative to introduce the effect of temperature in these systems.

Dissipative particle dynamics simulations of weak polyelectrolyte adsorption on charged and neutral surfaces as a function of the degree of ionization

The influence of the chain degree of ionization on the adsorption of weak polyelectrolytes on neutral and on oppositely and likely charged surfaces is investigated, by means of Monte Carlo simulations with a mesoscopic interaction model known as dissipative particle dynamics. The electrostatic interactions are calculated using the three-dimensional Ewald sum method, with an appropriate modification for confined systems. Effective wall forces confine the linear polyelectrolytes, and electric charges on the surfaces are included. The solvent, which is included explicitly, has a fluctuating particle number. The results show that the polyelectrolytes adsorb onto both neutral and charged surfaces, with the adsorption regulated by the chain degree of ionization, being larger at lower ionization degrees, where polyelectrolytes are less charged. Furthermore, polyelectrolyte adsorption is strongly modulated by the counterions screening of surface charges. These findings are supported by predictions of adsorption isotherms with varying ionization degrees. We obtain also the surface force mediated by adsorbed polyelectrolytes, which is calculated for the first time as a function of ionization degree. The adsorption and surface force isotherms obtained for weak polyelectrolytes are found to reproduce main trends in experiments, whenever those results are available, and provide additional insight into the role played by the competitive adsorption of the counterions and polyelectrolytes on the surfaces.

Coarse-grained simulations of the salt dependence of the radius of gyration of polyelectrolytes as models for biomolecules in aqueous solution.

The salt dependent radius of gyration of a polyelectrolyte in aqueous solution is calculated in an environment where the polyelectrolyte is surrounded by a permeable membrane that exchanges only solvent particles with the bulk. We obtain additionally the scaling exponent of the gyration radius as a function of the polymerization degree, and find that the polyelectrolyte retains a stretched conformation during the condensation and re-expansion process, indicating that these effects are of an electrostatic nature. The solvent quality is also shown to affect the polyelectrolyte conformation, especially for the poor solvent case. These results are obtained using a hybridized Monte Carlo technique with the coarse-grained, dissipative particle dynamics method with fluctuating number of solvent particles. The full range of the electrostatic interactions is included in the simulations, using the Ewald sum method, and the counterions and solvent molecules are included explicitly. In the complex systems mentioned above, the electrostatic interactions and the solvent quality play a key role in understanding phenomena that do not occur in uncharged systems. Our results are compared and validated with the behavior of some biomolecules under similar environments.

Colloidal stability dependence on polymer adsorption through disjoining pressure isotherms.

The disjoining pressure of polymers confined by colloidal walls was computed using dissipative particle dynamics simulations at constant chemical potential, volume, and temperature. The polymers are able to adsorb on the surfaces according to two models. In the so-called surface-modifying polymers, all monomers composing the chains have the same affinity for the substrate, whereas for the end-grafted polymer only the monomer at one of the ends of the polymer molecule adsorbs on the colloidal surface, resembling the behavior of dispersing agents. We find that these adsorption models yield markedly different disjoining pressure isotherms, which in turn predict different stability conditions for the colloidal dispersion. Our results show that for end-grafted polymers, a larger degree of polymerization at the same monomer concentration leads to better stability than for the surface-modifying ones. But also the unbound monomers of the surface-modifying type dominate over both kinds of polymers at large surface distances. The origin of these differences when the chemical nature of monomers is the same, and molecular weight and polymer concentration are used to characterize colloidal stability, is found to be mainly entropic.

Importance of Molecular Interactions in Colloidal Dispersions

We review briefly the concept of colloidal dispersions, their general properties and some of their most important applications, as well as the basic molecular interactions that give rise to their properties in equilibrium. Similarly, we revisit Brownian motion and hydrodynamic interactions associated with the concept of viscosity of colloidal dispersion. It is argued that the use of modern research tools, such as computer simulations, allows one to predict accurately some macroscopically measurable properties by solving relatively simple models of molecular interactions for a large number of particles. Lastly, as a case study, we report the prediction of rheological properties of polymer brushes using state of the art, coarse grained computer simulations, which are in excellent agreement with experiments.

On the computational modeling of the viscosity of colloidal dispersions and its relation with basic molecular interactions

The connection between fundamental interactions acting in molecules in a fluid and macroscopically measured properties, such as the viscosity between colloidal particles coated with polymers, is studied here. The role that hydrodynamic and Brownian forces play in colloidal dispersions is also discussed. It is argued that many body systems in which all these interactions take place can be accurately solved using computational simulation tools. One of those modern tools is the technique known as dissipative particle dynamics, which incorporates Brownian and hydrodynamic forces, as well as basic conservative interactions. A case study is reported, as an example of the applications of this technique, which consists of the prediction of the viscosity and friction between two opposing parallel surfaces covered with polymer chains, under the influence of a steady flow. This work is intended to serve as an introduction to the subject of colloidal dispersions and computer simulations, for last year undergraduate students and beginning graduate students who are interested in beginning research in soft matter systems. To that end, a computational code is included that students can use right away to study complex fluids in equilibrium.

Adsorption of polyelectrolytes on silica and gold surfaces

ABSTRACT The results of a study that helps understand the mechanisms of adsorption of polyelectrolytes on particles, using numerical simulation methods, specifically the one known as dissipative particle dynamics are reported here. The adsorption of cationic polyelectrolytes of two different polymerisation degrees interacting with two types of surfaces, one made of gold and the other of silica, is predicted and compared. We find that a more negatively charged wall does not necessarily adsorb more cationic polyelectrolytes because the electrostatic repulsion between the wall and the polyelectrolytes is stronger. Additionally, intra-chain repulsion plays an important role, because the largest polyelectrolyte chains have larger excluded volume than the shorter ones. In regard to the adsorption dependence on the polyelectrolyte polymerisation degree, we find that the excluded volume drives the adsorption throughout the intra-chain electrostatic repulsion, because the SiO2 surface is strongly negative. These results are expected to be useful for several nanotechnological applications of current interest, such as in gene therapy and in the improvement of drug delivering mechanisms.