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

Control of the Structural Charge Distribution and Hydration State upon Intercalation of CO2 into Expansive Clay Interlayers

Numerous experimental investigations indicated that expansive clays such as montmorillonite can intercalate CO₂ preferentially into their interlayers and therefore potentially act as a material for CO₂ separation, capture, and storage. However, an understanding of the energy-structure relationship during the intercalation of CO₂ into clay interlayers remains elusive. Here, we use metadynamics molecular dynamics simulations to elucidate the energy landscape associated with CO₂ intercalation. Our free energy calculations indicate that CO₂ favorably partitions into nanoconfined water in clay interlayers from a gas phase, leading to an increase in the CO₂/H₂O ratio in clay interlayers as compared to that in bulk water. CO₂ molecules prefer to be located at the centers of charge-neutral hydrophobic siloxane rings, whereas interlayer spaces close to structural charges tend to avoid CO₂ intercalation. The structural charge distribution significantly affects the amount of CO₂ intercalated in the interlayers. These results provide a mechanistic understanding of CO₂ intercalation in clays for CO₂ separation, capture, and storage.

Competitive Sorption of CO2 with Gas Mixtures in Nanoporous Shale for Enhanced Gas Recovery from Density Functional Theory

CO₂ competitive sorption with shale gas under various conditions from simple to complex pore characteristics is studied using a molecular density functional theory (DFT) that reduces to perturbed chain-statistical associating fluid theory in the bulk fluid region. The DFT model is first verified by grand canonical Monte Carlo simulation in graphite slit pores for pure and binary component systems at different temperatures, pressures, pore sizes, and bulk gas compositions for methane/ethane with CO₂. Then, the model is utilized in multicomponent systems that include CH₄, C₂H₆, and C3+ components of different compositions. It is shown that the selectivity of CO₂ decreases with increases in temperature, pressure, nanopore size, and average molecular weight of shale gas. Extending the model to more realistic situations, we consider the impact of water present in the pore and consider the effect of permeation of fluid molecules into the kerogen that forms the pore walls. The water-graphite interaction is calibrated with contact angle from molecular simulation data from the literature. The kerogen pore model prediction of gas absolute sorption is compared with experimental and molecular simulation values in the literature. It is shown that the presence of water reduces the CO₂ adsorption but improves the CO₂ selectivity. The dissolution of gases into the kerogen matrix also leads to the increase in CO₂ selectivity. The effect of kerogen type and maturity on the gas sorption amount and CO₂ selectivity is also studied. The associated mechanisms are discussed to provide fundamental understanding for gas recovery by CO₂.

Molecular simulation of methane adsorption in slit-like quartz pores

The GCMC and MD methods were used to investigate the structural properties and adsorption behaviors of methane in quartz nanopores.

Monte Carlo study of the adsorption and aggregation of alkyltrimethylammonium chloride on the montmorillonite-water interface

Organically modified clays exhibit adsorption capacities for cations, anions, and nonpolar organic compounds, which make them valuable for various environmental technical applications. To improve the understanding of the adsorption processes, the molecular-scale characterization of the structures of organic aggregates assembled on the external basal surfaces of clay particles is essential. The focus of this Monte Carlo simulation study was on the effects of the surface coverage and the alkyl chain length n on the structures of alkyltrimethylammonium chloride ((C(n)TMA)Cl) aggregates assembled on the montmorillonite-water interface. We found that the amount of adsorbed C(n)TMA(+) ions is independent of the alkyl chain length and increases with the C(n)TMA(+) surface coverage. The C(n)TMA(+) ions predominantly adsorb as inner-sphere complexes; the fraction of outer-sphere adsorbed ions equals only about 10%. The conformational order of the C(n)TMA(+) alkyl chains substantially decreases with decreasing alkyl chain length. In agreement with previous experiments, the amount of C(n)TMA(+) ions that are aggregated at the mineral surface increases with increasing chain length. The maximum value of 0.66 C(n)TMA(+) adsorption complex per unit cell area of the clay surface considerably exceeds the amount of cations required to compensate the negative charge of the montmorillonite surface. Furthermore, in most of the studied systems, fractions of Na(+) surface cations remain adsorbed on montmorillonite. The resulting interfacial positive charge excess is counterbalanced by coadsorbed chloride ions forming ion pairs with both C(n)TMA(+) and Na(+).

Aggregation of alkyltrimethylammonium ions at the cleaved muscovite mica-water interface: a Monte Carlo study

The precise molecular structure of organically modified mineral surfaces is still not well understood. To establish a relation between experimental observations and underlying molecular structure, we performed Monte Carlo simulations of the aggregation behavior of alkyltrimethylammonium surfactants (C(n)TMA(+)) at the interface between C(n)TMACl solution and cleaved K(+)-muscovite. The structures were examined with regard to the influence of varying alkyl chain length n (n = 8, 12, 16) and surface coverage of C(n)TMA(+) ions. The simulation results indicate that the water film structure at the muscovite surface is considerably influenced by the adsorption of C(n)TMA(+). A fraction of the C(n)TMA(+) ions forms inner-sphere and outer-sphere adsorption complexes with nitrogen-surface distances of 3.3-3.8 and 5.5-8.4 Å, respectively. The simulated monolayer aggregates exhibit thicknesses of 31-35, 22-27, and ∼18 Å for C(16)TMA(+), C(12)TMA(+), and C(8)TMA(+), respectively. C(16)TMA(+) and C(12)TMA(+) ions form bilayer aggregates, which show a strong interdigitation of the two opposing layers composing them. The aggregate thicknesses equal 35-39 and 30-35 Å, respectively, and are in agreement with available experimental data. In contrast, the short-chained C(8)TMA(+) ions do not form bilayer aggregates. In agreement with previous experimental studies, the alkyl chains of the aggregated ions show high conformational order markedly decreasing with decreasing chain length. We suggest that the simulated structures represent C(n)TMA(+) aggregates, which are formed on muscovite during the experimentally observed initial equilibration phase characterized by the presence of inorganic ions within the aggregates.

Effects of solid surfaces on hydrate kinetics and stability

Reservoirs of clathrate hydrates of natural gases (hydrates), found worldwide and containing huge amounts of bound natural gases (mostly methane), represent potentially vast and yet untapped energy resources. Since CO2‐containing hydrates are considerably more stable thermodynamically than methane hydrates, if we find a way to replace the original hydrate-bound hydrocarbons with the CO2, two goals can be accomplished at the same time: safe storage of carbon dioxide in hydrate reservoirs, and in situ release of hydrocarbon gas. We have applied the techniques of magnetic resonance imaging as a tool to visualize the conversion of CH4 hydrate within Bentheim sandstone matrix into the CO2 hydrate. Corresponding model systems have been simulated using the phase field theory approach. Our theoretical studies indicate that the kinetic behaviour of the systems closely resembles that of CO2 transport through an aqueous solution. We have interpreted this to mean that the hydrate and the matrix mineral surfaces are separated by liquid-containing channels. These channels will serve as escape routes for released natural gas, as well as distribution channels for injected CO2.