SIMULATION OF ADSORPTION-ISOTHERMS OF WATER ON IONIC SURFACES

The intention of this work is to aid the construction of a water model capable of describing bulk water in contact with ionic surfaces. The model should be sufficiently simple to be useful for simulating large systems. Here, the extended simple point charge (SPC/E) water model, w hich is one of the most commonly used descriptions of bulk water, was tested in a strongly heterogeneous system. Simulations were carried ou t on the NaCl(100) and the calcite(10 (1) over bar 4) surfaces, descri bed in terms of Lennard-Jones and Coulomb potentials. On the NaCl surf ace Molecular Dynamics simulations show that SPC/E water forms cluster s due to its high permanent dipole moment. Better agreement with the e xperimental adsorbate structure was achieved by reducing the water dip ole moment to its gas phase value. On the calcite surface adsorption i sotherms were simulated using a Molecular Dynamics-Monte Carlo hybrid technique, which consists of a normal Molecular Dynamic simulation wit h randomly chosen insertion and removal Monte Carlo steps repeated aft er a constant time period. The method is tested by calculating adsorpt ion isotherms for two Lennard-Jones systems, i.e. methane in contact w ith the graphite basal plane and with the zeolite silicalite (ZSM5), f or which we obtain excellent agreement with the experimental isotherms . For calcite, however, the simulated adsorption isotherms again show a strong dependence on the choice of the dipole moment. The SPC/E valu e yields a strongly exaggerated coverage, whereas the gas phase value severely underestimates the coverage in comparison with the experiment . From our results we draw the conclusion that a fixed charged model i s not recommendable for the simulation on ionic surfaces, and that it is highly desirable to include variable polarization.