Molecular motion and phase transition in hydroquinone clathrate compounds

Distinction between the host-guest and guest-guest interaction in the roles they play in determining the physical properties of the clathrate compounds is discussed. It is shown that the latter causes cooperative effects and phase transitions. Experimental data on the phase transitions and molecular motion in the hydroquinone clathrate compounds are reviewed. Dipole interaction between the guest molecules is shown to have a correct magnitude of energy to explain the experimentally found transition temperatures. Possibility of quantum effects in the clathrate property is discussed in relation to the free rotation and ortho-para conversion of the hydrogen sulphide and other guest molecules.Contribution No. 80 from Chemical Thermodynamics Laboratory     

Molecular hydrogen storage in binary THF-H2 clathrate hydrates

The hydrogen storage capacity of binary THF-H(2) clathrate hydrate has been determined as a function of formation pressure, THF composition, and time. The amount of hydrogen stored in the stoichiometric hydrate increases with pressure and exhibits asymptotic (Langmuir) behavior to approximately 1.0 wt % H(2). This hydrogen concentration corresponds to one hydrogen molecule occupying each of the small 5(12) cavities and one THF molecule in each large 5(12)6(4) cavity in the hydrate framework. Contrary to previous reports, hydrogen storage was not increased upon decreasing the THF concentration below the stoichiometric 5.6 mol % solution to 0.5 mol %, at constant pressure, even after one week. This provides strong evidence that THF preferentially occupies the large 5(12)6(4) cavity over hydrogen, for the range of experimental conditions tested. The maximum amount of hydrogen stored in this binary hydrate was about 1.0 wt % at moderate pressure (<60 MPa) and is independent of the initial THF concentration over the range of conditions tested.     

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.     

Thermodynamic studies of clathrate hydrates

Thermodynamic studies of clathrate hydrates, mainly of structures I and II, are considered in this review which is based on 147 references. There are two main subjects. The first is the host lattice energy and the guest-host interaction energy, both of these quantities being related to the enthalpy of dissociation and composition of the hydrates. The second subject concerns orientational ordering phenomena occurring in both host and guest, as reflected in the low temperature heat capacity. The classical theoretical treatment of clathrate formation has been reconsidered on the basis of recent experimental results. Particular emphasis has been given to orientational ordering since this topic is undoubtedly central to clarifying the nature of clathrate hydrates.Ausgehend von 147 Literaturangaben wurden in diesem Review thermodynamische Untersuchungen von Klathrathydraten hauptsächlich der Struktur I und II betrachtet. Es gibt zwei Hauptaugenmerke. Als erstes die Wirtsgitterenergie und die Gast-Wirt-Wechselwirkungsenergie, beide bezogen auf die Dissoziationsenthalpie und die Bildungsenthalpie der Hydrate. Das zweite Hauptaugenmerk betrifft Orientierungs-Konditionierungserscheinungen sowohl in Wirt als auch Gast, wie in den Wärmekapazitäten bei niedrigen Temperaturen widergespiegelt wird. Auf der Basis jüngster experimenteller Ergebnisse wurde die klassische theoretische Betrachtung über die Bildung von Klathraten überprüft. Der Orientierung-Konditionierung wurde besonderer Nachdruck verliehen, da dies zweifellos eine entscheidende Rolle bei der Klärung der Natur der Klathrathydrate spielt. 147 I II. . «» « — », . «» « », . . , .

---

Contribution No. 155 from the Chemical Thermodynamics Laboratory.     

Structural properties and phase transitions in a silica clathrate

Melanophlogite, a low-pressure silica polymorph, has been extensively studied at different temperatures and pressures by molecular dynamics simulations. While the high-temperature form is confirmed as cubic, the low-temperature phase is found to be slightly distorted, in agreement with experiments. With increasing pressure, the crystalline character is gradually lost. At 8 GPa, the radial distribution function is consistent with an amorphous state. Like pristine glass, the topology changes, plastic behavior, and permanent densification appear above ～12 GPa, triggered by Si coordination number changes. We predict that a partial crystalline and amorphous sample can be obtained by recovering the sample from a pressure of ～12-16 GPa.     

Molecular dynamics study of the stability of methane structure H clathrate hydrates

Molecular dynamics simulations are used to study the stability of structure H (sH) methane clathrate hydrates in a 3 x 3 x 3 sH unit cell replica. Simulations are performed at experimental conditions of 300 K and 2 GPa for three methane intermolecular potentials. The five small cages of the sH unit cell are assigned methane guest occupancies of one and large cage guest occupancies of one to five are considered. Radial distribution functions, unit cell volumes, and configurational energies are studied as a function of large cage CH(4) occupancy. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies. Large cage occupancy of five is the most stable configuration for a Lennard-Jones united-atom potential and the Tse-Klein-McDonald potential parametrized for condensed methane phases and two for the most stable configuation for the Murad and Gubbins potential.     

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.     

Dynamic and thermodynamic properties of clathrate hydrates

Vibrational spectra and thermodynamic properties of ices and the cubic structure I (CS-I) clathrate hydrate have been studied by the lattice dynamics method. The phonon density of states for the empty hydrate framework and for xenon hydrate have been determined; the vibrational frequencies of the guest molecules in large and small cavities have been found. The stability of the hydrate with respect to the external pressure at low temperatures and its thermodynamic stability at temperatures around 0°C have been studied. It has been found that the empty hydrate framework is unstable in certain temperature and pressure regions. A definite degree of occupation of the large cavities by the guest molecules is necessary for the hydrate to become stable. It has been found that there is a maximum of the critical temperature at which the hydrate exists, which is a function of the external pressure.Dedicated to Dr. W. Davidson in honor of his great contributions to the sciences of inclusion phenomena.     

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.     

Polymorphism in Br2 clathrate hydrates

The structure and composition of bromine clathrate hydrate has been controversial for more than 170 years due to the large variation of its observed stoichiometries. Several different crystal structures were proposed before 1997 when Udachin et al. (Udachin, K. A.; Enright, G. D.; Ratcliffe, C. I.; Ripmeester, J. A. J. Am. Chem. Soc. 1997, 119, 11481) concluded that Br2 forms only the tetragonal structure (TS-I). We show polymorphism in Br2 clathrate hydrates by identifying two distinct crystal structures through optical microscopy and resonant Raman spectroscopy on single crystals. After growing TS-I crystals from a liquid bromine-water solution, upon dropping the temperature slightly below -7 degrees C, new crystals of cubic morphology form. The new crystals, which have a limited thermal stability range, are assigned to the CS-II structure. The two structures are clearly distinguished by the resonant Raman spectra of the enclathrated Br2, which show long overtone progressions and allow the extraction of accurate vibrational parameters: omega(e) = 321.2 +/- 0.1 cm(-1) and omega(e)x(e) = 0.82 +/- 0.05 cm(-1) in TS-I and omega(e) = 317.5 +/- 0.1 cm(-1) and omega(e)x(e) = 0.70 +/- 0.1 cm(-1) in CS-II. On the basis of structural analysis, the discovery of the CS-II crystals implies stability of a large class of bromine hydrate structures and, therefore, polymorphism.     

Structural and dynamical properties of methane clathrate hydrates

Equilibrium molecular dynamics (MD) simulations have been performed in both the NVT and NPT ensembles to study the structural and dynamical properties of fully occupied methane clathrate hydrates at 50, 125, and 200 K. Five atomistic potential models were used for water, ranging from fully flexible to rigid polarizable and nonpolarizable. A flexible and a rigid model were utilized for methane. The phonon densities of states were evaluated and the localized rattling modes for the methane molecules were found to couple to the acoustic phonons of the host lattice. The calculated methane density of states was found to be in reasonable agreement with available experimental data.     

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.     

Thermodynamic modeling for clathrate hydrates of ozone

We report a theoretical study to predict the phase-equilibrium properties of ozone-containing clathrate hydrates based on the statistical thermodynamics model developed by van der Waals and Platteeuw. The Patel–Teja–Valderrama equation of state is employed for an accurate estimation of the properties of gas phase ozone. We determined the three parameters of the Kihara intermolecular potential for ozone as a = 6.815 · 10⁻² nm, σ = 2.9909 · 10⁻¹ nm, and ε · k_B⁻¹ = 184.00 K. An infinite set of ε–σ parameters for ozone were determined, reproducing the experimental phase equilibrium pressure–temperature data of the (O₃ + O₂ + CO₂) clathrate hydrate. A unique parameter pair was chosen based on the experimental ozone storage capacity data for the (O₃ + O₂ + CCl₄) hydrate that we reported previously. The prediction with the developed model showed good agreement with the experimental phase equilibrium data within ±2% of the average deviation of the pressure. The Kihara parameters of ozone showed slightly better suitability for the structure-I hydrate than CO₂, which was used as a help guest. Our model suggests the possibility of increasing the ozone storage capacity of clathrate hydrates (∼7% on a mass basis) from the previously reported experimental capacity (∼1%).     

Localized oscillators and heat conduction in clathrate hydrates

The glassy behaviour of the thermal conductivities of structure I xenon and structure II hydrate at high temperature are found to be described surprisingly well by a localized oscillator model. This observation leads to the suggestion of strong coupling between the localized vibrations of the guest with the lattice acoustic phonons. This conjecture is confirmed by a phenomenological calculation using the Anderson-Fano resonant scattering model.Published as NRCC 37248.     

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.     

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.