Stability of rare gas structure H clathrate hydrates

Analysing the stability of He-filled hydrates: how many He atoms fit in the sII crystal?

Clathrate hydrates have the ability to encapsulate atoms and molecules within their cavities, and thus they could be potentially large storage capacity materials. The present work studies the multiple cage occupancy effects in the recently discovered He@sII crystal. On the basis of previous theoretical and experimental findings, the stability of He(1/1)@sII, He(1/4)@sII and He(2/4)@sII crystals was analysed in terms of structural, mechanical and energetic properties. For this purpose, first-principles DFT/DFT-D computations were performed by using both semi-local and non-local functionals, not only to elucidate which configuration is the most energetically favoured, but also to scrutinize the relevance of the long-range dispersion interactions in these kinds of compounds. We have encountered that dispersion interactions play a fundamental role in describing the underlying interactions, and different tendencies were observed depending on the choice of the functional. We found that PW86PBE-XDM shows the best performance, while the non-local functionals tested here were not able to correctly account for them. The present results reveal that the most stable configuration is the one presenting singly occupied D cages and tetrahedrally occupied H cages (He(1/4)@sII) in line with the experimental observation.

Structure and stability of the sH binary hydrate cavity and host-guest versus guest-guest interactions therein: A DFT approach

The intrinsic ability of clathrate hydrates to encage gaseous molecules is explored. Encapsulation ability depends on the cavity size and the type of guest gaseous species in addition to the physical parameters, temperature and pressure. Here we have reported the structure, stability and nature of interaction in dissimilar guest occupied sH hydrate cavity. Diatomic gas molecules and small polyatomic hydrocarbons are considered as guests. The irregular icosahedron (5¹² 6⁸ ) cavity of sH hydrate is the host. Different thermodynamic parameters for the guest molecules encapsulation were calculated using three different hybrid DFT functionals, B3LYP, M05-2X, M06, and moreover using dispersion correction (PBE0-D3). With the consideration of large H-bonded systems the 6-31G* and cc-pVTZ basis sets were used for two set of computations. To disclose the nature of interaction between the host-guest systems as well as the interaction between the guest molecules inside the host the non-covalent interaction (NCI) indices and energy decomposition analysis (EDA) were done. Impact of host-guest and guest-guest interactions are discussed.

Phase equilibria molecular simulations of hydrogen hydrates via the direct phase coexistence approach

We report the three-phase (hydrate-liquid water-vapor) equilibrium conditions of the hydrogen-water binary system calculated with molecular dynamics simulations via the direct phase coexistence approach. A significant improvement of ∼10.5 K is obtained in the current study, over earlier simulation attempts, by using a combination of modifications related to the hydrogen model that include (i) hydrogen Lennard-Jones parameters that are a function of temperature and (ii) the water-guest energy interaction parameters optimized further by using the Lorentz-Berthelot combining rules, based on an improved description of the solubility of hydrogen in water.

From Infrared Spectra to Macroscopic Mechanical Properties of sH Gas Hydrates through Atomistic Calculations

The vibrational characteristics of gas hydrates are key identifying molecular features of their structure and chemical composition. Density functional theory (DFT)-based IR spectra are one of the efficient tools that can be used to distinguish the vibrational signatures of gas hydrates. In this work, ab initio DFT-based IR technique is applied to analyze the vibrational and mechanical features of structure-H (sH) gas hydrate. IR spectra of different sH hydrates are obtained at 0 K at equilibrium and under applied pressure. Information about the main vibrational modes of sH hydrates and the factors that affect them such as guest type and pressure are revealed. The obtained IR spectra of sH gas hydrates agree with experimental/computational literature values. Hydrogen bond’s vibrational frequencies are used to determine the hydrate’s Young’s modulus which confirms the role of these bonds in defining sH hydrate’s elasticity. Vibrational frequencies depend on pressure and hydrate’s O···O interatomic distance. OH vibrational frequency shifts are related to the OH covalent bond length and present an indication of sH hydrate’s hydrogen bond strength. This work presents a new route to determine mechanical properties for sH hydrate based on IR spectra and contributes to the relatively small database of gas hydrates’ physical and vibrational properties.

Gas hydrates in sustainable chemistry

This review includes the current state of the art understanding and advances in technical developments about various fields of gas hydrates, which are combined with expert perspectives and analyses.

A Chiral Gas-Hydrate Structure Common to the Carbon Dioxide-Water and Hydrogen-Water Systems

We present full in situ structural solutions of carbon dioxide hydrate-II and hydrogen hydrate C₀ at elevated pressures using neutron and X-ray diffraction. We find both hydrates adopt a common water network structure. The structure exhibits several features not previously found in hydrates; most notably it is chiral and has large open spiral channels along which the guest molecules are free to move. It has a network that is unrelated to any experimentally known ice, silica, or zeolite network but is instead related to two Zintl compounds. Both hydrates are found to be stable in electronic structure calculations, with hydration ratios in very good agreement with experiment.

Computational investigations of the thermodynamic properties of size-selected water and Ar-water clusters: high-pressure transitions

Classical parallel-tempering Monte Carlo simulations in the isothermal-isobaric ensemble were carried out for the (H2O)20 and Ar(H2O)20 clusters, over a wide range of temperatures (30-1000 K) and pressures (3 kPa-10 GPa) in order to study their thermodynamic properties and structural changes. The TIP4P/ice water model is employed for the water-water interactions, while both semiempirical and ab initio-based potentials are used to model the interaction between the rare-gas atoms and the water molecules. Temperature-pressure phase diagrams for these cluster systems were constructed by employing a two-dimensional multiple-histogram method. Structural changes were detected by analyzing the heat capacity landscape and the Pearson correlation coefficient profile for the interaction energy and volume. Those at high pressure correspond to solid-to-solid transitions and are found to be related to clathrate-like cages around the Ar atom. It is also shown that the formation and thermodynamic stability of such structures are determined by the intermolecular interaction between the rare-gas atoms and the host water molecules.

The structure of CO₂ hydrate between 0.7 and 1.0 GPa

A deuterated sample of CO2 structure I (sI) clathrate hydrate (CO2·8.3 D2O) has been formed and neutron diffraction experiments up to 1.0 GPa at 240 K were performed. The sI CO2 hydrate transformed at 0.7 GPa into the high pressure phase that had been observed previously by Hirai et al. [J. Phys. Chem. 133, 124511 (2010)] and Bollengier et al. [Geochim. Cosmochim. Acta 119, 322 (2013)], but which had not been structurally identified. The current neutron diffraction data were successfully fitted to a filled ice structure with CO2 molecules filling the water channels. This CO2+water system has also been investigated using classical molecular dynamics and density functional ab initio methods to provide additional characterization of the high pressure structure. Both models indicate the water network adapts a MH-III "like" filled ice structure with considerable disorder of the orientations of the CO2 molecule. Furthermore, the disorder appears to be a direct result of the level of proton disorder in the water network. In contrast to the conclusions of Bollengier et al., our neutron diffraction data show that the filled ice phase can be recovered to ambient pressure (0.1 MPa) at 96 K, and recrystallization to sI hydrate occurs upon subsequent heating to 150 K, possibly by first forming low density amorphous ice. Unlike other clathrate hydrate systems, which transform from the sI or sII structure to the hexagonal structure (sH) then to the filled ice structure, CO2 hydrate transforms directly from the sI form to the filled ice structure.

Inter-cage dynamics in structure I, II, and H fluoromethane hydrates as studied by NMR and molecular dynamics simulations

Prospective industrial applications of clathrate hydrates as materials for gas separation require further knowledge of cavity distortion, cavity selectivity, and defects induction by guest-host interactions. The results presented in this contribution show that under certain temperature conditions the guest combination of CH3F and a large polar molecule induces defects on the clathrate hydrate framework that allow intercage guest dynamics. (13)C NMR chemical shifts of a CH3F/CH4/TBME sH hydrate and a temperature analysis of the (2)H NMR powder lineshapes of a CD3F/THF sII and CD3F/TBME sH hydrate, displayed evidence that the populations of CH4 and CH3F in the D and D' cages were in a state of rapid exchange. A hydrogen bonding analysis using molecular dynamics simulations on the TBME/CH3F and TBME/CH4 sH hydrates showed that the presence of CH3F enhances the hydrogen bonding probability of the TBME molecule with the water molecules of the cavity. Similar results were obtained for THF/CH3F and THF/CH4 sII hydrates. The enhanced hydrogen bond formation leads to the formation of defects in the water hydrogen bonding lattice and this can enhance the migration of CH3F molecules between adjacent small cages.

Clathrate hydrates with hydrogen-bonding guests

Clathrate hydrates (CHs) are inclusion compounds in which "tetrahedrally" bonded H(2)O forms a crystalline host lattice composed of a periodic array of cages. The structure is stabilized by guest particles which occupy the cages and interact with cage walls via van der Waals interactions. A host of atoms or small molecules can act as guests; here the focus is on guests that are capable of strong to intermediate H-bonding to water (small ethers, H(2)S, etc.) but nevertheless "choose" this hydrate crystal form in which H-bonding is absent from the equilibrium crystal structure. These CHs can form by exposure of ice to guest molecules at temperatures as low as 100-150 K, at the (low) guest saturation pressure. This is in contrast to the "normal" CHs whose formation typically requires temperatures well above 200 K and at least moderate pressures. The experimental part of this study addresses formation kinetics of CHs with H-bonding guests, as well as transformation kinetics between different CH forms, studied by CH infrared spectroscopy. The accompanying computational study suggests that the unique properties of this family of CHs are due to exceptional richness of the host lattice in point defects, caused by defect stabilization by H-bonding of water to the guests.

Synthesis and characterization of a new structure of gas hydrate

Atoms and molecules <0.9 nm in diameter can be incorporated in the cages formed by hydrogen-bonded water molecules making up the crystalline solid clathrate hydrates. For these materials crystallographic structures generally fall into 3 categories, which are 2 cubic forms and a hexagonal form. A unique clathrate hydrate structure, previously known only hypothetically, has been synthesized at high pressure and recovered at 77 K and ambient pressure in these experiments. These samples contain Xe as a guest atom and the details of this previously unobserved structure are described here, most notably the host-guest ratio is similar to the cubic Xe clathrate starting material. After pressure quench recovery to 1 atmosphere the structure shows considerable metastability with increasing temperature (T <160 K) before reverting back to the cubic form. This evidence of structural complexity in compositionally similar clathrate compounds indicates that the reaction path may be an important determinant of the structure, and impacts upon the structures that might be encountered in nature.

Thermodynamic stability of type-I and type-II clathrate hydrates depending on the chemical species of the guest substances

The free energy differences are calculated for various type-I and type-II clathrate hydrates based on molecular-dynamics simulations, thereby evaluating the thermodynamic stability of the hydrates depending on the chemical species of the guest substances. The simulation systems consist of 27 unit cells, that is, 1242 water molecules and 216 guest molecules for type-I hydrates, and 3672 water molecules and 648 guest molecules for type-II hydrates. The water molecules are described by TIP4P potential, while the guest molecules are described by one-site Lennard-Jones potential, U=4epsilon{(sigma/r)12-(sigma/r)6}, where U is the potential energy, r is the particle distance, sigma is the particle diameter, and epsilon is the energy well depth. The optimal values of sigma that yield the minimum free energy (the best thermodynamic stability) were determined to be 0.39 nm for the type-I hydrates and 0.37 nm for the type-II hydrates.

Molecular dynamics study of the stability of methane structure H clathrate hydrates

Molecular dynamics simulations are used to study the stability of structure H (sH) methane clathrate hydrates in a 3 x 3 x 3 sH unit cell replica. Simulations are performed at experimental conditions of 300 K and 2 GPa for three methane intermolecular potentials. The five small cages of the sH unit cell are assigned methane guest occupancies of one and large cage guest occupancies of one to five are considered. Radial distribution functions, unit cell volumes, and configurational energies are studied as a function of large cage CH(4) occupancy. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies. Large cage occupancy of five is the most stable configuration for a Lennard-Jones united-atom potential and the Tse-Klein-McDonald potential parametrized for condensed methane phases and two for the most stable configuation for the Murad and Gubbins potential.