Geoffrey M. Bowers

Structural and dynamical relationships of Ca2+ and H2O in smectite/2H2O systems

Abstract We present an X-ray diffraction and multi-nuclear (2H and 43Ca) NMR study of Ca-exchanged hectorite (a smectite clay) that provides important new insight into molecular behavior at the smectite-H2O interface. Variable-temperature 43Ca MAS NMR and controlled humidity XRD indicate that Ca2+ occurs as proximity-restricted outer-sphere hydration complexes between -120 and +25 °C in a two-layer hydrate and at T ≤ -50 °C in a 2:1 water/solid paste. Changes in the 43Ca NMR peak width and position with temperature are more consistent with diffusion-related processes than with dynamics involving metal-surface interactions such as site exchange. The 2H NMR signal between -50 and +25 °C for a two-layer hydrate of Ca-hectorite is similar to that of Na- and other alkali metal hectorites and represents 2H2O molecules experiencing anisotropic motion describable using the 2H C2/C3 jump model we proposed previously. 2H T1 relaxation results for Ca- and Na-hectorite are well fit with a fast-exchange limit, rotational diffusion model for 2H2O dynamics, yielding GHz-scale rotational reorientation rates compatible with the C3 component of the C2/C3 hopping model. The apparent activation energy for 2H2O rotational diffusion in the two-layer hydrate is greater for Ca-hectorite than Na-hectorite (25.1 vs. 21.1 kJ/mol), consistent with the greater affinity of Ca2+ for H2O. The results support the general principle that the dynamic mechanisms of proximity-restricted H2O are only weakly influenced by the cation in alkali metal and alkaline earth metal smectites and provide critical evidence that the NMR resonances of charge-balancing cations in smectites become increasingly influenced by diffusion-like dynamic processes at low temperatures as the charge density of the unhydrated cation increases.

NMR and computational molecular modeling studies of mineral surfaces and interlayer galleries: A review.

Abstract This paper reviews experimental nuclear magnetic resonance (NMR) and computational molecular dynamics (MD) investigations of the structural and dynamical behavior of cations, anions, H2O, and CO2 on the surfaces and in the interlayer galleries of layer-structure minerals and their composites with polymers and natural organic matter (NOM). The interaction among mineral surfaces, chargebalancing cations or anions, H2O, CO2, and NOM are dominated by Coulombic, H-bond, and van der Waals interactions leading to statically and dynamically disordered systems and molecular-scale processes with characteristic room-temperature frequencies varying from at least as small as 102 to >1012 Hz. NMR spectroscopy provides local structural information about such systems through the chemical shift and quadrupolar interactions and dynamical information at frequencies from the subkilohertz to gigahertz ranges through the T1 and T2 relaxation rates and line shape analysis. It is often difficult to associate a specific structure or dynamical process to a given NMR observation, however, and computational molecular modeling is often effective in providing a much more detailed picture in this regard. The examples discussed here illustrate these capabilities of combining experimental NMR and computational modeling in mineralogically and geochemically important systems, including clay minerals and layered double hydroxides.

Molecular dynamics modeling of carbon dioxide, water and natural organic matter in Na-hectorite

Molecular dynamics (MD) modeling of systems containing a Na-exchanged smectite clay (hectorite) and model natural organic matter (NOM) molecules along with pure H2O, pure CO2, or a mixture of H2O and CO2 provides significant new insight into the molecular scale interactions among silicate surfaces, dissolved cations and organic molecules, H2O and CO2 relevant to geological C-sequestration strategies. The simulations for systems containing H2O show the following results; (1) Na(+) does not bridge between NOM molecules and the clay surface at protonation states comparable to near neutral pH conditions. (2) In systems without CO2 the NOM molecules retain charge balancing cations and drift away from the silicate surface. (3) In systems containing both H2O and CO2, the NOM molecules adopt equilibrium positions at the H2O-CO2 interface with the more hydrophilic structural elements facing the H2O and the more hydrophobic ones facing the CO2. In systems with only CO2, NOM and Na(+) ions are pinned to the clay surface with the hydrophilic structural elements of the NOM pointed toward the clay surface. Dynamically, in systems with only CO2, Na(+) diffusion is nearly eliminated, and in systems with a thin water film on the clay surface diffusion perpendicular the surface is greatly reduced relative to the system with bulk water. Energetically, the results for the systems with only H2O show that hydration of the net charge neutral Na-NOM molecule outweighs the sum of its Coulombic and dispersive interactions with the net charge-neutral Na-clay particle and the interactions of the water molecules with the hydrophobic structural elements of the NOM. The aggregation of NOM molecules in solution appears to be driven not by Na(+) bridging between the molecules but by hydrophobic interactions between them. In contrast, for the systems with only CO2 the interaction between the Na-NOM molecules and the CO2 is outweighed by the interaction of NOM with the clay particle. With both H2O and CO2 present, the energetic interactions leading to the hydration of the Na-clay surface and the hydrophilic structural elements of the Na-NOM molecule and the hydrophobic interactions between the CO2 and the hydrophobic aromatic and aliphatic structural elements of the NOM can both be satisfied, leading to the Na-NOM molecules migrating away from the surface and residing at the H2O-CO2 interface. The MD results suggest some alternative explanations for the previously observed (23)Na NMR behavior of Na-hectorite at elevated temperatures and CO2 pressures.

Tipping Point for Expansion of Layered Aluminosilicates in Weakly Polar Solvents: Supercritical CO2

Layered aluminosilicates play a dominant role in the mechanical and gas storage properties of the subsurface, are used in diverse industrial applications, and serve as model materials for understanding solvent-ion-support systems. Although expansion in the presence of H2O is well-known to be systematically correlated with the hydration free energy of the interlayer cation, particularly in environments dominated by nonpolar solvents (i.e., CO2), uptake into the interlayer is not well-understood. Using novel high-pressure capabilities, we investigated the interaction of dry supercritical CO2 with Na-, NH4-, and Cs-saturated montmorillonite, comparing results with predictions from molecular dynamics simulations. Despite the known trend in H2O and that cation solvation energies in CO2 suggest a stronger interaction with Na, both the NH4- and Cs-clays readily absorbed CO2 and expanded, while the Na-clay did not. The apparent inertness of the Na-clay was not due to kinetics, as experiments seeking a stable expanded state showed that none exists. Molecular dynamics simulations revealed a large endothermicity to CO2 intercalation in the Na-clay but little or no energy barrier for the NH4- and Cs-clays. Indeed, the combination of experiment and theory clearly demonstrate that CO2 intercalation of Na-montmorillonite clays is prohibited in the absence of H2O. Consequently, we have shown for the first time that in the presence of a low dielectric constant, gas swelling depends more on the strength of the interaction between the interlayer cation and aluminosilicate sheets and less on that with solvent. The finding suggests a distinct regime in layered aluminosilicate swelling behavior triggered by low solvent polarizability, with important implications in geomechanics, storage, and retention of volatile gases, and across industrial uses in gelling, decoloring, heterogeneous catalysis, and semipermeable reactive barriers.

Collboration: Interfacial Soil Chemistry of Radionuclides in the Unsaturated Zone

The principal goal of this project was to assess the molecular nature and stability of radionuclide immoblization during weathering reactions in bulk Hanford sediments and their high surface area clay mineral constituents. We focused on the unique aqueous geochemical conditions that are representative of waste-impacted locations in the Hanford site vadose zone; high ionic strength, high pH and high Al concentrations. The specific objectives of the work were to measure the coupling of clay mineral weathering and contaminant uptake kinetics of Cs, Sr and I; determine the molecular structure of contaminant binding sites and their change with weathering time during and after exposure to synthetic tank waste leachate; establish the stability of neoformed weathering products and their sequestered contaminbants upon exposure of the solids to more natural soil solutaions afer remofal of the caustic waste source; and integrate macroscopic, microscopic and spectroscopic data to distinguish labile from non-labile contaminant binding environments, including their dependence on system composition and weathering time.