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.     

氮气笼型水合物的激光拉曼光谱

在253K和16MPa的压力下,于实验室内合成了氮气水合物,用显微共焦拉曼光谱对其N-N和O-H键伸缩振动的光谱特征进行了研究．结果表明,氮气水合物中的N-N和O—H键的拉曼峰分别为2322．4和3092．1cm^-1,与天然的空气水合物中的数据十分接近．另外,还测定了液氮和溶解于水中的氮分子中N—N键的拉曼峰值,分别为2326．6和2325．0cm^-1．氮气笼型水合物分解的拉曼谱图表明,氮分子同时进入水合物的大笼和小笼中,但由于氮分子在大、小笼中的环境氛围十分接近,其拉曼位移相差不大,故拉曼谱图只能显示N—N键伸缩振动一个峰．     

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.     

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.     

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.     

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.     

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.     

Guest-Host Interaction Study in Clathrate Hydrates Using Lattice Dynamics Simulation

Lattice dynamics simulation of several gas hydrates （helium, argon, and methane） with different occupancy rates has been performed using TIP3P potential model. Results show that the coupling between the guest and host is not simple as depicted by the conventional viewpoints. For clathrate hydrate enclosing small guest, the small cages are dominantly responsible for the thermodynamic stability of clathrate hydrates. And the spectrum of methane hydrate is studied compared with argon hydrate, then as a result, shrink effect from positive hydrogen shell is proposed.     

Competitive cage occupancy of hydrogen and argon in structure-II hydrates

The cage occupancy of hydrogen in the single-crystals of simple hydrogen hydrates and hydrogen + argon mixed-gas hydrates was investigated by means of in situ Raman spectroscopy under the three-phase (hydrate + water + fluid) equilibrium condition. In the equilibrium pressure region higher than approximately 25 MPa, four hydrogen cluster and argon competitively occupied the large cages of structure-II hydrogen + argon mixed-gas hydrates. In addition, Raman spectroscopic analysis at liquid nitrogen temperature (77 K) supports that the clusters of two, three, or four hydrogen molecules occupy large cages.     

Phase relations and binary clathrate hydrate formation in the system H2-THF-H2O

Experimentally determined equilibrium phase relations are reported for the system H2-THF-H2O as a function of aqueous tetrahydrofuran (THF) concentration from 260 to 290 K at pressures up to 45 MPa. Data are consistent with the formation of cubic structure-II (CS-II) binary H2-THF clathrate hydrates with a stoichiometric THF-to-water ratio of 1:17, which can incorporate modest volumes of molecular hydrogen at elevated pressures. Direct compositional analyses of the clathrate phase, at both low (0.20 mol %) and stoichiometric (5.56 mol %) initial THF aqueous concentrations, are consistent with observed phase behavior, suggesting full occupancy of large hexakaidecahedral (51264) clathrate cavities by THF, coupled with largely complete (80-90%) filling of small dodecahedral (512) cages by single H2 molecules at pressures of >30 MPa, giving a clathrate formula of (H2) < or =2.THF.17H2O. Results should help to resolve the current controversy over binary H2-THF hydrate hydrogen contents; data confirm recent reports that suggest a maximum of approximately 1 mass % H2, this contradicting values of up to 4 mass % previously claimed for comparable conditions.     

High-pressure gas hydrates

It has long been known that crystalline hydrates are formed by many simple gases that do not interact strongly with water, and in most cases the gas molecules or atoms occupy 'cages' formed by a framework of water molecules. The majority of these gas hydrates adopt one of two cubic cage structures and are called clathrate hydrates. Notable exceptions are hydrogen and helium which form 'exotic' hydrates with structures based on ice structures, rather than clathrate hydrates, even at low pressures. Clathrate hydrates have been extensively studied because they occur widely in nature, have important industrial applications, and provide insight into water-guest hydrophobic interactions. Until recently, the expectation-based on calculations-had been that all clathrate hydrates were dissociated into ice and gas by the application of pressures of 1 GPa or so. However, over the past five years, studies have shown that this view is incorrect. Instead, all the systems so far studied undergo structural rearrangement to other, new types of hydrate structure that remain stable to much higher pressures than had been thought possible. In this paper we review work on gas hydrates at pressures above 0.5 GPa, identify common trends in transformations and structures, and note areas of uncertainty where further work is needed.     

Vibrational spectra of methane clathrate hydrates from molecular dynamics simulation

Molecular dynamics simulations were performed on methane clathrate hydrates at ambient conditions. Thermal expansion results over the temperature range 60-300 K show that the unit cell volume increases with temperature in agreement with experiment. Power spectra were obtained at 273 K from velocity autocorrelation functions for selected atoms, and normal modes were assigned. The spectra were further classified according to individual atom types, allowing the assignment of contributions from methane molecules located in small and large cages within the structure I unit cell. The symmetric C-H stretch of methane in the small cages occurs at a higher frequency than for methane located in the large cages, with a peak separation of 14 cm(-1). Additionally, we determined that the symmetric C-H stretch in methane gas occurs at the same frequency as methane in the large cages. Results of molecular dynamics simulations indicate the use of power spectra obtained from the velocity autocorrelation function is a reliable method to investigate the vibrational behavior of guest molecules in clathrate hydrates.     

Molecular vibrations of methane molecules in the structure I clathrate hydrate from ab initio molecular dynamics simulation

Vibrational frequencies of guest molecules in clathrate hydrates reflect the molecular environment and dynamical behavior of molecules. A detailed understanding of the mechanism for the vibrational frequency changes of the guest molecules in the clathrate hydrate cages is still incomplete. In this study, molecular vibrations of methane molecules in a structure I clathrate hydrate are calculated from ab initio molecular dynamics simulation. The vibrational spectra of methane are computed by Fourier transform of autocorrelation functions, which reveal distinct separation of each vibrational mode. Calculated symmetric and asymmetric stretching vibrational frequencies of methane molecules are lower in the large cages than in the small cages (8 and 16 cm(-1) for symmetric and asymmetric stretching, respectively). These changes are closely linked with the C-H bond length. The vibrational frequencies for the bending and rocking vibrational modes nearly overlap in each of the cages.     

CS-II binary clathrate hydrates at pressures of up to 15 kbar

The pressure dependence of the decomposition temperatures of binary clathrate hydrates of tetra-hydrofuran with xenon and methane as well as of chloroform and carbon tetrachloride clathrates with xenon has been studied. The absence of phase transitions at pressures from 1 to 15,000 bar indicates that the structure of all the hydrates remains constant (CS-II). The decomposition temperatures of the binary hydrates of tetrahydrofuran and carbon tetrachloride with xenon at 15 kbar (above 124ℴC) are exceedingly high for polyhedral clathrate hydrates because the guest molecules are highly complementary to the cavities of the clathrate lattice. The paper also considers the packing density effect in the crystal structure of hydrates on the behavior of the latter at elevated pressure. Translated fromZhurnal Strukturnoi Khimii, Vol. 41, No. 3, pp. 582-589, May–June, 2000.     

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.     

Molecular dynamics simulations of binary structure H hydrogen and methyl-tert-butylether clathrate hydrates

Binary structure H (sH) hydrogen and methyl-tert-butylether (MTBE) clathrate hydrates are studied with molecular dynamics simulations. Simulations on a 3 x 3 x 3 sH unit cell with up to 4.7 mass % hydrogen gas are run at pressures of 100 bars and 2 kbars at 100 and 273 K. For the small and medium cages of the sH unit cell, H2 guest molecule occupancies of 0, 1 (single occupancy), and 2 (double occupancy) are considered with the MTBE molecule occupying all of the large cages. An increase of the small and medium cage occupancies from 1 to 2 leads to a jump in the unit cell volume and configurational energy. Calculations are also set up with 13, 23, and 89 of the MTBE molecules in the large cages replaced by sets of three to six H2 molecules, and the effects on the configurational energy and volume of the simulation cell are determined. As MTBE molecules are replaced with sets of H2 guests in the large cages, the configurational energy of the unit cell increases. At the lower temperature, the energy and volume of the clathrate are not sensitive to the number of hydrogen guests in the large cages; however, at higher temperatures the repulsions among the H2 guest molecules in the large cages cause an increase in the system energy and volume.     

Phase equilibrium measurements and the tuning behavior of new sII clathrate hydrates

We suggest two types of new amine-type sII formers: pyrrolidine and piperidine. These guest compounds fail to form clathrate hydrate structures with host water, but instead have to combine with light gaseous guest molecules (methane) for enclathration. First, two binary clathrate hydrates of (pyrrolidine + methane) and (piperidine + methane) were synthesized at various amine concentrations. ¹³C NMR and Raman analysis were done to identify the clathrate hydrate structure and guest distribution over sII-S and sII-L cages. XRD was also used to find the exact structure and corresponding cell parameters. At a dilute pyrrolidine concentration of less than 5.56 mol%, the tuning phenomenon is observed such that methane molecules surprisingly occupy sII-L cages. At the critical guest concentration of about 0.1 mol%, the cage occupancy ratio reaches the maximum of approximately 0.5. At very dilute guest concentration below 0.1 mol%, the methane molecules fail to occupy large cages on account of their rarefied distribution in the network. Direct-release experiments were performed to determine the actual guest compositions in the clathrate hydrate phases. Finally, we measured the clathrate hydrate phase equilibria of (pyrrolidine + methane) and (piperidine + methane).     

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.     

Hydrogen NMR of H2-TDF-D2O clathrate

The three-component clathrate H2-TDF-D2O offers hydrogen storage at lower pressure, but with reduced weight fraction of H2, compared to H2-H2O clathrate. In H2-TDF-D2O, H2 resides exclusively and singly in the small cages of structure II, allowing the rotational behavior of H2 in this nominally uniform environment to be probed. Here we report NMR measurements of the H2 line shape and relaxation times T1, T2, and T1rho. The principal differences in the results, compared to the H2-D2O binary system, are the dips in T2 and T1rho near 28 K due to thermally activated reorientation of TDF molecules, line-narrowing and decreases in T2 and T1rho near 175 K due to D2O reorientations and diffusion, and the apparent absence of H2 diffusion between small cages.