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.     

On the thermodynamic stability of clathrate hydrates IV: double occupancy of cages

We have extended the van der Waals and Platteeuw theory to treat multiple occupancy of a single cage of clathrate hydrates, which has not been taken into account in the original theory but has been experimentally confirmed as a real entity. We propose a simple way to calculate the free energy of multiple cage occupancy and apply it to argon clathrate structure II in which a larger cage can be occupied by two argon atoms. The chemical potential of argon is calculated treating it as an imperfect gas, which is crucial to predict accurate pressure dependence of double occupancy expected at high pressure. It is found that double occupancy dominates over single occupancy when the guest pressure in equilibrium with the clathrate hydrate exceeds 270 MPa.     

On the thermodynamic stability of hydrogen clathrate hydrates

The cage occupancy of hydrogen clathrate hydrate has been examined by grand canonical Monte Carlo (GCMC) simulations for wide ranges of temperature and pressure. The simulations are carried out with a fixed number of water molecules and a fixed chemical potential of the guest species so that hydrogen molecules can be created or annihilated in the clathrate. Two types of the GCMC simulations are performed; in one the volume of the clathrate is fixed and in the other it is allowed to adjust itself under a preset pressure so as to take account of compression by a hydrostatic pressure and expansion due to multiple cage occupancy. It is found that the smaller cage in structure II is practically incapable of accommodating more than a single guest molecule even at pressures as high as 500 MPa, which agrees with the recent experimental investigations. The larger cage is found to encapsulate at most 4 hydrogen molecules, but its occupancy is dependent significantly on the pressure of hydrogen.     

Communication: Single crystal x-ray diffraction observation of hydrogen bonding between 1-propanol and water in a structure II clathrate hydrate

Single crystal x-ray crystallography is used to detect guest-host hydrogen bonding in structure II (sII) binary clathrate hydrate of 1-propanol and methane. X-ray structural analysis shows that the 1-propanol oxygen atom is at a distance of 2.749 and 2.788 ? from the closest clathrate hydrate water oxygen atoms from a hexagonal face of the large sII cage. The 1-propanol hydroxyl hydrogen atom is disordered and at distances of 1.956 and 2.035 ? from the closest cage water oxygen atoms. These distances are compatible with guest-water hydrogen bonding. The C-C-C-O torsional angle in 1-propanol in the cage is 91.47° which corresponds to a staggered conformation for the guest. Molecular dynamics studies of this system demonstrated guest-water hydrogen bonding in this hydrate. The molecular dynamics simulations predict most probable distances for the 1-propanol-water oxygen atoms to be 2.725 ?, and the average C-C-C-O torsional angle to be ~59° consistent with a gauche conformation. The individual cage distortions resulting from guest-host hydrogen bonding from the simulations are rather large, but due to the random nature of the hydrogen bonding of the guest with the 24 water molecules making up the hexagonal faces of the large sII cages, these distortions are not observed in the x-ray structure.     

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.     

How much carbon dioxide can be stored in the structure H clathrate hydrates?: a molecular dynamics study

The stability of structure H (sH) carbon dioxide clathrate hydrates at three temperature-pressure conditions are determined by molecular dynamics simulations on a 3x3x3 sH unit cell replica. Simulations are performed at 100 K at ambient pressure, 273 K at 100 bars and also 300 K and 5.0 kbars. The small and medium cages of the sH unit cell are occupied by a single carbon dioxide guest and large cage guest occupancies of 1-5 are considered. Radial distribution functions are given for guests in the large cages and unit cell volumes and configurational energies are studied as a function of large cage CO(2) occupancy. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies at three temperature/pressure conditions stated above. At the low temperature, large cage occupancy of 5 is the most stable while at the higher temperature, the occupancy of 3 is the most favored. Calculations are also performed to show that the CO(2) sH clathrate is more stable than the methane clathrate analog. Implications on CO(2) sequestration by clathrate formation are discussed.     

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.     

The interface between water and a hydrophobic gas

Classical molecular dynamics simulations have been performed to investigate the interface between liquid water and methane gas under methane hydrate forming conditions. The local environments of the water molecules were studied using order parameters which distinguish between liquid water, ice and methane hydrate phases. Bulk water and water/air interfaces were also studied to allow comparisons to be made between water molecules in the different environments and to determine the effects of the different methane densities studied. Good agreement between experimental and calculated surface tensions is obtained if long range corrections are included. The water surface is found to have a structure which is very similar to that of bulk water, but more tetrahedral, and more clathrate-like than ice-like. In these simulations the concentration of methane in water at the interface is shown to be appropriate for clathrates at higher gas densities (pressures). The orientation of water molecules around methane molecules in the interfacial region appears to depend only weakly on pressure and one of the difficulties in forming hydrate is the availability of water molecules tangential to the hydrate cage. At the interface, the water structure is more disordered than in the bulk water region with increased occurrence compared with the bulk of those angles and orientations found in the clathrate structure.     

Stability of rare gas structure H clathrate hydrates

Molecular dynamics simulations are used to study the stability of structure H (sH) clathrate hydrates with the rare gases Ne, Ar, Kr, and Xe. Simulations on a 3 x 3 x 3 sH unit cell replica are performed at ambient pressure at 40 and 100 K temperatures. The small and medium (s+m) cages of the sH unit cell are assigned rare gas guest occupancies of 1 and for large (l) cages guest occupancies of 1-6 are considered. Radial distribution functions for guest pairs with occupancies in the l-l, l-(s+m), and (s+m)-(s+m) cages are presented. The unit cell volumes and configurational energies are studied as a function of large cage occupancy for the rare gases. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies at 100 K and 1 bar and 20 kbar pressures. These studies show that the most stable argon clathrate has five guests in the large cages. For krypton and xenon the most stable configurations have three and two guests in the large cages, respectively.     

Phase changes of CO2 hydrate under high pressure and low temperature

High pressure and low temperature experiments with CO(2) hydrate were performed using diamond anvil cells and a helium-refrigeration cryostat in the pressure and temperature range of 0.2-3.0 GPa and 280-80 K, respectively. In situ x-ray diffractometry revealed that the phase boundary between CO(2) hydrate and water+CO(2) extended below the 280 K reported previously, toward a higher pressure and low temperature region. The results also showed the existence of a new high pressure phase above approximately 0.6 GPa and below 1.0 GPa at which the hydrate decomposed to dry ice and ice VI. In addition, in the lower temperature region of structure I, a small and abrupt lattice expansion was observed at approximately 210 K with decreasing temperature under fixed pressures. The expansion was accompanied by a release of water content from the sI structure as ice Ih, which indicates an increased cage occupancy. A similar lattice expansion was also described in another clathrate, SiO(2) clathrate, under high pressure. Such expansion with increasing cage occupancy might be a common manner to stabilize the clathrate structures under high pressure and low temperature.     

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 influence of pressure on the stability of clathrate hydrates of hydrogen and tetrahydrofuran

Changes in the Gibbs energy of hydration of molecular hydrogen and tetrahydrofuran (THF) at pressures of 0.1, 6.0, and 12.0 MPa over the temperature range 230–300 K were studied by the molecular dynamics method. The Gibbs energy of hydrogen in water-tetrahydrofuran-hydrogen solutions passed minima over the temperature range 235–265 K, which were indicative of a comparatively stable clathrate hydrate state. The Gibbs energy of the hydrogen molecule at the local minimum at 262 K was ∼4.5 kJ/mol; at atmospheric pressure and room temperature, it was ∼2 kJ/mol. An analysis of the radial distribution function and the coordination number of the THF molecule showed that, at 240–257 K, a clathrate hydrate of THF with the structure close to clathrate sII was predominantly formed.     

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.     

Rotation and diffusion of H2 in hydrogen-ice clathrate by 1H NMR

A recently reported hydrogen-ice clathrate carries up to four H(2) in each large cage and one H(2) in each small cage. We report pulsed proton NMR line shape measurements on H(2)-D(2)O clathrate formed at 1500 bar and 250 K. The behavior of the two-pulse spin-echo amplitude with respect to the nutation angle of the refocusing pulse shows that intramolecular dipolar broadening, modulated by H(2) molecular reorientations, dominates the line width of the ortho-H(2). Dipolar interaction between H(2) guests and host D atoms explains the echo variation with the relative phases of the pulses. From 12 to 120 K, the line width varies as 1/T, demonstrating that the three sublevels of J = 1 are split by a constant energy, epsilon. The splitting arises from distortion in the otherwise high-symmetry cages from frozen-out D(2)O orientational disorder. Above 120 K, further line-narrowing signals the onset of H(2) diffusion from cage to cage. At the lowest temperature, 1.9 K, the spectrum has Pake powder doublet-like features; the doublet is not fully developed, indicating a broad distribution of order parameters and energies epsilon.     

Clathrate-formation mediated adsorption of methane on Cu-complex crystals

We measured adsorption and desorption isotherms of methane on [Cu(4, 4'-bipyridine)2(BF4)2] (LPC) at 258, 273, and 303 K. Adsorption proceeds almost vertically at a definite pressure, which is named gate pressure. The lower the measurement temperature, the smaller the gate pressure. The temperature dependence of the gate pressure is expressed by the Clapeyron-Clausius equation, giving a thermodynamic evidence on the clathrate formation between the Cu complex and methane.     

The influence of gas phase fugacity and solubility on correlation of gas-hydrate formation pressure

The dissociation pressure for single gas-hydrate systems is correlated by van der Waals and Platteeuw's model with a Kihara spherical-core potential for the interaction between water and the guest molecule. By fitting to dissociation pressures along the hydrate-ice-gas-line, Kihara parameters are obtained independent of the mutual solubility of the gas and water. Further, the fugacity coefficients m that region are close to unity so that the Kihara parameters are rather insensitive to the choice of equation of state. By fitting to the ice-line only, we can investigate to what extent the equation of state fugacities influence the calculated dissociation pressures along the hydrate-water gas and hydrate-water-condensate equilibrium lines. By comparing the calculated equilibrium data to the experimental data we can conclude that an accurate prediction of hydrate dissociation pressures requires an equation of stat which gives an accurate correlation of fugacities rather than densities, and that it is necessary to take gas solubility into account even for non-polar gases like nitrogen and methane. Such investigation has not been reported previously. The Kihara parameters for the single gas-hydrate systems are estimated using the thermodynamically stable hydrate structure.     

Measurement and modeling of hydrocarbon dew points for five synthetic natural gas mixtures

The dew points of five synthetic natural gas (SNG) mixtures were measured using a custom made chilled mirror apparatus. The chilled mirror apparatus was designed to detect hydrocarbon dew points from low pressures up to the cricondenbar. The experimental temperature range was from 235 to 280 K and the pressure range from 0.3 to 10 MPa. The synthetic natural gases were comprised of methane and gravimetrically prepared fractions of ethane, propane, i-butane, n-butane and n-pentane. The experimental data were compared to calculations with the Soave–Redlich–Kwong (SRK) equation of state with classical mixing rule. However, considerable and increasing deviations between calculated and experimental dew points were observed as the pressure approached the cricondenbar. Therefore, a model was utilized based on the Redlich–Kwong (RK) equation of state. The Mathias and Copeman (MC) function was used to express the temperature dependence of the attractive term for all components. For methane, different sets of MC coefficients were used below and above the critical point. An optimization procedure was employed to fit the coefficients of supercritical methane to both pure component fugacity and experimental dew points. There is good agreement between experimental data and modeling results. For pressures higher than the pressure corresponding to the cricondentherm, the proposed model is better than the standard SRK equation of state. Good predictions with the model were obtained when comparing to bubble and dew points data from literature.     

High-pressure reactivity of clathrate hydrates by two-photon dissociation of water

It is shown that photoinduced reactions are observed at room temperature and pressure of few tenths of gigapascal in clathrate hydrates of CO and of model hydrocarbons under mild irradiation at 350 nm with power in the 50-610 mW range. The reactions are triggered by highly reactive OH radicals produced by two-photon excitation of the lowest electronic excited state of water having dissociative character. The formation of CO(2) is observed in all the reactions involving carbonaceous clathrate hydrates, and direct or indirect evidence for the formation of molecular hydrogen is obtained. The CO(2) produced in the reactions can be sequestered as a clathrate hydrate whose stability range seems to extend to room temperature at pressures of 0.5-0.6 GPa. Although the N(2) hydrate is stable up to 0.9 GPa under irradiation, a partial cleavage of the N-N triple bond is produced once the hydrate decomposes at 0.1 GPa.     

Free energies of carbon dioxide sequestration and methane recovery in clathrate hydrates

Classical molecular dynamics simulations are used to compare the stability of methane, carbon dioxide, nitrogen, and mixed CO(2)N(2) structure I (sI) clathrates under deep ocean seafloor temperature and pressure conditions (275 K and 30 MPa) which were considered suitable for CO(2) sequestration. Substitution of methane guests in both the small and large sI cages by CO(2) and N(2) fluids are considered separately to determine the separate contributions to the overall free energy of substitution. The structure I clathrate with methane in small cages and carbon dioxide in large cages is determined to be the most stable. Substitutions of methane in the small cages with CO(2) and N(2) have positive free energies. Substitution of methane with CO(2) in the large cages has a large negative free energy and substitution of the methane in the large cages with N(2) has a small positive free energy. The calculations show that under conditions where storage is being considered, carbon dioxide spontaneously replaces methane from sI clathrates, causing the release of methane. This process must be considered if there are methane clathrates present where CO(2) sequestration is to be attempted. The calculations also indicate that N(2) does not directly compete with CO(2) during methane substitution or clathrate formation and therefore can be used as a carrier gas or may be present as an impurity. Simulations further reveal that the replacement of methane with CO(2) in structure II (sII) cages also has a negative free energy. In cases where sII CO(2) clathrates are formed, only single occupancy of the large cages will be observed.