Clathrate formation in water-noble gas (Hydrogen) systems at high pressures

Phase equilibria in helium-water, neon-water, and hxdrogen-water svstems were studied at pressures up to 15 kbar. The results are compared with the data for the previously investigated water systems with argon, crypton, and xenon. It is concluded that classical polyhedral clathrate hydrates are formed in all the systems, the stability of the hydrates diminishing from xenon to neon. In all the systems, except the xenon system, the hydrates are based on the crystalline framework of ice II. Their formation demands high pressures; the larger the guest molecule, the higher the pressure required. The xenon molecule seems to be too large to fit the cage of the ice II framework; therefore, the xenon hydrate CS-I remains stable up to at least 15 kbar. Translated fromZhurnal Strukturnoi Khimii, Vol. 40, No. 5, pp. 974–980, September–October, 1999.     

Molecular Dynamics Simulation of Structure II Clathrate Hydrates of Xenon and Large Hydrocarbon Guest Molecules

Molecular dynamics (MD) simulations of structure II clathrate hydrates are performed under canonical (NVT) and isobaric–isothermal (NPT) ensembles. The guest molecule as a small help gas is xenon and gases such as cyclopropane, isobutane and propane are used as large hydrocarbon guest molecule (LHGM). The dynamics of structure II clathrate hydrate is considered in two cases: empty small cages and small cages containing xenon. Therefore, the MD results for structure II clathrate hydrates of LHGM and LHGM + Xe are obtained to clarify the effects of guest molecules on host lattice structure. To understand the characteristic configurations of structure II clathrate hydrate the radial distribution functions (RDFs) are calculated for the studied hydrate system. The obtained results indicate the significance of interactions of the guest molecules on stabilizing the hydrate host lattice and these results is consistent with most previous experimental and theoretical investigations.     

Simulation of thermobaric conditions of the formation,composition, and structure of mixed hydrates containing xenon and nitrous oxide

Structural, dynamic, and thermodynamic features of double hydrates of xenon and nitrous oxide are calculated. Thermodynamic stability regions of these hydrates are found. At the atmospheric pressure the xenon hydrate is in the equilibrium with the gas phase at temperatures up to 263 K, whereas at these pressures the nitrous oxide hydrate decomposes already at 218 K. A strong dependence of the equilibrium temperatures and pressures of the formation/decomposition of double nitrous oxide and xenon hydrates on the composition of their mixture in the gas phase is shown.     

Double clathrate hydrates of tetrabutylammonium fluoride + helium, neon, hydrogen and argon at high pressures

Decomposition curves of double ionic clathrate hydrates of tetrabutylammonium fluoride with helium, neon, hydrogen and argon were studied at pressures up to 800 MPa. Formation of double hydrates with helium, neon and hydrogen does not lead to any significant increase of the temperatures of decomposition of these hydrates; at high temperatures the hydrates may decompose even at lower temperatures than the hydrate of pure tetraalkylammonium salt does. Decomposition temperatures of double hydrates with argon in all cases were 4–8 °C higher in comparison with the decomposition temperature of ionic clathrate hydrates of tetrabutylammonium fluoride. We suppose that this behavior is caused by simultaneous effect of three factors on hydrate decomposition temperature: (1) partial filling of the small cavities in the framework of the hydrate with water molecules, (2) weakness of the van der Waals interactions between the gas molecules and the host water molecules, and (3) dissolution of helium, hydrogen and neon in the solution of tetrabutylammonium salt causing a decrease of melting temperatures of the hydrates formed from these solutions.     

A molecular dynamics study of ethanol-water hydrogen bonding in binary structure I clathrate hydrate with CO2

Guest-host hydrogen bonding in clathrate hydrates occurs when in addition to the hydrophilic moiety which causes the molecule to form hydrates under high pressure-low temperature conditions, the guests contain a hydrophilic, hydrogen bonding functional group. In the presence of carbon dioxide, ethanol clathrate hydrate has been synthesized with 10% of large structure I (sI) cages occupied by ethanol. In this work, we use molecular dynamics simulations to study hydrogen bonding structure and dynamics in this binary sI clathrate hydrate in the temperature range of 100-250 K. We observe that ethanol forms long-lived (>500 ps) proton-donating and accepting hydrogen bonds with cage water molecules from both hexagonal and pentagonal faces of the large cages while maintaining the general cage integrity of the sI clathrate hydrate. The presence of the nondipolar CO(2) molecules stabilizes the hydrate phase, despite the strong and prevalent alcohol-water hydrogen bonding. The distortions of the large cages from the ideal form, the radial distribution functions of the guest-host interactions, and the ethanol guest dynamics are characterized in this study. In previous work through dielectric and NMR relaxation time studies, single crystal x-ray diffraction, and molecular dynamics simulations we have observed guest-water hydrogen bonding in structure II and structure H clathrate hydrates. The present work extends the observation of hydrogen bonding to structure I hydrates.     

Effect of Guest–Host Hydrogen Bonding on the Structures and Properties of Clathrate Hydrates

To provide improved understanding of guest–host interactions in clathrate hydrates, we present some correlations between guest chemical structures and observations on the corresponding hydrate properties. From these correlations it is clear that directional interactions such as hydrogen bonding between guest and host are likely, although these have been ignored to greater or lesser degrees because there has been no direct structural evidence for such interactions. For the first time, single‐crystal X‐ray crystallography has been used to detect guest–host hydrogen bonding in structure II (sII) and structure H (sH) clathrate hydrates. The clathrates studied are the tert‐butylamine (tBA) sII clathrate with H₂S/Xe help gases and the pinacolone + H₂S binary sH clathrate. X‐ray structural analysis shows that the tBA nitrogen atom lies at a distance of 2.64 Å from the closest clathrate hydrate water oxygen atom, whereas the pinacolone oxygen atom is determined to lie at a distance of 2.96 Å from the closest water oxygen atom. These distances are compatible with guest–water hydrogen bonding. Results of molecular dynamics simulations on these systems are consistent with the X‐ray crystallographic observations. The tBA guest shows long‐lived guest–host hydrogen bonding with the nitrogen atom tethered to a water HO group that rotates towards the cage center to face the guest nitrogen atom. Pinacolone forms thermally activated guest–host hydrogen bonds with the lattice water molecules; these have been studied for temperatures in the range of 100–250 K. Guest–host hydrogen bonding leads to the formation of Bjerrum L‐defects in the clathrate water lattice between two adjacent water molecules, and these are implicated in the stabilities of the hydrate lattices, the water dynamics, and the dielectric properties. The reported stable hydrogen‐bonded guest–host structures also tend to blur the longstanding distinction between true clathrates and semiclathrates.     

Stability of rare gas structure H clathrate hydrates

Molecular dynamics simulations are used to study the stability of structure H (sH) clathrate hydrates with the rare gases Ne, Ar, Kr, and Xe. Simulations on a 3 x 3 x 3 sH unit cell replica are performed at ambient pressure at 40 and 100 K temperatures. The small and medium (s+m) cages of the sH unit cell are assigned rare gas guest occupancies of 1 and for large (l) cages guest occupancies of 1-6 are considered. Radial distribution functions for guest pairs with occupancies in the l-l, l-(s+m), and (s+m)-(s+m) cages are presented. The unit cell volumes and configurational energies are studied as a function of large cage occupancy for the rare gases. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies at 100 K and 1 bar and 20 kbar pressures. These studies show that the most stable argon clathrate has five guests in the large cages. For krypton and xenon the most stable configurations have three and two guests in the large cages, respectively.     

Using structural data for estimating the stability of water networks in clathrate and semiclathrate hydrates

Methods are discussed for estimating the energy of formation of a clathrate network from ice Ih based on the structural data for clathrate and semiclathrate hydrates. The energies of real networks in clathrate hydrates and of ideal undistorted networks (8 structures) are calculated from the crystal structure data using the Zimmerman-Pimentel quasiharmonic potential. For semiclathrate hydrates, two versions of estimation are suggested-Estimations for 7 hydrates of amines and quaternary ammonium base salts showed that formation of semiclathrate hydrates requires much energy for network deformation. Two methods are described for estimating the network energies according to topological information. The first technique involves calculating the number of cycles with different numbers of sides per unit cell; this method is applicable to networks involving strained cycles with 4, 7, and 8 sides. The second method represents the lattice energy via the partial energies of polyhedron units forming the lattice; it is applicable to the most widespread class of networks constructed from Allen’s polyhedra. The energies of the tetragonal and orthorhombic networks belonging to this class are estimated. A relative stability series involving 9 structures of empty undistorted networks is formed. Translated fromZhurnal Struktumoi Khimii, Vol. 40, No. 2, pp. 287–295, March–April, 1999.     

Phase equilibrium measurements of structure II clathrate hydrates of hydrogen with various promoters

Phase equilibrium measurements of single and mixed organic clathrate hydrates with hydrogen were determined within a pressure range of 2.0-14.0 MPa. The organic compounds studied were furan, 2,5-dihydrofuran, tetrahydropyran, 1,3-dioxolane and cyclopentane. These organic compounds are known to form structure II clathrate hydrates with water. It was found that the addition of hydrogen to form a mixed clathrate hydrate increases the stability compared to the single organic clathrate hydrates. Moreover, the mixed clathrate hydrate also has a much higher stability compared to a pure hydrogen structure II clathrate hydrate. Therefore, the organic compounds act as promoter materials. The stabilities of the single and mixed organic clathrate hydrates with hydrogen showed the following trend in increasing order: 1,3-dioxolane < 2,5-dihydrofuran < tetrahydropyran < furan < cyclopentane, indicating that both size and geometry of the organic compound determine the stability of the clathrate hydrates.     

On the contributions of NMR spectroscopy to clathrate science

This paper presents a short overview of the contribution that NMR spectroscopy has made to clathrate science. Both have been in the realm of the experimental scientist for about 50 years. Different degrees of development of NMR spectroscopy have led to increasing sophistication in the kind of information available, best exemplified by consideration of studies on clathrate hydrates. Initially most results related to guest and host dynamics, progressing to site and structure specific measurement of chemical shift tensors and isotropic shifts that led to the ability to measure sample compositions. Currently developments are leading to time-resolved information, providing new insights into processes such as hydrate formation, as well as magnetic resonance imaging, leading to space-resolved studies of hydrate formation and decomposition. B NRCC No. 42188. Translated fromZhurnal Strukturnoi Khimii, Vol. 40, No. 5, pp. 809–821, September–October, 1999.     

Experimental and theoretical analysis of the rotational Raman spectrum of hydrogen molecules in clathrate hydrates

The Raman spectra of H(2) and HD molecules in simple hydrogen and binary hydrogen-tetrahydrofuran clathrate hydrates have been measured at temperatures as low as 20 K. The rotational bands of trapped molecules in simple and binary hydrates have been analyzed, and the contributions originating from hydrogen molecules in the large cages have been separated from those in the small cages. A theoretical model, consisting in rigid cages enclosing interacting hydrogen molecules, has been exploited to calculate, on the basis of quantum mechanics, the Raman intensity of the rotational transitions for up to two interacting molecules in one cage. A comparison with experiment leads to a clear interpretation of sidebands appearing in the Raman rotational lines. The quantitative agreement between theory and experiment obtained in some cases clarifies the importance of the choice of the interaction potential, and of the proton disorder in the clathrate crystal.     

How much carbon dioxide can be stored in the structure H clathrate hydrates?: a molecular dynamics study

The stability of structure H (sH) carbon dioxide clathrate hydrates at three temperature-pressure conditions are determined by molecular dynamics simulations on a 3x3x3 sH unit cell replica. Simulations are performed at 100 K at ambient pressure, 273 K at 100 bars and also 300 K and 5.0 kbars. The small and medium cages of the sH unit cell are occupied by a single carbon dioxide guest and large cage guest occupancies of 1-5 are considered. Radial distribution functions are given for guests in the large cages and unit cell volumes and configurational energies are studied as a function of large cage CO(2) occupancy. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies at three temperature/pressure conditions stated above. At the low temperature, large cage occupancy of 5 is the most stable while at the higher temperature, the occupancy of 3 is the most favored. Calculations are also performed to show that the CO(2) sH clathrate is more stable than the methane clathrate analog. Implications on CO(2) sequestration by clathrate formation are discussed.     

The Cubic Superstructure-I of Tetrabutylammonium Fluoride (C4H9)4NF·29.7H2O Clathrate Hydrate

The crystal structure of (C₄H₉)₄NF·29.7H₂O clathrate hydrate (ionic clathrate) determined by X-Ray analysis is reported. The structure is cubic, I $\overline 4 3dThe crystal structure of (C₄H₉)₄NF·29.7H₂O clathrate hydrate (ionic clathrate) determined by X-Ray analysis is reported. The structure is cubic, I , a = 24.375(3) ? (150 K). Its idealized water framework is analogous to that of cubic structure-I of gas hydrates but with eight-fold unit cell, that is a superstructure of cubic structure-I. This is the last structure found in the binary system (C₄H₉)₄NF–H₂O which was not characterized by X-ray analysis earlier. The structure features of the compound under investigation and others existing in H₂O–(C₄H₉)₄NF binary system are discussed.     

Compressibility of Gas Hydrates

Experimental data on the pressure dependence of unit cell parameters for the gas hydrates of ethane (cubic structure I, pressure range 0–2 GPa), xenon (cubic structure I, pressure range 0–1.5 GPa) and the double hydrate of tetrahydrofuran+xenon (cubic structure II, pressure range 0–3 GPa) are presented. Approximation of the data using the cubic Birch–Murnaghan equation, P=1.5B₀[(V₀/V)^7/3?(V₀/V)^5/3], gave the following results: for ethane hydrate V₀=1781 ų, B₀=11.2 GPa; for xenon hydrate V₀=1726 ų, B₀=9.3 GPa; for the double hydrate of tetrahydrofuran+xenon V₀=5323 ų, B₀=8.8 GPa. In the last case, the approximation was performed within the pressure range 0–1.5 GPa; it is impossible to describe the results within a broader pressure range using the cubic Birch–Murnaghan equation. At the maximum pressure of the existence of the double hydrate of tetrahydrofuran+xenon (3.1 GPa), the unit cell volume was 86 % of the unit cell volume at zero pressure. Analysis of the experimental data obtained by us and data available from the literature showed that 1) the bulk modulus of gas hydrates with classical polyhedral structures, in most cases, are close to each other and 2) the bulk modulus is mainly determined by the elasticity of the hydrogen‐bonded water framework. Variable filling of the cavities with guest molecules also has a substantial effect on the bulk modulus. On the basis of the obtained results, we concluded that the bulk modulus of gas hydrates with classical polyhedral structures and existing at pressures up to 1.5 GPa was equal to (9±2) GPa. In cases when data on the equations of state for the hydrates were unavailable, the indicated values may be recommended as the most probable ones.     

Multiple H2 occupancy of cages of clathrate hydrate under mild conditions

Experiments were carried out by reacting H(2) gas with N(2) hydrate at a temperature of 243 K and a pressure of 15 MPa. The characterizations of the reaction products indicated that multiple H(2) molecules can be loaded into both large and small cages of structure II clathrate hydrates. The realization of multiple H(2) occupancy of hydrate cages under moderate conditions not only brings new insights into hydrogen clathrates but also refreshes the perspective of clathrate hydrates as hydrogen storage media.     

The nuclear magnetic resonance line shapes of Xe in the cages of clathrate hydrates

We report, for the first time, a prediction of the line shapes that would be observed in the (129)Xe nuclear magnetic resonance (NMR) spectrum of xenon in the cages of clathrate hydrates. We use the dimer tensor model to represent pairwise contributions to the intermolecular magnetic shielding tensor for Xe at a specific location in a clathrate cage. The individual tensor components from quantum mechanical calculations in clathrate hydrate structure I are represented by contributions from parallel and perpendicular tensor components of Xe-O and Xe-H dimers. Subsequently these dimer tensor components are used to reconstruct the full magnetic shielding tensor for Xe at an arbitrary location in a clathrate cage. The reconstructed tensors are employed in canonical Monte Carlo simulations to find the Xe shielding tensor component along a particular magnetic field direction. The shielding tensor component weighted according to the probability of finding a crystal fragment oriented along this direction in a polycrystalline sample leads to a predicted line shape. Using the same set of Xe-O and Xe-H shielding functions and the same Xe-O and Xe-H potential functions we calculate the Xe NMR spectra of Xe atom in 12 distinct cage types in clathrate hydrates structures I, II, H, and bromine hydrate. Agreement with experimental spectra in terms of the number of unique tensor components and their relative magnitudes is excellent. Agreement with absolute magnitudes of chemical shifts relative to free Xe atom is very good. We predict the Xe line shapes in two cages in which Xe has not yet been observed.     

Clathrate hydrates of hexagonal structure III at high pressures: Structures and phase diagrams

It was shown in this work that the clathrate hydrates of Hexagonal Structure III, formed in the ternary systems 1-methylpiperazine-help gas-water and iso-amyl alcohol-help gas-water are stable in a wide range of pressures. The decomposition curves of these hydrates were studied for the first time up to the pressures 1 GPa. Ar, Kr, Xe and CH₄ were utilized as the help gases. In a number of the systems studied, high pressure phases were revealed that presumably form due to the distortion of the corresponding low pressure hydrate structures.     

Equilibrium conditions of clathrate hydrates formed from carbon dioxide and aqueous acetone solutions

Equilibrium conditions of clathrate hydrates formed from carbon dioxide and aqueous acetone solutions were experimentally measured at temperatures between 269.2 and 281.4 K and pressures up to 3.98 MPa. The acetone concentrations in solutions were investigated from 0.04 to 0.40 mass fractions. The experimental results suggested a transition in hydrate structure from structure I to another structure for acetone solutions between 0.04 and 0.12 mass fractions of acetone. The hydrate structure was suggested to be structure II which was the most stable with a 0.16 mass fraction acetone solution. For more than 0.16 mass fraction of acetone, the equilibrium conditions of the hydrate were shifted to lower temperatures as acetone concentrations increased.     

Gas‐Hydrate Packing and Stability at High Pressures

The decomposition temperatures of double gas hydrates of tetrahydrofuran with noble gases from krypton to helium at pressures up to 15 kbar were found by differential thermal analysis. The stability of hydrates was shown to rise as their packing coefficient increases. Krypton and argon hydrates retain the original cubic structure II in the whole pressure range. In neon and helium systems, polyhedral double hydrates have upper stability limits at 7.4 and 6.0 kbar, respectively.     

Formation and spectra of clathrate hydrates of methanol and methanol-ether mixtures

Infrared spectra of mixed clathrate hydrates, with either ethylene oxide (EO) or tetrahydrofuran (THF) and methanol molecules as the guest species, have been obtained from thin films prepared by vapor deposition of D₂O mixtures in the 115–130 K range. Although methanol acts as a suppressant to the direct vapor deposition of a type I clathrate with EO, nearly complete conversion of 115 K amorphous codeposits, to the crystalline mixed clathrate, occurs upon warming near 150 K. By contrast, the type II clathrate of THF shows an increased crystalline quality when methanol is included in the vapor deposits of the mixed clathrate hydrate at 130 K. The observation of the O---D stretch-mode band of weakly bonded CD₃OD near 2575 cm⁻¹ is part of the evidence that the methanol molecules are encaged. However, as shown theoretically by Tanaka, the clathrate hydrates of methanol, even when mixed with an ether help gas, are not stable structures but form at low temperatures because of kinetic factors, only to decompose in the 140–160 K range. Attempts to prepare a simple type I or type II clathrate hydrate of methanol have produced mixed results. Limited amounts of clathrate hydrate form during deposition but annealing does not result in complete conversion to crystalline clathrates, particularly for host : guest ratios of 17 : 1.