Water cavities of sH clathrate hydrate stabilized by molecular hydrogen: phase equilibrium measurements

In this experimental phase equilibrium study, we show for the first time that it is possible to stabilize structure sH of hydrogen clathrate hydrate with the help of some selected promoters. It was established that the formation pressures of these systems are significantly higher than that of structure sII of hydrogen clathrate hydrate when tetrahydrofuran (THF) is used as a promoter. Although no experimental evidence is available yet, it is estimated that the hydrogen storage capacity of structure sH can be as high as 1.4 wt % of H2, which is about 40% higher compared to the hydrogen storage capacity in structure sII.     

Water cavities of sH clathrate hydrate stabilized by molecular hydrogen

X-ray diffraction and Raman spectroscopic measurements confirm that molecular hydrogen can be contained within the small water cavities of a binary sH clathrate hydrate using large guest molecules that stabilize the large cavity. The potential increase in hydrogen storage could be more than 40% when compared with binary sII hydrates. This work demonstrates the stabilization of hydrogen in a hydrate structure previously unknown for encapsulating molecular hydrogen, indicating the potential for other inclusion compound materials with even greater hydrogen storage capabilities.     

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.     

Structure transition and swapping pattern of clathrate hydrates driven by external guest molecules

We first report here that under strong surrounding gas of external CH4 guest molecules the sII and sH methane hydrates are structurally transformed to the crystalline framework of sI, leading to a favorable change of the lattice dimension of the host-guest networks. The high power decoupling 13C NMR and Raman spectroscopies were used to identify structure transitions of the mixed CH4 + C2H6 hydrates (sII) and hydrocarbons (methylcyclohexane, isopentane) + CH4 hydrates (sH). The present findings might be expected to provide rational evidences regarding the preponderant occurrence of naturally occurring sI methane hydrates in marine sediments. More importantly, we note that the unique and cage-specific swapping pattern of multiguests is expected to provide a new insight for better understanding the inclusion phenomena of clathrate materials.     

Cage occupancies in the high pressure structure H methane hydrate: a neutron diffraction study

A neutron diffraction study was performed on the CD(4) : D(2)O structure H clathrate hydrate to refine its CD(4) fractional cage occupancies. Samples of ice VII and hexagonal (sH) methane hydrate were produced in a Paris-Edinburgh press and in situ neutron diffraction data collected. The data were analyzed with the Rietveld method and yielded average cage occupancies of 3.1 CD(4) molecules in the large 20-hedron (5(12)6(8)) cages of the hydrate unit cell. Each of the pentagonal dodecahedron (5(12)) and 12-hedron (4(3)5(6)6(3)) cages in the sH unit cell are occupied with on average 0.89 and 0.90 CD(4) molecules, respectively. This experiment avoided the co-formation of Ice VI and sH hydrate, this mixture is more difficult to analyze due to the proclivity of ice VI to form highly textured crystals, and overlapping Bragg peaks of the two phases. These results provide essential information for the refinement of intermolecular potential parameters for the water-methane hydrophobic interaction in clathrate hydrates and related dense structures.     

Molecular dynamics study of structure H clathrate hydrates of methane and large guest molecules

Methane storage in structure H (sH) clathrate hydrates is attractive due to the relatively higher stability of sH as compared to structure I methane hydrate. The additional stability is gained without losing a significant amount of gas storage density as happens in the case of structure II (sII) methane clathrate. Our previous work has showed that the selection of a specific large molecule guest substance (LMGS) as the sH hydrate former is critical in obtaining the optimum conditions for crystallization kinetics, hydrate stability, and methane content. In this work, molecular dynamics simulations are employed to provide further insight regarding the dependence of methane occupancy on the type of the LMGS and pressure. Moreover, the preference of methane molecules to occupy the small (5(12)) or medium (4(3)5(6)6(3)) cages and the minimum cage occupancy required to maintain sH clathrate mechanical stability are examined. We found that thermodynamically, methane occupancy depends on pressure but not on the nature of the LMGS. The experimentally observed differences in methane occupancy for different LMGS may be attributed to the differences in crystallization kinetics and/or the nonequilibrium conditions during the formation. It is also predicted that full methane occupancies in both small and medium clathrate cages are preferred at higher pressures but these cages are not fully occupied at lower pressures. It was found that both small and medium cages are equally favored for occupancy by methane guests and at the same methane content, the system suffers a free energy penalty if only one type of cage is occupied. The simulations confirm the instability of the hydrate when the small and medium cages are empty. Hydrate decomposition was observed when less than 40% of the small and medium cages are occupied.     

Structure transition and tuning pattern in the double (tetramethylammonium hydroxide + gaseous guests) clathrate hydrates

In this study, we present an extraordinary structural transition accompanying the occurrence of more than two coexisting clathrate hydrate phases in the double (CH4 + tetramethylammonium hydroxide (Me(4)NOH)) and (H2 + Me(4)NOH) ionic clathrate hydrates using solid-state NMR spectroscopy (high-powered decoupling and CP/MAS) and powder X-ray diffraction. It was confirmed that structure-I (sI) and structure-II (sII) hydrates coexist as the water concentration increases. In the Me(4)NOH-depleted region, the unique tuning phenomenon was first observed at a chemical shift of -8.4 ppm where relatively small gaseous CH4 molecules partly occupy the sII large cages (sII-L), pulling out large cationic Me(4)N+ that is considered to be strongly bound with the surrounding host lattices. Moreover, we note that, while pure Me(4)NOH.16H(2)O clathrate hydrates melted at 249 K under atmospheric pressure conditions, the double (CH4 + Me(4)NOH) clathrate hydrate maintained a solid state up to approximately 283 K under 120 bar of CH4 with a conductivity of 0.065 S cm(-1), suggesting its potential use as a solid electrolyte. The present results indicate that ionic contributions must be taken into account for ionic clathrate hydrate systems because of their distinctive guest dynamic behavior and structural patterns. In particular, microscopic analyses of ionic clathrate hydrates for identifying physicochemical characteristics are expected to provide new insights into inclusion chemistry.     

Thermal expansivity of tetrahydrofuran clathrate hydrate with diatomic guest molecules

The guest dynamics and thermal behavior occurring in the cages of clathrate hydrates appear to be too complex to be clearly understood through various structural and spectroscopic approaches, even for the well-known structures of sI, sII, and sH. Neutron diffraction studies have recently been carried out to clarify the special role of guests in expanding the host water lattices and have contributed to revealing the influence factors on thermal expansivity. Through this letter we attempt to address three noteworthy features occurring in guest inclusion: (1) the effect of guest dimension on host water lattice expansion; (2) the effect of thermal history on host water lattice expansion; and (3) the effect of coherent/incoherent scattering cross sections on guest thermal patterns. The diatomic guests of H 2, D 2, N 2, and O 2 have been selected for study, and their size and mass dependence on the degree of lattice expansion have been examined, and four sII clathrate hydrates with tetrahydrofuran (THF) have been synthesized in order to determine their neutron powder diffraction patterns. After thermal cycling, the THF + H 2 clathrate hydrate is observed to exhibit an irreversible plastic deformation-like pattern, implying that the expanded lattices fail to recover the original state by contraction. The host-water cage dimension after degassing the guest molecules remains as it was expanded, and thus host-guest as well as guest-guest interactions will be altered if guest uptake reoccurs.     

Micro-Raman investigations of mixed gas hydrates

We report laser Raman spectroscopic measurements on mixed hydrates (clathrates), with guest molecules tetrahydrofuran (THF) and methane (CH₄), at ambient pressure and at temperatures from 175 to 280 K. Gas hydrates were synthesized with different concentrations of THF ranging from 5.88 to 1.46 mol%. In all cases THF molecules occupied the large cages of sII hydrate. The present studies demonstrate formation of sII clathrates with CH₄ molecules occupying unfilled cages for concentrations of THF ranging from 5.88 to 2.95 mol%. The Raman spectral signature of hydrates with 1.46 mol% THF are distinctly different; hydrate growth was non-uniform and structural transformation occurred from sII to sI prior to clathrate melting.     

Hydrogen storage in sH hydrates: a Monte Carlo study

Grand canonical Monte Carlo simulations are performed to evaluate the hydrogen-storage capacity of the recently discovered hydrogen hydrates of the sH type, at 274 K and up to 500 MPa. First, the pure H2 hydrate is investigated in order to determine the upper limit of H 2 content in sH hydrates. It is found that the storage capacity of the hypothetical pure H2 hydrate could reach 3.6 wt % at 500 MPa. Depending on pressure, the large cavity of this hydrate can accommodate up to eight H2 molecules, while the small and medium ones are singly occupied even at pressures as high as 500 MPa. Next, the binary H2-methylcyclohexane sH hydrate is examined. In this case, the small and medium cavities are again singly occupied, resulting in a maximum H2 uptake of 1.4 wt %. Finally, the results from simulations on pure H2 and binary hydrates are utilized to investigate the potential of H2 storage in sH hydrates where the promoter molecules occupy the medium instead of the large cavities.     

Spectroscopic signatures of halogens in clathrate hydrate cages. 2. Iodine

UV-vis and Raman spectroscopy were used to study iodine molecules trapped in sII clathrate hydrate structures stabilized by THF, CH(2)Cl(2), or CHCl(3). The spectra show that the environment for iodine inside the water cage is significantly less perturbed than either in aqueous solution or in amorphous water-ice. The resonance Raman progression of I(2) in THF clathrate hydrate can be observed up to v = 6 when excited at 532 nm. The extracted vibrational frequency omega e = 214 +/- 1 cm(-1) is the same as that of the free molecule to within experimental error. At the same time, the UV-vis absorption spectrum of I(2) in the sII hydrate exhibits a relatively large, 1440 cm(-1), blue-shift. This is mainly ascribed to the differential solvation of the I(2) electronic states. We conclude that iodine in sII hydrate resides in a 5(12)6(4) cavity, in which the ground-state I(2) potential is not significantly perturbed by the hydrate lattice. In contrast, in water and in ice, the valence absorption band of I(2) is dramatically broadened and blue-shifted by 3000 cm(-1), and the resonance Raman scattering is effectively quenched. These observations are shown to be consistent with a strong interaction between water molecule and iodine through the lone pair of electrons on water as in the case of bromine in the same media. The results presented here, and the stability of other halogen hydrates, were used to test the predictions of simple models and force-field calculations of the host cage-guest association energy.     

Time-resolved in situ neutron diffraction studies of gas hydrate: transformation of structure II (sII) to structure I (sI).

We report the in situ observation from diffraction data of the conversion of a gas hydrate with the structure II (sII) lattice to one with the structure I (sI) lattice. Initially, the in situ formation, dissociation, and reactivity of argon gas clathrate hydrate was investigated by time-of-flight neutron powder diffraction at temperatures ranging from 230 to 263 K and pressures up to 5000 psi (34.5 MPa). These samples were prepared from deuterated ice crystals and transformed to hydrate by pressurizing the system with argon gas. Complete transformation from D(2)O ice to sII Ar hydrate was observed as the sample temperature was slowly increased through the D(2)O ice melting point. The transformation of sII argon hydrate to sI hydrate was achieved by removing excess Ar gas and exposing the hydrate to liquid CO(2) by pressurizing the Ar hydrate with CO(2). Results suggest the sI hydrate formed from CO(2) exchange in argon sII hydrate is a mixed Ar/CO(2) hydrate. The proposed exchange mechanism is consistent with clathrate hydrate being an equilibrium system in which guest molecules are exchanging between encapsulated molecules in the solid hydrate and free molecules in the surrounding gas or liquid phase.     

Prediction of hydrogen hydrate equilibrium by integrating ab initio calculations with statistical thermodynamics

This paper addresses a new calculation approach for the prediction of hydrogen hydrate equilibrium by introducing the concept of a single hydrogen cluster in one cavity. By integrating ab initio calculations with classical statistical thermodynamics, this approach enables the van der Waals model to predict the dissociation pressure of hydrogen hydrates. Compared to hydrates formed by light hydrocarbon gases, structure II (sII) hydrogen hydrates stably encage two and four hydrogen molecules in the small and large cavities, respectively. By treating two hydrogen molecules or four hydrogen molecules as one rigid body cluster, we determine ab initio binding energies between water molecules and hydrogen clusters at the MP2 level with the 6-31++G(2d,2p) basis set. These binding energies will be used to determine the parameters of the Exp-6 potential function from which the smooth cell potential and the Langmuir constant of each cluster are calculated. Then, the dissociation pressure is determined using the Zele-Lee-Holder cell distortion model: 105, 625, and 2000 bar at 150, 200, and 250 K, respectively.     

13C NMR studies of hydrocarbon guests in synthetic structure H gas hydrates: experiment and computation

(13)C NMR chemical shifts were measured for pure (neat) liquids and synthetic binary hydrate samples (with methane help gas) for 2-methylbutane, 2,2-dimethylbutane, 2,3-dimethylbutane, 2-methylpentane, 3-methylpentane, methylcyclopentane, and methylcyclohexane and ternary structure H (sH) clathrate hydrates of n-pentane and n-hexane with methane and 2,2-dimethylbutane, all of which form sH hydrates. The (13)C chemical shifts of the guest atoms in the hydrate are different from those in the free form, with some carbon atoms shifting specifically upfield. Such changes can be attributed to conformational changes upon fitting the large guest molecules in hydrate cages and/or interactions between the guests and the water molecules of the hydrate cages. In addition, powder X-ray diffraction measurements revealed that for the hexagonal unit cell, the lattice parameter along the a-axis changes with guest hydrate former molecule size and shape (in the range of 0.1 ?) but a much smaller change in the c-axis (in the range of 0.01 ?) is observed. The (13)C NMR chemical shifts for the pure hydrocarbons and all conformers were calculated using the gauge invariant atomic orbital method at the MP2/6-311+G(2d,p) level of theory to quantify the variation of the chemical shifts with the dihedral angles of the guest molecules. Calculated and measured chemical shifts are compared to determine the relative contribution of changes in the conformation and guest-water interactions to the change in chemical shift of the guest upon clathrate hydrate formation. Understanding factors that affect experimental chemical shifts for the enclathrated hydrocarbons will help in assigning spectra for complex hydrates recovered from natural sites.     

Molecular hydrogen occupancy in binary THF-H2 clathrate hydrates by high resolution neutron diffraction

We have determined the time-space average filling of hydrogen molecules in a binary tetrahydrofuran (THF)-d(8) + D(2) sII clathrate hydrate using high resolution neutron diffraction. The filling of hydrogen in the lattice of a THF-d(8) clathrate hydrate occurred upon pressurization. The hydrogen molecules were localized in the small dodecahedral cavities at 20 K, with nuclear density from the hydrogen approximately spherically distributed and centered in the small cavity. With a formation pressure of 70 MPa, molecular hydrogen was found to only singly occupy the sII small cavity. This result helps explain discrepancies about the hydrogen occupancy in the THF binary hydrate system.     

Gas hydrate phase equilibria for propane in the presence of additive components

A proper discussion on the possibility and feasibility of technological applications for gas hydrates requires knowledge of the phase behaviour and its relation to the gas hydrate structure and its occupation. This paper presents experimental data on gas hydrate phase equilibria for the system water+propane and for various systems of the kind water+propane+additive. The additives considered are tetrahydropyran, cyclobutanone and cyclohexane, which are assumed to occupy the large cavity of structure II (sII) hydrate, and methylcyclohexane that is a typical structure H (sH) hydrate former. All additives have in common that they are very poorly soluble in water and, therefore, an additional liquid phase is present in these systems. The pressure for the equilibrium hydrate–liquid water–vapour (H–L_w–V) in the system water+propane is reduced upon addition of each of these components. Simultaneously, the hydrate equilibrium hydrate–liquid water–liquid propane (H–L_w–L_C3H8) is shifted to lower temperatures. These observations can be explained in terms of mutual miscibility of propane and the additive component. However, it cannot be excluded that propane molecules are exchanged by additive molecules in occupying the large cavity of sII.     

Thermal expansivity for sI and sII clathrate hydrates

Knowledge of thermal expansivity can aid in the understanding of both microscopic and macroscopic behavior of clathrate hydrates. Diffraction studies have shown that hydrate volume changes significantly (as much as 1.5% over 50 K) as a function of temperature. It has been demonstrated previously via statistical mechanics that a minor change in hydrate volume (e.g., a 1.5% change in volume or 0.5% change in lattice parameter) can lead to a major change in the predicted hydrate formation pressure (e.g., >15% at >100 MPa for methane). Because of this sensitivity, hydrate thermal expansivity measurements, for both Structures I and II with various guests, are needed help quantify volume distortions in hydrate lattices to ensure accurate hydrate phase equilibria predictions. In addition to macroscopic phase equilibria, the thermal expansion of different hydrates can give information about the interactions between the guest molecules and the host lattice. In this work, the hydrate lattice parameters for four Structure I (C2H6, CO2, 47% C2H6 + 53% CO2, and 85% CH4 + 15% CO2) and seven Structure II (C3H8, 60% CH4 + 40% C3H8, 30% C2H6 + 70% C3H8, 18% CO2 + 82% C3H8, 87.6% CH4 + 12.4% i-C4H10, 95% CH4 + 5% C5H10O, and a natural gas mixture) systems were measured as a function of temperature. The lattice parameter measurements were combined with existing literature values. Both sI and sII hydrates, with a few exceptions, had a common thermal expansivity, independent of hydrate guest. Many guest-dependent correlations for linear thermal expansivity have been proposed. However, we present two guest-independent, structure-dependent correlations for sI and sII lattices, which have been developed to express the normalized hydrate lattice parameters (and therefore volume) as a function of temperature.     

Unusual crystalline and polycrystalline structures in methane hydrates

We present previously unreported crystalline and polycrystalline structures for methane hydrates obtained at relatively low pressures. These structures, which contain unusual cages with 12 pentagonal faces and 3 hexagonal faces, were observed during the atomistic simulations of the crystal growth of sI and sII methane hydrates. These 51263 cages have a significant impact on the structure of the resulting crystal and could explain several experimental observations regarding in-situ transformations between sI and sII hydrates. We document a previously unidentified structure of methane hydrates which we designate structure sK. Additionally, we predict a polycrystalline structure consisting of this new hydrate and sI and suggest a mechanism for the formation of a polycrystalline structure consisting of sequences of sI and sII hydrates.     

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.     

常见客体分子对笼型水合物晶格常数的影响

Natural gas hydrates are considered as ideal alternative energy resources for the future, and the relevant basic and applied research has become more attractive in recent years. The influence of guest molecules on the hydrate crystal lattice parameters is of great significances to the understanding of hydrate structural characteristics, hydrate formation/decomposition mechanisms, and phase stability behaviors. In this study, we test a series of artificial hydrate samples containing different guest molecules (e.g. methane, ethane, propane, iso-butane, carbon dioxide, tetrahydrofuran, methane + 2, 2-dimethylbutane, and methane + methyl cyclohexane) by a low-temperature powder X-ray diffraction (PXRD). Results show that PXRD effectively elucidates structural characteristics of the natural gas hydrate samples, including crystal lattice parameters and structure types. The relationships between guest molecule sizes and crystal lattice parameters reveal that different guest molecules have different controlling behaviors on the hydrate types and crystal lattice constants. First, a positive correlation between the lattice constants and the van der Waals diameters of homologous hydrocarbon gases was observed in the single-guest-component hydrates. Small hydrocarbon homologous gases, such as methane and ethane, tended to form sI hydrates, whereas relatively larger molecules, such as propane and iso-butane, generated sⅡ hydrates. The hydrate crystal lattice constants increased with increasing guest molecule size. The types of hydrates composed of oxygen-containing guest molecules (such as CO₂ and THF) were also controlled by the van der Waals diameters. However, no positive correlation between the lattice constants and the van der Waals diameters of guest molecules in hydrocarbon hydrates was observed for CO₂ hydrate and THF hydrate, probably due to the special interactions between the guest oxygen atoms and hydrate "cages". Furthermore, the influences of the macromolecules and auxiliary small molecules on the lengths of the different crystal axes of the sH hydrates showed inverse trends. Compared to the methane + 2, 2-dimethylbutane hydrate sample, the length of the a-axis direction of the methane + methyl cyclohexane hydrate sample was slightly smaller, whereas the length of the c-axis direction was slightly longer. The crystal a-axis length of the sH hydrate sample formed with nitrogen molecules was slightly longer, whereas the c-axis was shorter than that of the methane + 2, 2-dimethylbutane hydrate sample at the same temperature.