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.     

Gas hydrates of argon and methane synthesized at high pressures: composition, thermal expansion, and self-preservation

For the first time, the compositions of argon and methane high-pressure gas hydrates have been directly determined. The studied samples of the gas hydrates were prepared under high-pressure conditions and quenched at 77 K. The composition of the argon hydrate (structure H, stable at 460-770 MPa) was found to be Ar.(3.27 +/- 0.17)H(2)O. This result shows a good agreement with the refinement of the argon hydrate structure using neutron powder diffraction data and helps to rationalize the evolution of hydrate structures in the Ar-H(2)O system at high pressures. The quenched argon hydrate was found to dissociate in two steps. The first step (170-190 K) corresponds to a partial dissociation of the hydrate and the self-preservation of a residual part of the hydrate with an ice cover. Presumably, significant amounts of ice Ic form at this stage. The second step (210-230 K) corresponds to the dissociation of the residual part of the hydrate. The composition of the methane hydrate (cubic structure I, stable up to 620 MPa) was found to be CH(4).5.76H(2)O. Temperature dependence of the unit cell parameters for both hydrates has been also studied. Calculated from these results, the thermal expansivities for the structure H argon hydrate are alpha(a) = 76.6 K(-1) and alpha(c) = 77.4 K(-1) (in the 100-250 K temperature range) and for the cubic structure I methane hydrate are alpha(a) = 32.2 K(-1), alpha(a) = 53.0 K(-1), and alpha(a) = 73.5 K(-1) at 100, 150, and 200 K, respectively.     

On the thermodynamic stability and structural transition of clathrate hydrates

Gas mixtures of methane and ethane form structure II clathrate hydrates despite the fact that each of pure methane and pure ethane gases forms the structure I hydrate. Optimization of the interaction potential parameters for methane and ethane is attempted so as to reproduce the dissociation pressures of each simple hydrate containing either methane or ethane alone. An account for the structural transitions between type I and type II hydrates upon changing the mole fraction of the gas mixture is given on the basis of the van der Waals and Platteeuw theory with these optimized potentials. Cage occupancies of the two kinds of hydrates are also calculated as functions of the mole fraction at the dissociation pressure and at a fixed pressure well above the dissociation pressure.     

The conversion process of hydrocarbon hydrates into CO2 hydrates and vice versa: thermodynamic considerations

Microscopy, confocal Raman spectroscopy and powder X-ray diffraction (PXRD) were used for in situ investigations of the CO(2)-hydrocarbon exchange process in gas hydrates and its driving forces. The study comprises the exposure of simple structure I CH(4) hydrate and mixed structure II CH(4)-C(2)H(6) and CH(4)-C(3)H(8) hydrates to gaseous CO(2) as well as the reverse reaction, i.e., the conversion of CO(2)-rich structure I hydrate into structure II mixed hydrate. In the case of CH(4)-C(3)H(8) hydrates, a conversion in the presence of gaseous CO(2) from a supposedly more stable structure II hydrate to a less stable structure I CO(2)-rich hydrate was observed. PXRD data show that the reverse process requires longer initiation times, and structural changes seem to be less complete. Generally, the exchange process can be described as a decomposition and reformation process, in terms of a rearrangement of molecules, and is primarily induced by the chemical potential gradient between hydrate phase and the provided gas phase. The results show furthermore the dependency of the conversion rate on the surface area of the hydrate phase, the thermodynamic stability of the original and resulting hydrate phase, as well as the mobility of guest molecules and formation kinetics of the resulting hydrate phase.     

Investigation of hydrate formation in the system H2-CH4-H2O at a pressure up to 250 MPa

Phase equilibria in the system H2-CH4-H2O are investigated by means of differential thermal analysis within hydrogen concentration range 0-70 mol % and at a pressure up to 250 MPa. All the experiments were carried out under the conditions of gas excess. With an increase in hydrogen concentration in the initial gas mixture, decomposition temperature of the formed hydrates decreased. X-ray diffraction patterns and Raman spectra of the quenched hydrate samples obtained at a pressure of 20 MPA from a gas mixture containing 40 mol % hydrogen were recorded. It turned out that the hydrate has cubic structure I under these conditions. The Raman spectra showed that hydrogen molecules are not detected in the hydrate within the sensitivity of the method, that is, almost pure methane hydrate is formed. The general view of the phase diagram of the investigated system is proposed. A thermodynamic model was proposed to explain a decrease in hydrate decomposition temperature in the system with an increase in the concentration of hydrogen in the initial mixture.     

Dissociation enthalpies of synthesized multicomponent gas hydrates with respect to the guest composition and cage occupancy

This study presents the influences of additional guest molecules such as C2H6, C3H8, and CO2 on methane hydrates regarding their thermal behavior. For this purpose, the onset temperatures of decomposition as well as the enthalpies of dissociation were determined for synthesized multicomponent gas hydrates in the range of 173-290 K at atmospheric pressure using a Calvet heat-flow calorimeter. Furthermore, the structures and the compositions of the hydrates were obtained using X-ray diffraction and Raman spectroscopy as well as hydrate prediction program calculations. It is shown that the onset temperature of decomposition of both sI and sII hydrates tends to increase with an increasing number of larger guest molecules than methane occupying the large cavities. The results of the calorimetric measurements also indicate that the molar dissociation enthalpy depends on the guest-to-cavity size ratio and the actual concentration of the guest occupying the large cavities of the hydrate. To our knowledge, this is the first study that observes this behavior using calorimetrical measurements on mixed gas hydrates at these temperature and pressure conditions.     

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.     

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.     

Some structural and thermodynamic studies of clathrate hydrates

X-ray and neutron diffraction studies show argon and krypton to preferentially form clathrate hydrates of structure II, rather than structure I as previously assumed; methane and hydrogen sulphide do form structure I. Re-examination of solid-solution thermodynamic theory shows that structure II is basically the more stable; structure I is generally formed only when the guest molecule is in the size range that favours occupancy of the 14-hedral over the 12-hedral cages. For molecules too large to enter the 12-hedra the relative stability of structure II is greatest at 0°C, in agreement with the observed sequence of change of stability of cyclopropane hydrate: I to II at –16° and II to I at 1.5°. Carbon dioxide hydrate is observed to decompose on prolonged standing at 105K in accord with the low-temperature instability predicted by Miller.     

Application of the cell potential method to predict phase equilibria of multicomponent gas hydrate systems

We present the application of a mathematical method reported earlier by which the van der Waals-Platteeuw statistical mechanical model with the Lennard-Jones and Devonshire approximation can be posed as an integral equation with the unknown function being the intermolecular potential between the guest molecules and the host molecules. This method allows us to solve for the potential directly for hydrates for which the Langmuir constants are computed, either from experimental data or from ab initio data. Given the assumptions made in the van der Waals-Platteeuw model with the spherical-cell approximation, there are an infinite number of solutions; however, the only solution without cusps is a unique central-well solution in which the potential is at a finite minimum at the center to the cage. From this central-well solution, we have found the potential well depths and volumes of negative energy for 16 single-component hydrate systems: ethane (C2H6), cyclopropane (C3H6), methane (CH4), argon (Ar), and chlorodifluoromethane (R-22) in structure I; and ethane (C2H6), cyclopropane (C3H6), propane (C3H8), isobutane (C4H10), methane (CH4), argon (Ar), trichlorofluoromethane (R-11), dichlorodifluoromethane (R-12), bromotrifluoromethane (R-13B1), chloroform (CHCl3), and 1,1,1,2-tetrafluoroethane (R-134a) in structure II. This method and the calculated cell potentials were validated by predicting existing mixed hydrate phase equilibrium data without any fitting parameters and calculating mixture phase diagrams for methane, ethane, isobutane, and cyclopropane mixtures. Several structural transitions that have been determined experimentally as well as some structural transitions that have not been examined experimentally were also predicted. In the methane-cyclopropane hydrate system, a structural transition from structure I to structure II and back to structure I is predicted to occur outside of the known structure II range for the cyclopropane hydrate. Quintuple (L(w)-sI-sII-L(hc)-V) points have been predicted for the ethane-propane-water (277.3 K, 12.28 bar, and x(eth,waterfree) = 0.676) and ethane-isobutane-water (274.7 K, 7.18 bar, and x(eth,waterfree) = 0.81) systems.     

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.     

I型和II型结构气体水合物的记忆效应

利用水合物二次生成实验装置, 采用“定容法”对I型(甲烷、二氧化碳)和II型(丙烷)结构气体水合物的二次生成进行了实验, 研究了不同结构水合物(I型、II型)彼此间的记忆效应, 发现水合物生成过程存在明显的诱导期, I型结构水合物间在二次生成过程中存在着记忆效应. I型与II型结构水合物之间在相互二次生成过程中存在着显著的记忆效应.     

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

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.     

用分子动力学模拟水合物储氢

用分子动力学(MD)模拟方法系统研究了结构Ⅰ型(SⅠ)和结构H型(SH)氢气水合物中氢气的占据情况并确定了氢气水合物的稳定结构: SⅠ水合物氢气分子数小胞中为2, 大胞中为3; SH水合物氢气分子数小胞中为2, 中胞中为2, 大胞中为11. 分析了稳定情况下水合物各胞腔内氢气分子之间的径向分布函数(RDF), 得出了氢气分子在各胞腔内的稳定位置. 由稳定位置得到了稳定结构下氢气水合物的储氢质量分数: SⅠ为5.085%, SH为6.467%. 与实验对比得出结论: SH水合物稳定结构下的储氢能力最强.     

Gas hydrate formation in the system C2H6–H2–H2O at pressures up to 250 MPa

Decomposition curves of gas hydrates formed in the ethane–hydrogen–water system were studied in the pressure interval 2–250 MPa. Gas hydrates synthesized at low (up to 5 MPa) pressures were also studied with use of X-ray powder diffraction and Raman spectroscopy. It was shown that ethane–hydrogen mixtures with hydrogen contents 0–30 mol.% form cubic structure I gas hydrates. Higher hydrogen concentration most probably results in appearance of another hydrate phase. We speculate that the gas mixtures with the hydrogen content above 60 mol.% form cubic structure II double hydrate of hydrogen and ethane at temperatures below ≈280 K and pressures above 25 MPa.     

Crystal lattice size and stability of type H clathrate hydrates with various large-molecule guest substances

To gain a better understanding of the effects of guest molecules on the lattice and stability of type H hydrates, we performed powder X-ray diffraction (PXRD) measurements and semiempirical molecular orbital calculations. The unit cell parameters and cohesive energies of various type H hydrates that contain methane (CH4) were analyzed. PXRD measurements indicated that an increase in the large-molecule guest volume caused the unit cell volume to increase. It was also indicated that a large-molecule guest substance caused the a-axis-direction of the unit cell to increase with little decrease in the c-axis direction. Calculations of cohesive energy by means of a semiempirical molecular orbital method indicated that the functional group and configuration of large-molecule guest substances affects the stability of type H hydrates. It was concluded that the icosahedron (5(12)6(8)) cages do not easily increase in length along the c-axis direction when larger guest molecules are used to form the hydrate, but the 5(12)6(8) cage and the layer of dodecahedron (5(12)) cages can easily increase in length along the a-axis direction due to interactions of the guest-host molecule.     

Crystal polymorphs of 6‐[(phenylcarbamoyl)oxy]hexa‐2,4‐diyn‐1‐yl isonicotinate

The title compound, C₁₉H₁₄N₂O₄, was found to have two crystal polymorphs, in which the molecular structures of the diacetylenic compound are broadly similar. The main structural difference between the polymorphs concerns the intermolecular hydrogen‐bonding motifs adopted, namely a one‐dimensional zigzag polymer linked by N—H…N(py) (py is pyridine) interactions in polymorph I and a centrosymmetric dimeric motif formed by N—H…O=C interactions in polymorph II. The diacetylene cores of the molecules stack along the a and b axes in polymorphs I and II, respectively. It was found that only the molecular arrangement in polymorph II satisfies Baughman's criterion to afford polydiacetylenes (PDAs) by thermal annealing or irradiation with light. This predicted polymerization activity was confirmed by experiment.     

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.     

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.     

First principles molecular dynamics study of filled ice hydrogen hydrate

We investigated structural changes, phase diagram, and vibrational properties of hydrogen hydrate in filled-ice phase C(2) by using first principles molecular dynamics simulation. It was found that the experimentally reported "cubic" structure is unstable at low temperature and∕or high pressure: The "cubic" structure reflects the symmetry at high (room) temperature where the hydrogen bond network is disordered and the hydrogen molecules are orientationally disordered due to thermal rotation. In this sense, the "cubic" symmetry would definitely be lowered at low temperature where the hydrogen bond network and the hydrogen molecules are expected to be ordered. At room temperature and below 30 GPa, it is the thermal effects that play an essential role in stabilizing the structure in "cubic" symmetry. Above 60 GPa, the hydrogen bonds in the framework would be symmetrized and the hydrogen bond order-disorder transition would disappear. These results also suggest the phase behavior of other filled-ice hydrates. In the case of rare gas hydrate, there would be no guest molecules' rotation-nonrotation transition since the guest molecules keep their spherical symmetry at any temperature. On the contrary methane hydrate MH-III would show complex transitions due to the lower symmetry of the guest molecule. These results would encourage further experimental studies, especially nuclear magnetic resonance spectroscopy and neutron scattering, on the phases of filled-ice hydrates at high pressures and∕or low temperatures.