Acid-Base Properties of 2:1 Clays. I. Modeling the Role of Electrostatics

We present a theoretical investigation of the titratable charge of clays with various structural charge (sigma(b)): pyrophyllite (sigma(b) = 0 e x nm(-2)), montmorillonite (sigma(b) = -0.7 e x nm(-2)) and illite (sigma(b) = -1.2 e x nm(-2)). The calculations were carried out using a Monte Carlo method in the Grand Canonical ensemble and in the framework of the primitive model. The clay particle was modeled as a perfect hexagonal platelet, with an "ideal" crystal structure. The only fitting parameters used are the intrinsic equilibrium constants (pK(0)) for the protonation/deprotonation reactions of the broken-bond sites on the lateral faces of the clay particles, silanol, =SiO(-) + H(+) --=…