Molecular dynamics simulations of sorption of organic compounds at the clay mineral/aqueous solution interface

Abstract

The adsorption of trichloroethene, C2HCl3, on clay mineral surfaces in the presence of water has been modeled as an example describing a general program that uses molecular dynamics simulations to study the sorption of organic materials at the clay mineral/aqueous solution interface. Surfaces of the clay minerals kaolinite and pyrophyllite were hydrated at different water levels corresponding to partial and complete monolayers of water. In agreement with experimental trends, water was found to outcompete C2HCl3 for clay surface sites. The simulations suggest that at least three distinct mechanisms coexist for C2HCl3 on clay minerals in the environment. The most stable interaction of C2HCl3 with clay surfaces is by full molecular contact, coplanar with the basal surface. This kind of interaction is suppressed by increasing water loads. A second less stable and more reversible interaction involves adsorption through single-atom contact between one Cl atom and the surface. In a third mechanism, adsorbed C2HCl3 never contacts the clay directly but sorbs onto the first water layer. To test the efficacy of existing force field parameters of organic compounds in solid state simulations, molecular dynamics simulations of several representative organic crystals were also performed and compared with the experimental crystal structures. These investigations show that, in general, in condensed-phase studies, parameter evaluations are realistic only when thermal motion effects are included in the simulations. For chlorohydrocarbons in particular, further explorations are needed of atomic point charge assignments. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 144–153, 1998