Martin Turesson

Calcium Mediated Polyelectrolyte Adsorption on Like-Charged Surfaces

Monte Carlo simulations within the primitive model of calcium-mediated adsorption of linear and comb polyelectrolytes onto like-charged surfaces are described, focusing on the effect of calcium and polyion concentrations as well as on the ion pairing between polymers and calcium ions. We use a combination of Monte Carlo simulations and experimental data from titration and calcium binding to quantify the ion pairing. The polymer adsorption is shown to occur as a result of surface overcharging by Ca(2+) and ion pairing between charged monomers and Ca(2+). In agreement with experimental observations, the simulations predict that the polymer adsorption isotherm goes through a maximum as the calcium or the polymer concentration is increased. The non-Langmuir isotherms are rationalized in terms of charge-charge correlations.

Coarse-Graining Intermolecular Interactions in Dispersions of Highly Charged Colloids

Effective pair potentials between charged colloids, obtained from Monte Carlo simulations of two single colloids in a closed cell at the primitive model level, are shown to reproduce accurately the structure of aqueous salt-free colloidal dispersions, as determined from full primitive model simulations by Linse et al. (Linse, P.; Lobaskin, V. Electrostatic Attraction and Phase Separation in Solutions of Like-Charged Colloidal Particles. Phys. Rev. Lett.1999, 83, 4208). Excellent agreement is obtained even when ion-ion correlations are important and is in principle not limited to spherical particles, providing a potential route to coarse-grained colloidal interactions in more complex systems.

Stability of Negatively Charged Platelets in Calcium-Rich Anionic Copolymer Solutions

Controlling the stability of anisotropic particles is key to the development of advanced materials. Here, we report an investigation, by means of mesoscale molecular dynamics simulations, of the stability and structural change of calcium-rich dispersions containing negatively charged nanoplatelets, neutralized by calcium counterions, in the presence of either comb copolymers composed of anionic backbones with attached neutral side chains or anionic-neutral linear block copolymers. In agreement with experimental observations, small stacks of platelets (tactoids) are formed, which are greatly stabilized in the presence of copolymers. In the absence of polymers, tactoids will grow and aggregate strongly due to large attractive Ca(2+)-Ca(2+) correlation forces. Unlike comb copolymers which only adsorb on the external surfaces, block copolymers are found to intercalate between the platelets. The present results show that the stabilization is the result of a free energy barrier induced by the excluded volume of hydrophilic chains, while the intercalation is due to bridging forces. More generally, the results shed new light on the recent finding of the first hybrid cementitious mesocrystal.

Structure and Yielding of Colloidal Silica Gels Varying the Range of Interparticle Interactions

The relationship between interaction range, structure, fluid-gel transition, and viscoelastic properties of silica dispersions at intermediate volume fraction, Φv ≈ 0.1 and in alkaline conditions, pH = 9 was investigated. For this purpose, rheological, physicochemical, and structural (synchrotron-SAXS) analyses were combined. The range of interaction and the aggregation state of the dispersions were tuned by adding either divalent counterions (Ca(2+)) or polycounterions (PDDA). With increasing calcium chloride concentration, a progressive aggregation was observed which precludes a fluid-gel transition at above 75 mM of calcium chloride. In this case, the aggregation mechanism is driven by short-range ion-ion correlations. Upon addition of PDDA, a fluid-gel transition, at a much lower concentration, followed by a reentrant gel-fluid transition was observed. The gel formation with PDDA was induced by charge neutralization and longer range polymer bridging interactions. The refluidification at high PDDA concentrations was explained by the overcompensation of the charge of the silica particles and by the steric repulsions induced by the polycation chains. Rheological measurements on the so-obtained gels reveal broad yielding transition with two steps when the size of the silica particle clusters exceeds ≈0.5 μm.

Classical density functional theory & simulations on a coarse-grained model of aromatic ionic liquids

A new classical density functional approach is developed to accurately treat a coarse-grained model of room temperature aromatic ionic liquids. Our major innovation is the introduction of charge-charge correlations, which are treated in a simple phenomenological way. We test this theory on a generic coarse-grained model for aromatic RTILs with oligomeric forms for both cations and anions, approximating 1-alkyl-3-methyl imidazoliums and BF₄⁻, respectively. We find that predictions by the new density functional theory for fluid structures at charged surfaces are very accurate, as compared with molecular dynamics simulations, across a range of surface charge densities and lengths of the alkyl chain. Predictions of interactions between charged surfaces are also presented.

Aggregation of Calcium Silicate Hydrate Nanoplatelets

We study the aggregation of calcium silicate hydrate nanoplatelets on a surface by means of Monte Carlo and molecular dynamics simulations at thermodynamic equilibrium. Calcium silicate hydrate (C-S-H) is the main component formed in cement and is responsible for the strength of the material. The hydrate is formed in early cement paste and grows to form platelets on the nanoscale, which aggregate either on dissolving cement particles or on auxiliary particles. The general result is that the experimentally observed variations in these dynamic processes generically called growth can be rationalized from interaction free energies, that is, from pure thermodynamic arguments. We further show that the surface charge density of the particles determines the aggregate structures formed by C-S-H and thus their growth modes.

Interaction and aggregation of charged platelets in electrolyte solutions: a coarse-graining approach.

A coarse-graining approach has been developed to replace the effect of explicit ions with an effective pair potential between charged sites in anisotropic colloidal particles by optimizing a potential of mean force against the results of simulations of two such colloidal particles with all ions in a cell model. More specifically, effective pair potentials were obtained for charged platelets in electrolyte solutions by simulating two rotating parallel platelets with ions at the primitive model level, enclosed in a cylindrical cell. One-component bulk simulations of many platelets interacting via the effective pair potentials are in excellent agreement with the corresponding bulk simulations with all mobile charges present. The bulk simulations were mainly used to study the effects of platelet size, flexibility, and surface charge density on platelet aggregation in an aqueous 2:1 electrolyte, but systems in a 1:1 electrolyte were also investigated.

Block polyelectrolytes and colloidal stability

We simulate interactions between charged flat surfaces in the presence of block polymers, where the end blocks carry a charge opposite to the surfaces. Using a recently developed simulation technique, we allow full equilibrium, i.e. the chemical potential of the polyelectrolyte is retained as the separation is changed. In general, the block polyions will adsorb strongly enough to overcharge the surfaces. This results in a double layer repulsion at large separation, with a concomitant free energy barrier. At short separations, the surfaces are pulled together by bridging forces. We make some efforts to theoretically design the polymers to be efficient flocculants, i.e. minimize the free energy barrier and still allow for a long-ranged bridging attraction. Here, we also take into account the possibility of nonequilibrium circumstances, which may be relevant in practice. Our results suggest that short chains, with small charged end blocks and a (relatively speaking) long neutral mid block, are likely to promote rapid flocculation.

Simulating equilibrium surface forces in polymer solutions using a canonical grid method

A new simulation method for nonuniform polymer solutions between planar surfaces at full chemical equilibrium is described. The technique uses a grid of points in a two-dimensional thermodynamic space, labeled by surface area and surface separations. Free energy differences between these points are determined via Bennetts optimized rates method in the canonical ensemble. Subsequently, loci of constant chemical potential are determined within the grid via simple numerical interpolation. In this way, a series of free energy versus separation curves are determined for a number of different chemical potentials. The method is applied to the case of hard sphere polymers between attractive surfaces, and its veracity is confirmed via comparisons with established alternative simulation techniques, namely, the grand canonical ensemble and isotension ensemble methods. The former method is shown to fail when the degree of polymerization is too large. An interesting interplay between repulsive steric interactions and attractive bridging forces occurs as the surface attraction and bulk monomer density are varied. This behavior is further explored using polymer density functional theory, which is shown to be in good agreement with the simulations. Our results are also discussed in light of recent self-consistent field calculations which correct the original deGennes results for infinitely long polymers. In particular, we look at the role of chain ends by investigating the behavior of ring polymers.

Theoretical predictions of structures in dispersions containing charged colloidal particles and non-adsorbing polymers

We develop a theoretical model to describe structural effects on a specific system of charged colloidal polystyrene particles, upon the addition of non-adsorbing PEG polymers. This system has previously been investigated experimentally, by scattering methods, so we are able to quantitatively compare predicted structure factors with corresponding experimental data. Our aim is to construct a model that is coarse-grained enough to be computationally manageable, yet detailed enough to capture the important physics. To this end, we utilize classical polymer density functional theory, wherein all possible polymer configurations are accounted for, subject to a mean-field Boltzmann weight. We make efforts to counteract drawbacks with this mean-field approach, resulting in structural predictions that agree very well with computationally more demanding simulations. Electrostatic interactions are handled at the fully non-linear Poisson-Boltzmann level, and we demonstrate that a linearization leads to less accurate predictions. The particle charge is an experimentally unknown parameter. We define the surface charge such that the experimental and theoretical gel point at equal polymer concentration coincide. Assuming a fixed surface charge for a certain salt concentration, we find very good agreements between measured and predicted structure factors across a wide range of polymer concentrations. We also present predictions for other structural quantities, such as radial distribution functions, and cluster size distributions. Finally, we demonstrate that our model predicts the occurrence of equilibrium clusters at high polymer concentrations, but low particle volume fractions and salt levels.