Christophe Labbez

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.

Mechanism of aluminium incorporation into C–S–H from ab initio calculations

Blended cements have great potential to reduce the CO2 footprint due to cement production. C(alcium)–S(ilicate)–H(ydrate) in these novel materials is known to incorporate a considerable amount of Al. We have for the first time applied large-scale first principles calculations to address the mechanism of Al incorporation into low C/S ratio C–S–H. In agreement with state-of-the-art NMR information, our calculations show that Al substitutes Si in bridging tetrahedra only, and that substitutions in pairing tetrahedra are strongly disfavoured in a wide range of conditions. In broad terms, the energy penalty for having an Al atom in a pairing position is of about 20 kcal mol−1. Al in bridging tetrahedra is therefore the thermodynamically favoured state, rather than merely a kinetically trapped one in a solid–liquid equilibrium known experimentally to take a very long time to reach. A systematic investigation of Al–Al and defect–Al correlations shows that having two Al atoms as next-neighbours is particularly unfavourable, which gives clues on the limit of Al incorporation into C–S–H. All in all, the current work supports the model and methodology employed to pursue further studies in such materials (e.g., higher C/S ratio systems), in the context of what is still the open question of the structure of C–S–H.

Adsorption of monovalent and divalent cations on planar water-silica interfaces studied by optical reflectivity and Monte Carlo simulations

Adsorption on planar silica substrates of various monovalent and divalent cations from aqueous solution is studied by optical reflectivity. The adsorbed amount is extracted by means of a thin slab model. The experimental data are compared with grand canonical Monte Carlo titration simulations at the primitive model level. The surface excess of charge due to adsorbed cations is found to increase with pH and salt concentration as a result of the progressive dissociation of silanol groups. The simulations predict, in agreement with experiments, that the surface excess of charge from divalent ions is much larger than from monovalent ions. Ion-ion correlations explain quantitatively the enhancement of surface ionization by multivalent cations. On the other hand, the combination of experimental and simulation results strongly suggests the existence of a second ionizable site in the acidic region. Variation of the distance of closest approach between the ions and surface sites captures ion specificity of water-silica interfaces in an approximate fashion.

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.

Monte Carlo simulations of a clay inspired model suspension: the role of rim charge

We present a theoretical investigation of a model clay dispersion in 1-1 salt solutions by varying the particle volume fraction and ionic strength as well as the charge distribution on the clay platelets. The platelets are modeled as discs with charged sites distributed on a hexagonal lattice. The edge sites can be positively charged while the remaining sites are negative giving rise to a strong charge anisotropy. Simulations are carried out using a Monte Carlo method in the canonical ensemble. The interactions between the platelet sites are described with a screened Coulomb potential plus a short range repulsive potential. Simulations show a complex phase behavior. When the charge anisotropy is strong, i.e. all edge sites are positively charged, a fluid phase dominated by repulsion is found at low volume fraction and ionic strength. When increasing the latter an attractive liquid phase forms. At these volume fractions the platelets aggregate in an “Overlapping Coins” configuration. With increasing volume fraction the dispersion becomes unstable and the pressure goes through a van der Waals loop. A liquid crystalline phase, Smectic B, forms in the thermodynamically unstable region. On the other side of the van der Waals loop a stable gel phase is found. A phase separation between a fluid and a gel is thus predicted. The threshold value of the volume fraction at which the phase separation occurs is found to increase with the salt concentration. The gel structure is a mixture of “Overlapping Coins” and “House of Cards” configurations. When the charge anisotropy is intermediate, no phase separation occurs. Instead, a gel forms from a sol of clusters of individual particles randomly oriented that progressively grow with the volume fraction. These results are discussed in light of experimental observations on clay suspensions.

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.

Identification of binding peptides on calcium silicate hydrate: a novel view on cement additives.

Phage display experiments on industrially important calcium silicate hydrates (C-S-H), the main hydration product of ordinary Portland cement, suggest fundamentally different specific binding motifs compared to hitherto existing commercial cement additives. According to that, a strong and specific adsorbing additive on C-S-H should have three features which are a negative charge, H-bond formers (especially amide functions) and a hydrophobic part.

Liquid Crystal Phases in Suspensions of Charged Plate-Like Particles.

Anisotropic interactions in colloidal suspensions have recently emerged as a route for the design of new soft materials. Nonisotropic particles can form nematic, smectic, hexatic, and columnar liquid crystals. Although the formation of these phases is well rationalized when excluded volume is solely at play, the role of electrostatic interactions still remains unclear and even less so when particles present a charge heterogeneity, for example, clays. Here, we use Monte Carlo simulations of concentrated suspensions of charged disk-like particles to reveal the role of Coulomb interactions and charge anisotropy underlying liquid crystal formation and structures. We observe a vast zoo of exotic structures, going from hexatic to columnar phases, which are shown to be controlled by the charge anisotropy. The particle volume fraction at which these phases start to form is found to decrease with increasing Coulomb interactions and charge anisotropy, which suggests a route to tune the structure of aqueous liquid crystals.

Monte carlo simulations of parallel charged platelets as an approach to tactoid formation in clay.

The free energy of interaction between parallel charged platelets with divalent counterions has been calculated using Monte Carlo simulations to investigate the electrostatic effects on aggregation. The platelets are primarily intended to represent clay particles. With divalent counterions, the free energy for two platelets or two tactoids (clusters of parallel platelets) shows a minimum at a short separation due to the attraction caused by ion-ion correlations. In a salt-free system, the free energy of interaction has a long-range repulsive tail beyond the minimum. The repulsion increases for tactoids with larger aggregation numbers, whereas the depth of the free-energy minimum is gradually reduced. For large enough aggregation numbers, the repulsion is dominating and the minimum is no longer a global free-energy minimum. This is an effect of the depletion of counterions free in solution (outside tactoids) as counterions and platelets aggregate into tactoids and the resulting redistribution of counterions in the system changes the effective interactions between platelets and tactoids. The difference in tactoid-tactoid interactions as a function of aggregation number can be removed by adding enough salt to mask the depletion. Adding salt also reduces the repulsive tail of the free energy of interaction and enhances the minimum. No dependence on the aggregation number suggests that an isodesmic model with a monotonically decaying distribution of aggregation numbers can be used to describe a clay system. This may help to explain the experimental observations of low average numbers of platelets in tactoids, although factors not included in the simulation model may also play an important role.

A new Monte Carlo method for the titration of molecules and minerals

The charge state of molecules and solid/liquid interfaces is of paramount importance in the understanding of the reactivity and the physico-chemical properties of many systems. In this work, we porpose a new Monte Carlo method in the grand canonical ensemble using the primitive model, which allows us to simulate the titration behavior of macromolecules or solids at constant pH. The method is applied to the charging process of colloidal silica particles dispersed in a sodium salt solution for various concentrations and calcium silicate hydrate nano-particles in a calcium hydroxide solution. An excellent agreement is found between the experimental and simulated results.