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.