We report a Monte Carlo and molecular dynamics simulations study of carbon dioxide in hydrated sodium montmorillonite, including thermodynamical, structural and dynamical properties. In order to simulate the behaviour of a clay caprock in contact with a CO 2 reservoir, we consider clays in equilibrium with H 2 O−CO 2 mixtures under conditions close to relevant ones for geological storage, namely a temperature T =348 K, and pressures P=25 and 125 bar, and under which two bulk phases coexist: H 2 O-rich liquid on the one hand and CO 2-rich gas (P=25 bar) or supercritical fluid (P=125 bar) on the other hand. We first use grand-canonical MC simulations to determine the number of stable states in clay, their composition and the corresponding equilibrium interlayer distances. The vertical, horizontal and radial distribution functions of the confined mixture, subsequently obtained using molecular dynamics, reveal some structural feature induced by the presence of CO 2. Finally, the simulations indicate that carbon dioxide considerably influences the diffusion of mobile species in clays. We discuss these results by comparing them with those obtained for the bulk mixtures, as well as for Namontmorillonite in equilibrium with a pure water reservoir water at the same temperature and pressure.

Departures from ideality in electrolytes are described in the framework of the primitive model of ionic solutions in which the solvent is a dielectric continuum, using the mean spherical approximation (MSA). To include solvation and solvent concentration effects, we consider that the permittivity of the solvent and the sizes of the ions are concentration-dependent parameters. New expressions are derived for the activity coefficients and the osmotic coefficient. They are applied to pure ionic aqueous solutions of 18 salts, taking simple functions for the adjusted parameters. Good fittings are obtained in the concentration range 0−6 mol/kg.

We report a molecular dynamics study of the structure and dynamics of water at a clay surface. The negative charge of the surface and the presence of surface oxygen atoms perturbs water over two to three molecular layers, while the nature of the counterions (Na(+)or Cs(+)) has only a small effect. In the first molecular layer, approximately half of the water molecules are H-bonded to the surface. We also analyze the H-bond network between surface water molecules. The diffusion of water molecules along the surface is slowed down compared to the bulk case. As far as the orientational order and dynamics of the water dipole are concerned, only the component normal to the clay surface is perturbed. We investigate the surface H-bond formation and dissociation dynamics and their coupling to the release of molecules from the first molecular layer. We introduce a simple kinetic model in the spirit of Luzar and Chandler [Nature, 1996, 379, 55] to allow for a comparison with bulk water dynamics. This model semi-quantitatively reproduces the molecular simulation results and suggests that H-bond formation is faster with the surface than in the bulk, while H-bond dissociation is slower.