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

Mesoscale Simulations Reveal How Salt Influences Clay Particles Agglomeration in Aqueous Dispersions

The aggregation of clay particles is an everyday phenomenon of scientific and industrial relevance. However, it is a complex multiscale process that depends delicately on the nature of the particle-particle and particle-solvent interactions. Toward understanding how to control such phenomena, a multiscale computational approach is developed, building from molecular simulations conducted at atomic resolution to calculate the potential of mean force (PMF) profiles in both pure and saline water environments. We document how it is possible to use such a model to develop a fundamental understanding concerning the mechanism of particle aggregation. For example, using molecular dynamics simulations conducted at the mesoscale in implicit solvents, it is possible to quantify the size and shape of clay aggregates as a function of system conditions. The approach is used to emphasize the role of salt concentration, which directly affects the potentials of the mean forces between kaolinite particles. While particle agglomeration in pure water yields large aggregates, the presence of sodium chloride in the aqueous brine leads instead to a large number of small aggregates. These results are consistent with macroscopic experimental observations, suggesting that the simulation protocol developed could be relevant for preventing pore blocking in heterogeneous porous matrixes.

Aggregation of Plate-like Colloids Induced by Charged Polymer Chains: Organization at the Nanometer Scale Tuned by Polymer Charge Density

We study the aggregation of charged plate-like colloids, Na-montmorillonite clays, in the presence of ionenes, oppositely charged polymer chains. The choice of the charged polymer allows tuning its linear charge density to match/mismatch the average charge separation on the clay surfaces. We assess the nanoscale structure of the aggregates formed by small-angle X-ray and neutron scattering. The nanoscale features of the formed clay aggregates are dominated by the presence of a stacking peak, giving clear evidence for the formation of clay tactoids, that is, a face-to-face aggregation geometry of the clay platelets. The chain charge density of ionenes influences not only the stacking repeat distance within the clay tactoids but also the extent of stacking and abundance of the tactoids. We may distinguish two regimes as a function of clay and ionene polymer charge densities (ρ_c and ρ_p, respectively). The first regime applies to ρ_p > ρ_c and ρ_p ≈ ρ_c, that is, for highly and "matching" charged chains. Under these conditions, the intercalated chains lie in a flat conformation within the tactoids, irrespective of the ionic strength (within the range studied, i.e., up to 0.05 M NaBr). For weakly charged chains, ρ_p < ρ_c, undulation of the ionene chains within the tactoid is seen. The degree of undulation increases with ionic strength due to the decreasing persistence length of the ionene chains. The extent of stacking (5-10 platelets per tactoid) is a general feature of all the systems, and its origin remains unknown. The system corresponding to the closest match in charge separations on the clay surface and on the polymer chain (ρ_p ≈ ρ_c) features the highest abundance of tactoids. This coincides with the highest macroscopic density as deduced from simple visual inspection of sediment volumes. This leads to the open question regarding the link between the density at the nanoscale and the macroscopic density and sedimentation behavior of the aggregate.

Microstructure and Soft Glassy Dynamics of an Aqueous Laponite Dispersion

Synthetic hectorite clay Laponite RD/XLG is composed of disk-shaped nanoparticles that acquire dissimilar charges when suspended in an aqueous medium. Owing to their property to spontaneously self-assemble, Laponite is used as a rheology modifier in a variety of commercial water-based products. In particular, an aqueous dispersion of Laponite undergoes a liquid-to-solid transition at about 1 vol % concentration. The evolution of the physical properties as the dispersion transforms to the solid state is reminiscent of physical aging in molecular as well as colloidal glasses. The corresponding soft glassy dynamics of an aqueous Laponite dispersion, including the rheological behavior, has been extensively studied in the literature. In this feature article, we take an overview of recent advances in understanding soft glassy dynamics and various efforts taken to understand the peculiar rheological behavior. Furthermore, the continuously developing microstructure that is responsible for the eventual formation of a soft solid state that supports its own weight against gravity has also been a topic of intense debate and discussion. In particularly, extensive experimental and theoretical studies lead to two types of microstructures for this system: an attractive gel-like or a repulsive glass-like structure. We carefully examine and critically analyze the literature and propose a state (phase) diagram that suggests an aqueous Laponite dispersion to be present in an attractive gel state.

The effect of the relative permittivity on the tactoid formation in nanoplatelet systems. A combined computer simulation, SAXS, and osmotic pressure study

The structural properties, and the intracrystalline swelling of Na⁺-, and Ca²⁺-montmorillonite (Na-, and Ca-mmt) have been investigated as an effect of decreasing the relative permittivity of the solvent, i.e. from water to ethanol (EtOH), utilizing the experimental techniques; small angle X-ray scattering (SAXS) and osmotic pressure measurements. The experimental data were compared with the continuum model, utilizing coarse-grained molecular dynamics bulk simulations, Monte Carlo simulations of two infinite parallel surfaces corresponding to two clay platelets, and the strong coupling theory. It was found that it is possible to tune the electrostatic interactions to obtain a transition from a repulsive to an attractive system for the Na-mmt by increasing the EtOH concentration, i.e. the Bjerrum length increases, and hence, the attractive ion-ion correlation forces are enhanced. A qualitative agreement was observed between the simulations and the experimental results. Moreover, a non-monotonic behavior of the intracrystalline swelling of Ca-mmt as a function of EtOH concentration was captured experimentally, where an increase in the osmotic pressure, and hence, an increase in the d-spacing was found at low concentrations, indicating that repulsive short-ranged interactions dominate in the system. Theoretically, the non-monotonic behavior could not be captured with the continuum model, probably due to the limitation that the electrostatic interactions solely enters the Hamiltonian via the Bjerrum length.

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 grand potential between calcium-silicate-hydrate nanoparticles at the molecular level

Calcium-silicate-hydrate (or C-S-H), an inosilicate, is the major binding phase in cement pastes and concretes and a porous hydrated material made up of a percolated and dense network of crystalline nanoparticles of a mean apparent spherical diameter of ∼5 nm that are each stacks of multiple C-S-H layers. Interaction forces between these nanoparticles are at the origin of C-S-H chemical, physical, and mechanical properties at the meso- and macroscales. These particle interactions and the resulting properties may be affected significantly by nanoparticle density and environmental conditions such as the temperature, relative humidity, or concentration of chemical species in the bulk solution. In this study, we combined grand canonical Monte Carlo simulations and an extension of the mean force integration method to derive the pair potentials. This approach enables realistic simulation of the physical environment surrounding the C-S-H particles. We thus constructed the pair potentials for C-S-H nanoparticles of defined chemical stoichiometry at 10% relative humidity (RH), varying the relative crystallographic orientations at a constant particle density of ρpart ∼ 2.21 mmol L(-1). We found that cohesion between nanoparticles is affected strongly by both the aspect ratio and the crystallographic misorientation of interacting particles. This method and the findings underscore the importance of accounting for relative dimensions and orientation among C-S-H nanoparticles in descriptions of physical and simulated multiparticle aggregates or mesoscale systems.

Reactivity, swelling and aggregation of mixed-size silicate nanoplatelets

Montmorillonite is a key ingredient in a number of technical applications. However, little is known regarding the microstructure and the forces between silicate platelets. The size of montmorillonite platelets from different natural sources can vary significantly. This has an influence on their swelling behavior in water as well as in salt solutions, particularly when tactoid formation occurs, that is when divalent counterions are present in the system. A tactoid consists of a limited number of platelets aggregated in a parallel arrangement with a constant separation. The tactoid size increases with platelet size and with very small nanoplatelets, ∼30 nm, no tactoids are observed irrespectively of the platelet origin and concentration of divalent ions. The formation and dissociation of tactoids seem to be reversible processes. A large proportion of small nanoplatelets in a mixed-size system affects the tactoid formation, reduces the aggregation number and increases the extra-lamellar swelling in the system.

Thermodynamics of the self-assembly of non-ionic chromonic molecules using atomistic simulations. The case of TP6EO2M in aqueous solution

Atomistic molecular dynamic simulations have been performed for the non-ionic chromonic liquid crystal 2,3,6,7,10,11-hexa-(1,4,7-trioxa-octyl)-triphenylene (TP6EO2M) in aqueous solution. TP6EO2M molecules consist of a central poly-aromatic core (a triphenylene ring) functionalized by six hydrophilic ethyleneoxy (EO) chains, and have a strong tendency to aggregate face-to-face into stacks even in very dilute solution. We have studied self-assembly of the molecules in the low concentration range corresponding to an isotropic solution of aggregates, using two force fields GAFF and OPLS. Our results reveal that the GAFF force field, even though it was successfully used previously for modelling of ionic chromonics, overestimates the attraction of TP6EO2M molecules in water. This results in an aggregation free energy which is too high, a reduced hydration of EO chains and, therefore, molecular self-assembly into compact disordered clusters instead of stacks. In contrast, use of the OPLS force field, leads to self-assembly into ordered stacks in agreement with earlier experimental studies of triphenylene-based chromonics. The free energy of association follows a "quasi-isodesmic" pattern, where the binding free energy of two molecules to form a dimer is of the order of 2.5 RT larger than the corresponding energy of addition of a molecule into a stack. The obtained value for the binding free energy, ΔG=-12 RT, is found to be in line with the published values for typical ionic chromonics (-7 to -12 RT), and agrees reasonably well with the experimental results for this system. The calculated interlayer distance between the molecules in a stack is 0.37 nm, which is at the top of the range found for typical chromonics (0.33-0.37 nm). We suggest that the relatively large layer spacing can be attributed to the repulsion between EO side chains.

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.