Non-bonded organo-mineral interactions and sorption of organic compounds on soil surfaces: a model approach

Abstract

Molecular mechanics and molecular dynamics calculations have been performed on organo-mineral composites that model the sorption of high-molecular-weight humic polymers on mineral surfaces and the sorption of low-molecular-weight organic contaminants on both mineral and organic surfaces in soil. Muscovite mica was chosen as a mineral model; an oxidized topological lignin-carbohydrate complex was chosen as a humic model; benzene, sodium benzoate, atrazine, and DDT represent different classes of contaminants. Sorption energies were estimated based on molecular mechanics calculations. Flexible linear polymers undergo drastic conformational changes when approaching the mineral surface, to ensure a gain in the interaction energy that outweighs a loss in the conformational energy proper. Therefore, the gas-phase conformation composi tion of environmental organic polymers is not directly related to their spatial organization in soil composites. Molecular dynamics simulation suggests high stability of the organic polymer coatings of mineral surfaces in the environment. Low-molecular-weight organic molecules demonstrate much less affinity for the mineral surface, which implies unhindered exchanges between the surface and its near environment. Ionizable compounds, e.g. salts of organic acids, are different, because they can form strong associations with a mineral surface through cation bridges. Sorption energies are compound-specific and depend on the sorbate-sorbent orientation. The calculations suggest some preference for the edges of a model muscovite sheet in comparison with the basal oxygen surface as a sorption site. Coating of mineral surfaces with organic polymers does not hinder the sorption of organic molecules except in the special case of organic ions.