Modeling oceanic sedimentary methane hydrate growth through molecular dynamics simulation

The crystallization process of methane hydrates in a confined geometry resembling seabed porous silica sedimentary conditions has been studied using molecular dynamics simulations. With this objective in mind, a fully atomistic quartz silica slit pore has been designed, and the temperature stability of a methane hydrate crystalline seed in the presence of water and guest molecule methane has been analyzed. NaCl ion pairs have been added in different concentrations, simulating salinity conditions up to values higher than average oceanic conditions. The structure obtained when the hydrate crystallizes inside the pore is discussed, paying special attention to the presence of ionic doping inside the hydrate and the subsequent induced structural distortion. The shift in the hydrate stability conditions due to the increasing water salinity is discussed and compared with the case of unconfined hydrate, concluding that the influence of the confinement geometry and pore hydrophilicity produces a larger deviation in the confined hydrate phase equilibria.

Upper storage-capacity limit and multiple occupancy phenomena in H2hydroquinone clathrates using Monte Carlo and DFT simulations

The upper hydrogen-storage capacity limit of the β-hydroquinone clathrate has been investigated using hybrid Grand-Canonical Monte Carlo/Molecular Dynamics simulations, for temperatures ranging from 77 K to 300 K. The evolution with pressure of the cage occupancies has been monitored in detail, describing the progressive nature of the uptake process. It is found that the storage capacity of the pure β-HQ + H2 clathrate could reach 0.6 wt% (weight percentage) only for pressures above 1400 bar, at ambient temperature. The enhancement of the storage capacities by the multiple occupancy phenomenom was accordingly shown to be very limited by the need for extreme conditions. Following this observation, an unmodified version of the van der Waals & Platteeuw theory was applied allowing for the prediction of experimentally accessible formation pressures. Density functional theory calculations were addittionnaly performed to comprehensively characterize the hydrogen diffusion process within the clathrate crystalline structure, considering different occupancy scenarios.

Modeling of Structure H Carbon Dioxide Clathrate Hydrates: Guest-Lattice Energies, Crystal Structure, and Pressure Dependencies

We performed first-principles computations to investigate the complex interplay of molecular interaction energies in determining the lattice structure and stability of CO₂@sH clathrate hydrates. Density functional theory computations using periodic boundary conditions were employed to characterize energetics and the key structural properties of the sH clathrate crystal under pressure, such as equilibrium lattice volume and bulk modulus. The performance of exchange-correlation functionals together with recently developed dispersion-corrected schemes was evaluated in describing interactions in both short-range and long-range regions of the potential. Structural relaxations of the fully CO₂-filled and empty sH unit cells yield crystal structure and lattice energies, while their compressibility parameters were derived by including the pressure dependencies. The present quantum chemistry computations suggest anisotropy in the compressibility of the sH clathrate hydrates, with the crystal being less compressible along the a-axis direction than along the c-axis one, in distinction from nearly isotropic sI and sII structures. The detailed results presented here give insight into the complex nature of the underlying guest-host interactions, checking earlier assumptions, providing critical tests, and improving estimates. Such entries may eventually lead to better predictions of thermodynamic properties and formation conditions, with a direct impact on emerging hydrate-based technologies.

Exploring CO2 @sI Clathrate Hydrates as CO2 Storage Agents by Computational Density Functional Approaches

The formation of specific clathrate hydrates and their transformation at given thermodynamic conditions depends on the interactions between the guest molecule/s and the host water lattice. Understanding their structural stability is essential to control structure-property relations involved in different technological applications. Thus, the energetic aspects relative to CO₂ @sI clathrate hydrate are investigated through the computation of the underlying interactions, dominated by hydrogen bonds and van der Waals forces, from first-principles electronic structure approaches. The stability of the CO₂ @sI clathrate is evaluated by combining bottom-up and top-down approaches. Guest-free and CO₂ guest-filled aperiodic cages, up to the gradually CO₂ occupation of the entire sI periodic unit cells were considered. Saturation, cohesive and binding energies for the systems are determined by employing a variety of density functionals and their performance is assessed. The dispersion corrections on the non-covalent interactions are found to be important in the stabilization of the CO₂ @sI energies, with the encapsulation of the CO₂ into guest-free/empty cage/lattice being always an energetically favorable process for most of the functionals studied. The PW86PBE functional with XDM or D3(BJ) dispersion corrections predicts a lattice constant in accord to the experimental values available, and simultaneously provides a reliable description for the guest-host interactions in the periodic CO₂ @sI crystal, as well as the energetics of its progressive single cage occupancy process. It has been found that the preferential orientation of the single CO₂ in the large sI crystal cages has a stabilizing effect on the hydrate, concluding that the CO₂ @sI structure is favored either by considering the individual building block cages or the complete sI unit cell crystal. Such benchmark and methodology cross-check studies benefit new data-driven model research by providing high-quality training information, with new insights that indicate the underlying factors governing their structure-driven stability, and triggering further investigations for controlling the stabilization of these promising long-term CO₂ storage materials.

Computational density-functional approaches on finite-size and guest-lattice effects in CO2@sII clathrate hydrate

We performed first-principles computations to investigate guest-host/host-host effects on the encapsulation of the CO₂ molecule in sII clathrate hydrates from finite-size clusters up to periodic 3D crystal lattice systems. Structural and energetic properties were first computed for the individual and first-neighbors clathrate-like sII cages, where highly accurate ab initio quantum chemical methods are available nowadays, allowing in this way the assessment of the density functional (DFT) theoretical approaches employed. The performance of exchange-correlation functionals together with recently developed dispersion-corrected schemes was evaluated in describing interactions in both short-range and long-range regions of the potential. On this basis, structural relaxations of the CO₂-filled and empty sII unit cells yield lattice and compressibility parameters comparable to experimental and previous theoretical values available for sII hydrates. According to these data, the CO₂ enclathration in the sII clathrate cages is a stabilizing process, either by considering both guest-host and host-host interactions in the complete unit cell or only the guest-water energies for the individual clathrate-like sII cages. CO₂@sII clathrates are predicted to be stable whatever the dispersion correction applied and in the case of single cage occupancy are found to be more stable than the CO₂@sI structures. Our results reveal that DFT approaches could provide a good reasonable description of the underlying interactions, enabling the investigation of formation and transformation processes as a function of temperature and pressure.

Local structure and distortions of mixed methane-carbon dioxide hydrates

A vast source of methane is found in gas hydrate deposits, which form naturally dispersed throughout ocean sediments and arctic permafrost. Methane may be obtained from hydrates by exchange with hydrocarbon byproduct carbon dioxide. It is imperative for the development of safe methane extraction and carbon dioxide sequestration to understand how methane and carbon dioxide co-occupy the same hydrate structure. Pair distribution functions (PDFs) provide atomic-scale structural insight into intermolecular interactions in methane and carbon dioxide hydrates. We present experimental neutron PDFs of methane, carbon dioxide and mixed methane-carbon dioxide hydrates at 10 K analyzed with complementing classical molecular dynamics simulations and Reverse Monte Carlo fitting. Mixed hydrate, which forms during the exchange process, is more locally disordered than methane or carbon dioxide hydrates. The behavior of mixed gas species cannot be interpolated from properties of pure compounds, and PDF measurements provide important understanding of how the guest composition impacts overall order in the hydrate structure.

Structural Stability of the CO2 @sI Hydrate: a Bottom-Up Quantum Chemistry Approach on the Guest-Cage and Inter-Cage Interactions

Through reliable first-principles computations, we have demonstrated the impact of CO₂ molecules enclathration on the stability of sI clathrate hydrates. Given the delicate balance between the interaction energy components (van der Waals, hydrogen bonds) present on such systems, we follow a systematic bottom-up approach starting from the individual 5¹² and 5¹² 6² sI cages, up to all existing combinations of two-adjacent sI crystal cages to evaluate how such clathrate-like models perform on the evaluation of the guest-host and first-neighbors inter-cage effects, respectively. Interaction and binding energies of the CO₂ occupation of the sI cages were computed using DF-MP2 and different DFT/DFT-D electronic structure methodologies. The performance of selected DFT functionals, together with various semi-classical dispersion corrections schemes, were validated by comparison with reference ab initio DF-MP2 data, as well as experimental data from x-ray and neutron diffraction studies available. Our investigation confirms that the inclusion of the CO₂ in the cage/s is an energetically favorable process, with the CO₂ molecule preferring to occupy the large 5¹² 6² sI cages compared to the 5¹² ones. Further, the present results conclude on the rigidity of the water cages arrangements, showing the importance of the inter-cage couplings in the cluster models under study. In particular, the guest-cage interaction is the key factor for the preferential orientation of the captured CO₂ molecules in the sI cages, while the inter-cage interactions seems to cause minor distortions with the CO₂ guest neighbors interactions do not extending beyond the large 5¹² 6² sI cages. Such findings on these clathrate-like model systems are in accord with experimental observations, drawing a direct relevance to the structural stability of CO₂ @sI clathrates.

Cage fusion from bi-cages to tri-cages during nucleation of methane hydrate: a DFT-D simulation

Water-cage clusters encapsulating guest molecules are the basic components of hydrate crystal structures. Herein, we investigated the fusion process from bi-cages to tri-cages to probe the nucleation mechanism at the initial stage of CH4 hydrate formation by employing dispersion-corrected density functional theory. We found that tri-cages possess high stability by sharing three, rather than two, polygonal faces. In addition, any mono-cage combined with a nonstandard 4151062 cage could achieve considerable stability regardless of which face is shared; this finding illustrates that 4151062 cages are more likely to appear at the early stages of CH4 hydrate nucleation than other nonstandard cages. We then simulated the Raman spectra of CH4 molecules in water-cage to characterize the spectral characteristics of the CH4 hydrate. The C-H symmetric stretching frequency of encapsulated CH4 molecules red-shifted with increasing mono-cage size, which follows the prediction of the "loose cage-tight cage" model. The symmetric stretching vibrational frequencies of trapped CH4 molecules in the tri-cage revealed a clear red-shift compared with those in the component mono- and bi-cages. The cage fusion process and spectroscopic properties described in this work are expected to provide new atomistic insights into CH4 hydrates at the initial nucleation stage.

A Systematic Protocol for Benchmarking Guest-Host Interactions by First-Principles Computations: Capturing CO2 in Clathrate Hydrates

Clathrate hydrates of CO₂ have been proposed as potential molecular materials in tackling important environmental problems related to greenhouse gases capture and storage. Despite the increasing interest in such hydrates and their technological applications, a molecular-level understanding of their formation and properties is still far from complete. Modeling interactions is a challenging and computationally demanding task, essential to reliably determine molecular properties. First-principles calculations for the CO₂ guest in all sI, sII, and sH clathrate cages were performed, and the nature of the guest-host interactions, dominated by both hydrogen-bond and van der Waals forces, was systematically investigated. Different families of density functionals, as well as pairwise CO₂ @H₂ O model potentials versus wavefunction-based quantum approaches were studied for CO₂ clathrate-like systems. Benchmark energies for new distance-dependent datasets, consisting of potential energy curves sampling representative configurations of the systems at the repulsive, near-equilibrium, and asymptotic/long-range regions of the full-dimensional surface, were generated, and a general protocol was proposed to assess the accuracy of such conventional and modern approaches at minimum and non-minimum orientations. Our results show that dispersion interactions are important in the guest-host stabilization energies of such clathrate cages, and the encapsulation of the CO₂ into guest-free clathrate cages is always energetically favorable. In addition, the orientation of CO₂ inside each cage was explored, and the ability of current promising approaches to accurately describe non-covalent CO₂ @H₂ O guest-host interactions in sI, sII, and sH clathrates was discussed, providing information for their applicability to future multiscale computer simulations.

Diffusion of gas mixtures in the sI hydrate structure

Replacing methane with carbon dioxide in gas hydrates has been suggested as a way of harvesting methane, while at the same time storing carbon dioxide. Experimental evidence suggests that this process is facilitated if gas mixtures are used instead of pure carbon dioxide. We studied the free energy barriers for diffusion of methane, carbon dioxide, nitrogen, and hydrogen in the sI hydrate structure using molecular simulation techniques. Cage hops between neighboring cages were considered with and without a water vacancy and with a potential inclusion of an additional gas molecule in either the initial or final cage. Our results give little evidence for enhanced methane and carbon dioxide diffusion if nitrogen is present as well. However, the inclusion of hydrogen seems to have a substantial effect as it diffuses rapidly and can easily enter occupied cages, which reduces the barriers of diffusion for the gas molecules that co-occupy a cage with hydrogen.

Nitrogen doped nanohoops as promising CO2 capturing devices

The impact of climate change in the face of steady or increasing emissions has made the capture and storage of CO2 a priority issue. Supramolecular chemistry is one of the tools that can be used for this task, due to the possibility of tuning intermolecular interactions for the capture of this gas in a selective and efficient way. In this context, this work presents a novel approach for the capture of CO2 based on n-cycloparaphenylenes ([n]-CPPs) doped with nitrogen atoms. This is the first time that the potential of these structures for the capture of polluting gases has been evaluated. Among all the structures analysed, the one yielding the best results (complexation energy of -32.80 kJ mol-1) contains 4 nitrogen atoms per monomer. The topology of the electron density of the host-guest complex and the nature of its non-covalent interactions have been analyzed in this work in order to explain this high binding energy and identify potential structural modifications to improve it. The capability of this system to be used as a sensing device for CO2 using vibrational spectroscopy is also explored.

CO2 Complexes with Five-Membered Heterocycles: Structure, Topology, and Spectroscopic Characterization

In a first step toward the rational design of macrocyclic structures optimized for CO₂ capture, we systematically explored the potential of 30 five-membered aromatic heterocycles to establish coordinating complexes with this pollutant. The interactions between the two moieties were studied in several orientations, and the obtained complexes were analyzed in terms of electron density and vibrational fingerprint. The former is an aid to provide an in-depth knowledge of the interaction, whereas the latter should help to select structural motifs that have not only good complexation properties but also diagnostic spectroscopic signals.

Computational study of the interplay between intermolecular interactions and CO2 orientations in type I hydrates

Carbon dioxide (CO₂) molecules show a rich orientation landscape when they are enclathrated in type I hydrates. Previous studies have described experimentally their preferential orientations, and some theoretical works have explained, but only partially, these experimental results. In the present paper, we use classical molecular dynamics and electronic density functional theory to advance in the theoretical description of CO₂ orientations within type I hydrates. Our results are fully compatible with those previously reported, both theoretical and experimental, the geometric shape of the cavities in hydrate being, and therefore, the steric constraints, responsible for some (but not all) preferential angles. In addition, our calculations also show that guest-guest interactions in neighbouring cages are a key factor to explain the remaining experimental angles. Besides the implication concerning equation of state hydrate modeling approximations, the conclusion is that these guest-guest interactions should not be neglected, contrary to the usual practice.

Direct transition mechanism for molecular diffusion in gas hydrates

In this work, we use dual cage explicit atomic systems to demonstrate theoretically that direct transitions are feasible through hexagonal and pentagonal faces in type I hydrate without compromising the overall structure integrity.