Structural Investigations of Argon Hydrates at Pressures up to 10 kbar

Neutron diffraction patterns of three argon hydrates which exist at the pressures up to 10 kbar has been studied; Rietveld refinement of their structures has been done. The phase which is stable from 1 bar to 4.6 kbar appears to be typical cubic structure II gas hydrate with variable degree of filling of the large cavities. Stoichiometry of this compound under high-pressure conditions has been determined for the first time and appears to be ArW4.5H₂O and Ar·4H₂O at 3.4 and 4.3 kbar, respectively. Gas hydrate existing in the pressure range of 4.6–7.7 kbar has a hexagonal structure (hexagonal structure III, so-called structure H). Refinement of the structure has shown that the best agreement between calculated and experimental pattern can be reached in the case of accommodation of five (!) argon atoms in the large cavity. Indexing of the neutron diffraction pattern of the hydrate stable in the 7.7–9.5 kbar range leads to the primitive tetragonal unit cell with parameters a = 6.342 Å, c = 10.610 Å at 9.2 kbar, which does not correspond to any known type of gas hydrates. The water framework of this structure was found by idealizing the structure of pinacol semiclathrate hydrate. This hydrate belongs to a new, earlier unknown, tetragonal structural type of gas hydrates. It contains only one type of polyhedral cavities with 14 faces. This type of polyhedrons are space-filling; two argon atoms occupy each cavity. This structure gives the first example of the gas hydrate water framework which contains only one type of polyhedral cavities.     

First determination of volume changes and enthalpies of the high-pressure decomposition reaction of the structure H methane hydrate to the cubic structure I methane hydrate and fluid methane

Pressure-temperature (P-T) conditions of the decomposition reaction of the structure H high-pressure methane hydrate to the cubic structure I methane hydrate and fluid methane were studied with a piston-cylinder apparatus at room temperature. For the first time, volume changes accompanying this reaction were determined. With the use of the Clausius-Clapeyron equation the enthalpies of the decomposition reaction of the structure H high-pressure methane hydrate to the cubic structure I methane hydrate and fluid methane have been calculated.     

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.     

Hydrate kinetics study in the presence of nonaqueous liquid by nuclear magnetic resonance spectroscopy and imaging

The dynamics of methane hydrate growth and decomposition were studied by nuclear magnetic resonance (NMR) spectroscopy and imaging (MRI). Three well-known large molecule guest substances (LMGS) were used as structure H hydrate formers: 2,2-dimethylbutane (NH), methylcyclohexane (MCH), tert-butyl methyl ether (TBME). In addition, the impact of a non-hydrate former (n-heptane/nC7) was studied. The methane diffusion and hydrate growth were monitored by recording the 2H NMR spectra at 253 K and approximately 4.5 MPa for 20 h. The results revealed that methane diffuses faster in TBME and NH, slower in nC7, and slowest in MCH. The TBME system gives the fastest hydrate formation kinetics followed by NH, MCH, and nC7. The conversion of water into hydrate was also observed. The imaging study showed that TBME has a strong affinity toward ice, which is not the case for the NH and MCH systems. The degree of ice packing was also found to affect the LMGS distribution between ice particles. Highly packed ice increases the mass transfer resistance and hence limits the contact between LMGS and ice. It was also found that "temperature ramping" above the ice point improves the conversion significantly. Finally, hydrates were found to dissociate quickly within the first hour at atmospheric pressure and subsequently at a much slower rate. Methane dissolved in LMGS was also seen. The residual methane in hydrate phase and dissolved in LMGS phase explain the faster kinetics during hydrate re-formation.     

Experimental and modeling study on decomposition kinetics of methane hydrates in different media

The decomposition kinetic behaviors of methane hydrates formed in 5 cm3 porous wet activated carbon were studied experimentally in a closed system in the temperature range of 275.8-264.4 K. The decomposition rates of methane hydrates formed from 5 cm3 of pure free water and an aqueous solution of 650 g x m(-3) sodium dodecyl sulfate (SDS) were also measured for comparison. The decomposition rates of methane hydrates in seven different cases were compared. The results showed that the methane hydrates dissociate more rapidly in porous activated carbon than in free systems. A mathematical model was developed for describing the decomposition kinetic behavior of methane hydrates below ice point based on an ice-shielding mechanism in which a porous ice layer was assumed to be formed during the decomposition of hydrate, and the diffusion of methane molecules through it was assumed to be one of the control steps. The parameters of the model were determined by correlating the decomposition rate data, and the activation energies were further determined with respect to three different media. The model was found to well describe the decomposition kinetic behavior of methane hydrate in different media.     

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.     

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.     

Activation energy of methyl radical decay in methane hydrate

The thermal stability of gamma-ray-induced methyl radicals in methane hydrate was studied using the ESR method at atmospheric pressure and 210-260 K. The methyl radical decay proceeded with the second-order reaction, and ethane molecules were generated from the dimerization process. The methyl radical decay proceeds by two different temperature-dependent processes, that is, the respective activation energies of these processes are 20.0 +/- 1.6 kJ/mol for the lower temperature region of 210-230 K and 54.8 +/- 5.7 kJ/mol for the higher temperature region of 235-260 K. The former agrees well with the enthalpy change of methane hydrate dissociation into ice and gaseous methane, while the latter agrees well with the enthalpy change into liquid water and gaseous methane. The present findings reveal that methane hydrates dissociate into liquid (supercooled) water and gaseous methane in the temperature range of 235-260 K.     

Mechanisms for thermal conduction in hydrogen hydrate

Extensive equilibrium molecular dynamics simulations have been performed to investigate thermal conduction mechanisms via the Green-Kubo approach for (type II) hydrogen hydrate, at 0.05 kbar and between 30 and 250 K, for both lightly filled H(2) hydrates (1s4l) and for more densely filled H(2) systems (2s4l), in which four H(2) molecules are present in the large cavities, with respective single- and double-occupation of the small cages. The TIP4P water model was used in conjunction with a fully atomistic hydrogen potential along with long-range Ewald electrostatics. It was found that substantially less damping in guest-host energy transfer is present in hydrogen hydrate as is observed in common type I clathrates (e.g., methane hydrate), but more akin in to previous results for type II and H methane hydrate polymorphs. This gives rise to larger thermal conductivities relative to common type I hydrates, and also larger than type II and H methane hydrate polymorphs, and a more crystal-like temperature dependence of the thermal conductivity.     

Phase Transition of a Structure II Cubic Clathrate Hydrate to a Tetragonal Form

The crystal structure and phase transition of cubic structure II (sII) binary clathrate hydrates of methane (CH₄) and propanol are reported from powder X‐ray diffraction measurements. The deformation of host water cages at the cubic–tetragonal phase transition of 2‐propanol+CH₄ hydrate, but not 1‐propanol+CH₄ hydrate, was observed below about 110 K. It is shown that the deformation of the host water cages of 2‐propanol+CH₄ hydrate can be explained by the restriction of the motion of 2‐propanol within the 5¹²6⁴ host water cages. This result provides a low‐temperature structure due to a temperature‐induced symmetry‐lowering transition of clathrate hydrate. This is the first example of a cubic structure of the common clathrate hydrate families at a fixed composition.     

Molecular simulation of gas hydrate nanoclusters in water shell: Structure and phase transitions

Gas hydrate nanoclusters surrounded by water shells are studied by the molecular dynamic method. Hydrates of methane (sI structures) and krypton (sII structures), as well an ice nanocluster in a supercooled water shell, are considered. The main attention was focused on studying the local structure and phase transitions. Variations in local partial densities with an increase in temperature are monitored. Melting points of nanosized samples of gas hydrates are determined using caloric curves. Additional information on the behavior of the considered systems is obtained from the temperature dependences of diffusion coefficients and the Lindemann criterion. Two-phase transitions are revealed for gas hydrate nanoclusters. The first phase transition at 210 K can be assigned to the melting of the ice shell. The second transition at 230–235 K is identified as the phase transition in the hydrate core. The melting of ice cluster is observed at 215 K, which corresponds to the melting point of bulk crystal upon the use of the SPC/E water model.     

Modulated DSC for gas hydrates analysis

Modulated DSC has been applied to the study of methane, ethane and propane hydrates at different hydrate and ice concentrations. The reversing component of the TMDSC curves, makes it possible to characterize such hydrates. Methane and ethane hydrates show the melting-decomposition peak at a temperatures higher than the ice contained in the sample, while propane hydrate melts and decomposes at lower temperature than the ice present in the sample. The hydrate peaks tend to disappear if the hydrate is stored at atmospheric pressure. Guest size and cavity occupation fix the heat of dissociation and stability of the hydrates, as confirmed by parallel tests on tetrahydrofurane hydrates.     

Raman spectra of gas hydrates--differences and analogies to ice 1h and (gas saturated) water

It is generally accepted that Raman spectroscopic investigations of gas hydrates provide vital information regarding the structure of the hydrate, hydrate composition and cage occupancies, but most research is focused on the vibrational spectra of the guest molecules. We show that the shape and position of the Raman signals of the host molecules (H(2)O) also contain useful additional information. In this study, Raman spectra (200-4000 cm(-1)) of (mixed) gas hydrates with variable compositions and different structures are presented. The bands in the OH stretching region (3000-3800 cm(-1)), the O-H bending region (1600-1700 cm(-1)) and the O-O hydrogen bonded stretching region (100-400 cm(-1)) are compared with the corresponding bands in Raman spectra of ice Ih and liquid water. The interpretation of the differences and similarities with respect to the crystal structure and possible interactions between guest and host molecules are presented.     

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.     

Thermodynamic calculations in the system CH4-H2O and methane hydrate phase equilibria

Using the Gibbs function of reaction, equilibrium pressure, temperature conditions for the formation of methane clathrate hydrate have been calculated from the thermodynamic properties of phases in the system CH4-H2O. The thermodynamic model accurately reproduces the published phase-equilibria data to within +/-2 K of the observed equilibrium boundaries in the range 0.08-117 MPa and 190-307 K. The model also provides an estimate of the third-law entropy of methane hydrate at 273.15 K, 0.1 MPa of 56.2 J mol(-1) K(-1) for 1/nCH4.H2O, where n is the hydrate number. Agreement between the calculated and published phase-equilibria data is optimized when the hydrate composition is fixed and independent of the pressure and temperature for the conditions modeled.     

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.     

Monte Carlo and molecular dynamics simulations of methane in potassium montmorillonite clay hydrates at elevated pressures and temperatures

The structure and dynamics of methane in hydrated potassium montmorillonite clay have been studied under conditions encountered in sedimentary basin and compared to those of hydrated sodium montmorillonite clay using computer simulation techniques. The simulated systems contain two molecular layers of water and followed gradients of 150 bar km(-1) and 30 K km(-1) up to a maximum burial depth of 6 km. Methane particle is coordinated to about 19 oxygen atoms, with 6 of these coming from the clay surface oxygen. Potassium ions tend to move away from the center towards the clay surface, in contrast to the behavior observed with the hydrated sodium form. The clay surface affinity for methane was found to be higher in the hydrated K-form. Methane diffusion in the two-layer hydrated K-montmorillonite increases from 0.39 x 10(-9) m2 s(-1) at 280 K to 3.27 x 10(-9) m2 s(-1) at 460 K compared to 0.36 x 10(-9) m2 s(-1) at 280 K to 4.26 x 10(-9) m2 s(-1) at 460 K in Na-montmorillonite hydrate. The distributions of the potassium ions were found to vary in the hydrates when compared to those of sodium form. Water molecules were also found to be very mobile in the potassium clay hydrates compared to sodium clay hydrates.     

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.     

The specific surface area of methane hydrate formed in different conditions and manners

The specific surface area of methane hydrates, formed both in the presence and absence of sodium dodecyl sulfate (SDS) and processed in different manners (stirring, compacting, holding the hydrates at the formation conditions for different periods of time, cooling the hydrates for different periods of time before depressurizing them), was measured under atmospheric pressure and temperatures below ice point. It was found that the specific surface area of hydrate increased with the decreasing temperature. The methane hydrate in the presence of SDS was shown to be of bigger specific surface areas than pure methane hydrates. The experimental results further demonstrated that the manners of forming and processing hydrates affected the specific surface area of hydrate samples. Stirring or compacting made the hydrate become finer and led to a bigger specific surface area. Supported by the National natural Science Foundation of China (Grant Nos.20490207, 2076145, uo633003), Program for New Century Excellent Talents in University and National The National High Technology Research and Development Program of China Project.     

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.