Molecular dynamics simulation of methane in sodium montmorillonite clay hydrates at elevated pressures and temperatures

Computer simulation has been used to study the structure and dynamics of methane in hydrated sodium montmorillonite clays under conditions encountered in sedimentary basins. Systems containing approximately one, two, three and four molecular layers of water have followed gradients of 150 bar km-1 and 30Kkm-1, to a maximum burial depth of 6 km (900 bar and 460 K). Methane is coordinated to approximately 19 oxygen atoms, of which typically 6 are provided by the clay surface. Only in the three- and four-layer hydrates is methane able to leave the clay surface. Diffusion depends strongly on the porosity (water content) and burial depth: self-diffusion coefficients are in the range 0.12 × 10-9m2s-1 for water and 0.04 × 10−9m2s−1 < D < 8.64 × 10−9m2s−1 for methane. Bearing in mind that porosity decreases with burial depth, it is estimated that maximum diffusion occurs at around 3 km. This is in good agreement with the known location of methane reservoirs in sedimentary basins.

Studies in clay chemistry

A study of clay chemistry has been approached with three aims: - to modify the conducting properties by intercalation of tetrathiafulvalene, - to study the electrochemistry of redox-active coordination compounds immobilised on clay coated electrodes, and - to study the role of clays as reagents in inorganic glass forming reactions using mainly solid-state magic-angle-spinning NMR. TTF was intercalated by smectites containing different interlayer and lattice cations. Evidence from ESR and 57Fe Mossbauer indicated charge-transfer from TTF to structural iron in natural montmorillonite, and to interlayer Cu2+ in Cu2+ exchanged laponite. No charge transfer was observed for laponite (Na+ form) itself. Ion exchange of TTF3(BF4)2 with laponite was found to proceed quantitatively. The intercalated species were believed to be (TTF)2+ dimers. Conductivity data showed an order of magnitude increase for the intercalated clays. The mechanism is thought to be ionic rather than CT as Na+ laponite showed a similar enhancement in conductivity. Mechanically robust colloidal clay films were prepared on platinum electrodes. After immersion in solutions containing redox active complexes [Co(bpy)3]3+ and [Cr(bpy)3]3+, the films became electroactive when a potential was applied. Cyclic voltammograms obtained for both complexes were found to be of the diffusion controlled type. For [Co(bpy)3]3+ immobilised on clay coated electrodes, a one-step oxidation and four-step reduction wave was observed corresponding to a one electron stepwise reversible reduction of Co(III), through Co(II), Co(I), Co(O) to Co(I) oxidation state. For [Cr(bpy)3]3+ the electrochemistry was complicated by the presence of additional waves corresponding to the dissociation of [Cr(bpy)3]3+ into the diaquo complex. ESR and diffuse reflectance data supported such a mechanism. 29Si, 27Al and 23Na MAS NMR spectroscopy, supported by powder XRD and FTIR, was used to probe the role of clays as reagents in glass forming reactions. 29Si MAS NMR was found to be a very sensitive technique for identifying the presence and relative abundance of crystalline and non-crystalline phases. In thermal reactions of laponite formation of new mineral phases such as forsterite, akermanite, sillimanite and diopside were detected. The relative abundance of each phase was dependent on thermal history, chemical nature and concentration of the modifier oxide present. In continuing work, the effect of selected oxides on the glass forming reactions of a model feldspar composition was investigated using solid state NMR alone. Addition of network modifying oxides generally produced less negative 29Si chemical shifts and larger linewidths corresponding to a wider distribution of Si-O-Si bond angles and lengths, and a dominant aluminosilicate phase with a less polymerised structure than the starting material. 29Si linewidths and 27Al chemical shifts were respectively correlated with cationic potential and Lewis acidity of the oxide cations. Anomalous Al(4) chemical shifts were thought to be due to precipitation of aluminate phases rather than a breakdown in Lowenstein's aluminium avoidance principle.

Monte Carlo and molecular dynamics simulations of methane in potassium montmorillonite clay hydrates at elevated pressures and temperatures

The structure and dynamics of methane in hydrated potassium montmorillonite clay have been studied under conditions encountered in sedimentary basin and compared to those of hydrated sodium montmorillonite clay using computer simulation techniques. The simulated systems contain two molecular layers of water and followed gradients of 150 barkm-1 and 30 Kkm-1 up to a maximum burial depth of 6 km. Methane particle is coordinated to about 19 oxygen atoms, with 6 of these coming from the clay surface oxygen. Potassium ions tend to move away from the center towards the clay surface, in contrast to the behavior observed with the hydrated sodium form. The clay surface affinity for methane was found to be higher in the hydrated K-form. Methane diffusion in the two-layer hydrated K-montmorillonite increases from 0.39×10-9 m2s-1 at 280 K to 3.27×10-9 m2s-1 at 460 K compared to 0.36×10-9 m2s-1 at 280 K to 4.26×10-9 m2s-1 at 460 K in Na-montmorillonite hydrate. The distributions of the potassium ions were found to vary in the hydrates when compared to those of sodium form. Water molecules were also found to be very mobile in the potassium clay hydrates compared to sodium clay hydrates. © 2004 Elsevier Inc. All All rights reserved.