A critical examination of published data obtained primarily from recent nuclear magnetic resonance spectroscopy, X-ray absorption near-edge structure spectroscopy, electrospray ionization-mass spectrometry, and pyrolysis studies reveals an evolving new view of the molecular structure of soil humic substances. According to the new view, humic substances are collections of diverse, relatively low molecular mass components forming dynamic associations stabilized by hydrophobic interactions and hydrogen bonds. These associations are capable of organizing into micellar structures in suitable aqueous environments. Humic components display contrasting molecular motional behavior and may be spatially segregated on a scale of nanometers. Within this new structural context, these components comprise any molecules intimately associated with a humic substance, such that they cannot be separated effectively by chemical or physical methods. Thus biomolecules strongly bound within humic fractions are by definition humic components, a conclusion that necessarily calls into question key biogeochemical pathways traditionally thought to be required for the formation of humic substances. Further research is needed to elucidate the intermolecular interactions that link humic components into supramolecular associations and to establish the pathways by which these associations emerge from the degradation of organic litter.

In this paper, we propose a structure for organo-mineral associations in soils based on recent insights concerning the molecular structure of soil organic matter (SOM), and on extensive published evidence from empirical studies of organo-mineral interfaces. Our conceptual model assumes that SOM consists of a heterogeneous mixture of compounds that display a range of amphiphilic or surfactant-like properties, and are capable of self-organization in aqueous solution. An extension of this self-organizational behavior in solution, we suggest that SOM sorbs to mineral surfaces in a discrete zonal sequence. In the contact zone, the formation of particularly strong organo-mineral associations appears to be favored by situations where either (i) polar organic functional groups of amphiphiles interact via ligand exchange with singly coordinated mineral hydroxyls, forming stable inner-sphere complexes, or (ii) proteinaceous materials unfold upon adsorption, thus increasing adhesive strength by adding hydrophobic interactions to electrostatic binding. Entropic considerations dictate that exposed hydrophobic portions of amphiphilic molecules adsorbed directly to mineral surfaces be shielded from the polar aqueous phase through association with hydrophobic moi

Clay minerals are layer type aluminosilicates that figure in terrestrial biogeochemical cycles, in the buffering capacity of the oceans, and in the containment of toxic waste materials. They are also used as lubricants in petroleum extraction and as industrial catalysts for the synthesis of many organic compounds. These applications derive fundamentally from the colloidal size and permanent structural charge of clay mineral particles, which endow them with significant surface reactivity. Unraveling the surface geochemistry of hydrated clay minerals is an abiding, if difficult, topic in earth sciences research. Recent experimental and computational studies that take advantage of new methodologies and basic insights derived from the study of concentrated ionic solutions have begun to clarify the structure of electrical double layers formed on hydrated clay mineral surfaces, particularly those in the interlayer region of swelling 2:1 layer type clay minerals. One emerging trend is that the coordination of interlayer cations with water molecules and clay mineral surface oxygens is governed largely by cation size and charge, similarly to a concentrated ionic solution, but the location of structural charge within a clay layer and the existence of hydrophobic patches on its surface provide important modulations. The larger the interlayer cation, the greater the inf luence of clay mineral structure and hydrophobicity on the configurations of adsorbed water molecules. This picture extends readily to hydrophobic molecules adsorbed within an interlayer region, with important implications for clay-hydrocarbon interactions and the design of catalysts for organic synthesis.

No abstract

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.

No abstract

Abstract--Monte Carlo simulations based on tested water-water, cation-water, and water-clay potential functions were applied to calculate radial distribution functions for O-O, O-H and H-H spatial correlations in the interlayer region of the two-layer hydrates of Na-and K-montmorillonite. The simulated radial distribution functions then were used to compute the total radial distribution function for interlayer water, a physical quantity that can be determined experimentally by H/D isotopic-difference neutron diffraction. The simulated total radial distribution functions were compared with that for bulk liquid water, and with a total radial distribution function determined experimentally for the two-layer hydrate of Na-montmorillonite by Powell et al. (1997). This comparison suggested that water molecules in the two-layer hydrate of montmorillonite have nearest-neighbor configurations which differ significantly from the tetrahedral ordering of nearest neighbors that characterizes bulk liquid water.

Computer animation of center of mass coordinates obtained from 800 ps molecular dynamics simulations of Cs-smectite hydrates (1/3 and 2/3 water monolayers) provided information concerning the structure and dynamics of the interlayer region that could not be obtained through traditional simulation analysis methods. Cs 1 formed inner sphere complexes with the mineral surface, and could be seen to jump from one attracting location near a layer charge site to the next, while water molecules were observed to migrate from the hydration shell of one ion to that of another. Neighboring ions maintained a partial hydration shell by sharing water molecules, such that a single water molecule hydrated two ions simultaneously for hundreds of picoseconds. Cs-montmorillonite hydrates featured the largest extent of this sharing interaction, because interlayer ions were able to inhabit positions near surface cavities as well as at their edges, close to oxygen triads. The greater positional freedom of Cs 1 within the montmorillonite interlayer, a result of structural hydroxyl orientation and low tetrahedral charge, promoted the optimization of distances between cations and water molecules required for water sharing. Preference of Cs 1 for locations near oxygen triads was observed within interlayer beidellite and hectorite. Water molecules also could be seen to interact directly with the mineral surface, entering its surface cavities to approach attracting charge sites and structural hydroxyls. With increasing water content, water molecules exhibited increased frequency and duration of both cavity habitation and water sharing interactions. Competition between Cs 1 and water molecules for surface sites was evident. These important cooperative and competitive features of interlayer molecular behavior were uniquely revealed by animation of an otherwise highly complex simulation output.