Adsorption isotherms of water in Li-, Na-, and K-montmorillonite by molecular simulation

A biased Monte Carlo method for the insertion of water in dense clay-water systems is presented. The use of this algorithm results in a considerable i ncrease of the success rate of insertion attempts. It allows us to compute water adsorption isotherms up to high water densities, where the convention al Monte Carlo scheme fails. The isotherms were calculated by a combination of molecular dynamics and grand-canonical Monte Carlo simulation for Li-, Na-, and K-montmorillonite at a fixed d(001) spacing of 12.0 Angstrom. At l ow water pressure, the degree of clay hydration is governed by the type of counterion, Li-montmorillonite having the highest water content. Hydrogen b onding between water molecules is absent. Li+ and Na+ are small enough to b e organized in two layers close to the clay mineral surfaces, whereas K+ is mainly located in the midplane. In both cases, the water molecules primari ly reside in the midplane of the interlayer. Increasing the water pressure leads to water adsorption at higher energy sites closer to the surface, i.e ., coordinating to the structural OH groups in the hexagonal cavities. A hy drogen bond network is formed in the clay interlayer. This points to water condensation and leads to a sharp increase in the clay water content. (C) 2 001 American Institute of Physics.