Molecular simulation of a model of dissolved organic matter

A series of atomistic simulations was performed to assess the ability of the Schulten dissolved organic matter (DOM) molecule, a well-established model humic molecule, to reproduce the physical and chemical behavior of natural humic substances. The unhydrated DOM molecule had a bulk density value appropriate to humic matter, but its Hildebrand solubility parameter was lower than the range of current experimental estimates. Under hydrated conditions, the DOM molecule went through conformational adjustments that resulted in disruption of intramolecular hydrogen bonds (H-bonds), although few water molecules penetrated the organic interior. The radius of gyration of the hydrated DOM molecule was similar to those measured for aquatic humic substances. To simulate humic materials under aqueous conditions with varying pH levels, carboxyl groups were deprotonated, and hydrated Na+ or Ca2+ were added to balance the resulting negative charge. Because of intrusion of the cation hydrates, the model metal-humic structures were more porous, had greater solvent-accessible surface areas, and formed more H-bonds with water than the protonated, hydrated DOM molecule. Relative to Na+, Ca2+ was both more strongly bound to carboxylate groups and more fully hydrated. This difference was attributed to the higher charge of the divalent cation. The Ca-DOM hydrate, however, featured fewer H-bonds than the Na-DOM hydrate, perhaps because of the reduced orientational freedom of organic moieties and water molecules imposed by Ca2+. The present work is, to our knowledge, the first rigorous computational exploration regarding the behavior of a model humic molecule under a range of physical conditions typical of soil and water systems.