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.     

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.     

Abnormal Thermal Expansion of Clathrate Hydrates Induced by Asymmetric Guest Molecules

We investigated for the first time the abnormal thermal expansion induced by an asymmetric guest structure using high‐resolution neutron powder diffraction. Three dihydrogen molecules (H₂, D₂, and HD) were tested to explore the guest dynamics and thermal behavior of hydrogen‐doped clathrate hydrates. We confirmed the restricted spatial distribution and doughnut‐like motion of the HD guest in the center of anisotropic sII‐S (sII‐S=small cages of structure II hydrates). However, we failed to observe a mass‐dependent relationship when comparing D₂ with HD. The use of asymmetric guest molecules can significantly contribute to tuning the cage dimension and thus can improve the stable inclusion of small gaseous molecules in confined cages.     

Thermal expansion of structure-H clathrate hydrates

The temperature dependence of the unit cell parameters of two newly identified hexagonal structure clathrate hydrates of hexamethylethane (HME) and 2,2-dimethylbutane (DMB) have been measured by X-ray powder diffraction. The thermal expansion of the two distinct crystallographic axes was found to be inequivalent. However, the coefficients of cubic expansion are comparable to that in the cubic structure I and II hydrates. The larger thermal expansivity in the clathrate hydrates relative to ice is attributed to the weakening of the host lattice due to the internal pressure generated by the rattling motions of the encaged guests.Dedicated to Dr D. W. Davidson in honor of his great contributions to the sciences of inclusion phenomena.     

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.     

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.     

Controlling nonclassical content of clathrate hydrates through the choice of molecular guests and temperature

Low-temperature, low-pressure studies of clathrate hydrates (CHs) have revealed that small ether and other proton-acceptor guests greatly enhance rates of clathrate hydrate nucleation and growth; rapid formation and transformations are enabled at temperatures as low as 110 K, and cool moist vapors containing small ether molecules convert to mixed-gas CHs on a subsecond time scale. More recently, FTIR spectroscopic studies of the tetrahydrofuran (THF)-HCN double clathrate hydrate revealed a sizable frequency shift accompanied by a four-fold intensification of the C-N stretch-mode absorption of the small cage HCN, behavior that is enhanced by cooling and which correlates precisely with similar significant changes of the ether C-O/C-C stretch modes. These temperature-dependent correlated changes in the infrared spectra have been attributed to equilibrated extensive hydrogen bonding of neighboring large- and small-cage guest molecules with water molecules of the intervening wall. An ether guest functions as a proton acceptor, particularly so when complemented by the action of a proton-donor (HCN)/electron-acceptor (SO(2)) small-cage guest. Because guest molecules of the classic clathrate hydrates do not participate in hydrogen bonds with the host water, this H-bonding of guests has been labeled "nonclassical". The present study has been enriched by comparing observed FTIR spectra with high-level molecular orbital computational results for guests and hydrogen-bonded guest-water dimers. Vibrational frequency shifts, from heterodimerization of ethers and water, correlate well with the corresponding observed classical to nonclassical shifts. The new spectroscopic data reveal that the nonclassical structures can contribute at observable levels to CH infrared spectra for a remarkable range of temperatures and choice of guest molecules. By the choice of guest molecules, it is now possible to select the abundance levels of nonclassical configurations, ranging from ～0 to 100%, for a given temperature. This ability is expected to hasten understanding of the role of guest-induced nonclassical structures in the acceleration or inhibition of the rates of CH formation and transformation.     

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.     

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.     

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.     

A new method for screening potential sII and sH hydrogen clathrate hydrate promoters with model potentials

A new predictive computational method for classifying clathrate hydrate promoter molecules is presented, based on the interaction energies between potential promoters and the water networks of sII and sH clathrates. The motivation for this work is identifying promoters for storing hydrogen compactly in clathrate hydrates. As a first step towards achieving this goal, we have developed a general method aimed at distinguishing between molecules that form sII clathrate hydrates and molecules that can-together with a weakly interacting help gas-form sH clathrate hydrates. The new computational method calculates differences in estimated formation energies of the sII and the sH clathrate hydrate. Model interaction potentials have been used, including the electrostatic interactions with newly calculated partial charges for all the considered potential promoter molecules. The methodology can discriminate between the clathrate structure types (sII or sH) formed by each potential promoter with good selectivity, i.e., better than achieved with a simple van der Waals diameter criterion.     

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.     

Selective preparation of beta-cyclodextrin clathrates by solid-phase exchange of included tetrahydrofurane for volatile guests in absence of water

Solid-phase guest-exchange products, prepared from dehydrated clathrate of beta-cyclodextrin (bCD) with tetrahydrofurane (THF) by its saturation with vapor of second guest, were studied using thermal analysis by thermogravimetry combined with mass-spectrometric detection of evolved vapors. This guest-exchange procedure was found to be effective for inclusion of volatile guests, which otherwise require a difficult optimization of preparation conditions. Besides, a performed solid-phase exchange without liquid/solid-phase contact is a standard, technologically friendly procedure of clathrate preparation, which does not require further drying to provide an end product. An observed exchange of THF in the absence of water is rather selective, with some hydrophobic guests being unable to replace THF in its dried clathrate with bCD. This selectivity together with low toxicity of THF may be an advantage for practical applications of this guest-exchange method.     

Spectroscopic Observation of the Hydroxy Position in Butanol Hydrates and Its Effect on Hydrate Stability

In this study, we investigate the crystal structures and phase equilibria of butanols+CH₄+H₂O systems to reveal the hydroxy group positioning and its effects on hydrate stability. Four clathrate hydrates formed by structural butanol isomers are identified with powder X‐ray diffraction (PXRD). In addition, Raman spectroscopy is used to analyze the guest distributions and inclusion behaviors of large alcohol molecules in these hydrate systems. The existence of a free OH indicates that guest molecules can be captured in the large cages of structure II hydrates without any hydrogen‐bonding interactions between the hydroxy group of the guests and the water‐host framework. However, Raman spectra of the binary (1‐butanol+CH₄) hydrate do not show the free OH signal, indicating that there could be possible hydrogen‐bonding interactions between the guests and hosts. We also measure the four‐phase equilibrium conditions of the butanols+CH₄+H₂O systems.     

Effect of small cage guests on hydrogen bonding of tetrahydrofuran in binary structure II clathrate hydrates

Molecular dynamics simulations of the pure structure II tetrahydrofuran clathrate hydrate and binary structure II tetrahydrofuran clathrate hydrate with CO(2), CH(4), H(2)S, and Xe small cage guests are performed to study the effect of the shape, size, and intermolecular forces of the small cages guests on the structure and dynamics of the hydrate. The simulations show that the number and nature of the guest in the small cage affects the probability of hydrogen bonding of the tetrahydrofuran guest with the large cage water molecules. The effect on hydrogen bonding of tetrahydrofuran occurs despite the fact that the guests in the small cage do not themselves form hydrogen bonds with water. These results indicate that nearest neighbour guest-guest interactions (mediated through the water lattice framework) can affect the clathrate structure and stability. The implications of these subtle small guest effects on clathrate hydrate stability are discussed.     

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.     

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.     

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.     

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.     

Molecular dynamics simulation of NMR powder lineshapes of linear guests in structure I clathrate hydrates

We perform molecular dynamics simulations (up to 6 ns) for the structure I clathrate hydrates of linear molecules CS, CS(2), OCS, and C(2)H(2) in large cages at different temperatures in the stability range to determine the angular distribution and dynamics of the guests in the large cages. The long axes of linear guest molecules in the oblate large structure I clathrate hydrate cages are primarily confined near the equatorial plane of the cage rather than axial regions. This non-uniform spatial distribution leads to well-known anisotropic lineshapes in the solid-state NMR spectra of the guest species. We use the dynamic distribution of guest orientations in the cages during the MD simulations at different temperatures to predict the (13)C NMR powder lineshapes of the guests in the large cages. The length of the guests and intermolecular interactions of the guests in the water cages determine the angular distribution and the mobility of the guests in the sI large cages at different temperatures. At low temperatures the range of motion of the guests in the cages are limited and this is reflected in the skew of the predicted (13)C lineshapes. As the guest molecules reach the fast motion limit at higher temperatures, the lineshapes for CS, OCS, and C(2)H(2) are predicted to have the "standard" powder lineshapes of guest molecules.