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.     

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.     

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.     

Molecular-dynamics study of structure II hydrogen clathrates

Molecular-dynamics simulations are used to study the stability of structure II hydrogen clathrates with different H2 guest occupancies. Simulations are done at pressures of 2.5 kbars and 1.013 bars and for temperatures ranging from 100 to 250 K. For a structure II unit cell with 136 water molecules, H2 guest molecule occupancies of 0-64 are studied with uniform occupancies among each type of cage. The simulations show that at 100 K and 2.5 kbars, the most stable configurations have single occupancy in the small cages and quadruple occupancy in the large cages. The optimum occupancy for the large cages decreases as the temperature is raised. Double occupancy in the small cages increases the energy of the structures and causes tetragonal distortion in the unit cell. The spatial distribution of the hydrogen guest molecules in the cages is determined by studying the guest-water and guest-guest radial distribution functions at various temperatures.     

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.     

Simulations of structure II H2 and D2 clathrates: potentials incorporating quantum corrections

Molecular dynamics simulations are used to study the stability of structure II H(2) and D(2) clathrates with different large and small guest occupancies at 160 and 250 K and 2.0 kbars. Simulations are performed with the recently proposed anisotropic site-site potentials of Wang for H2 and D2 [J. Quant. Spectrosc. Radiat. Transf. 76, 23 (2003)] which are parameterized to account for quantum corrections of order variant Planck's over 2pi(2) in the second virial coefficient. Occupancies of 0-2 in the small cages and 2-5 in the large cages are considered. Thermodynamic integration is used to determine the most stable guest occupancy at each temperature. Since lattice free energy and configurational energy differences are small for a number of different combinations of cage occupancies, one must expect that in bulk samples various combinations will indeed be observed. Special attention is given to the differences between H(2) and D(2) guests and implications on the hydrogen storage capacity of the clathrates are discussed.     

NMR shielding constants for hydrogen guest molecules in structure II clathrates

Proton NMR shielding constants and chemical shifts for hydrogen guests in small and large cages of structure II clathrates are calculated using density-functional theory and the gauge-invariant atomic-orbital method. Shielding constants are calculated at the B3LYP level with the 6-311++G(d,p) basis set. The calculated chemical shifts are corrected with a linear regression to reproduce the experimental chemical shifts of a set of standard molecules. The calculated chemical shifts of single hydrogen molecules in the small and large structure II cages are 4.94 and 4.84 ppm, respectively, which show that within the error range of the method the H2 guest molecules in the small and large cages cannot be distinguished. Chemical shifts are also calculated for double occupancy of the hydrogen guests in small cages, and double, triple, and quadruple occupancy in large cages. Multiple occupancy changes the chemical shift of the hydrogen guests by approximately 0.2 ppm. The relative effects of other guest molecules and the cage on the chemical shift are studied for the cages with multiple occupancies.     

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.     

Effect of small cage guests on hydrogen bonding of tetrahydrofuran in binary structure II clathrate hydrates

Molecular dynamics simulations of the pure structure II tetrahydrofuran clathrate hydrate and binary structure II tetrahydrofuran clathrate hydrate with CO(2), CH(4), H(2)S, and Xe small cage guests are performed to study the effect of the shape, size, and intermolecular forces of the small cages guests on the structure and dynamics of the hydrate. The simulations show that the number and nature of the guest in the small cage affects the probability of hydrogen bonding of the tetrahydrofuran guest with the large cage water molecules. The effect on hydrogen bonding of tetrahydrofuran occurs despite the fact that the guests in the small cage do not themselves form hydrogen bonds with water. These results indicate that nearest neighbour guest-guest interactions (mediated through the water lattice framework) can affect the clathrate structure and stability. The implications of these subtle small guest effects on clathrate hydrate stability are discussed.     

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.     

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.     

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.     

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.     

Molecular dynamics simulations for the growth of CH_4-CO_2 mixed hydrate

Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xCO2= 75%, xCO2= 50%, and xCO2= 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results show that the growth rate of CH4-CO2 mixed hydrate increases as the CO2 concentration in the initial solution phase increases and the temperature decreases. Via hydrate formation, the composition of CO2 in hydrate phase is higher than that in initial solution phase and the encaging capacity of CO2 in hydrates increases with the decrease in temperature. By analysis of the cage occupancy ratio of CH4 molecules and CO2 molecules in large cages to small cages, we find that CO2 molecules are preferably encaged into the large cages of the hydrate crystal as compared with CH4 molecules. Interestingly, CH4 molecules and CO2 molecules frequently replace with each other in some particular cage sites adjacent to hydrate/solution interface during the crystal growth process. These two species of guest molecules eventually act to stabilize the newly formed hydrates, with CO2 molecules occupying large cages and CH4 molecules occupying small cages in hydrate.     

Molecular dynamics simulations for the growth of CH4-CO2 mixed hydrate

Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xco2 = 75%, xco2 = 50%, and zco2 = 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results show that the growth rate of CH4-CO2 mixed hydrate increases as the CO2 concentration in the initial solution phase increases and the temperature decreases. Via hydrate formation, the composition of CO2 in hydrate phase is higher than that in initial solution phase and the encaging capacity of CO2 in hydrates increases with the decrease in temperature. By analysis of the cage occupancy ratio of CH4 molecules and CO2 molecules in large cages to small cages, we find that CO2 molecules are preferably encaged into the large cages of the hydrate crystal as compared with CH4 molecules. Interestingly, CH4 molecules and CO2 molecules frequently replace with each other in some particular cage sites adjacent to hydrate/solution interface during the crystal growth process. These two species of guest molecules eventually act to stabilize the newly formed hydrates, with CO2 molecules occupying large cages and CH4 molecules occupying small cages in hydrate.     

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.     

Theoretical Study of Clathrate Hydrates with Multiple Occupation

In this work, the electronic, structural, dynamic andthermodynamic properties of structure II, H and tetragonalAr clathrate hydrates have been calculated and the effectof multiple occupancy on their stability has been examinedusing first-principles and lattice dynamics calculations.The dynamic properties of these clathrates have beeninvestigated depending on the number of guest moleculesin a clathrate cage. It has been found that selectedhydrate structures are dynamically stable. The calculatedcell parameters are in agreement with experimental data.We also report the results of a systematic investigationof cage-like water structures using first-principles calculations. Ithas been observed that Ar clusters can be stabilized indifferent water cages and the stability is strongly dependenton the number of argon atoms inside the cages.     

Hydrogen encapsulation in a silicon clathrate type I structure: Na5.5(H2)2.15Si46: synthesis and characterization

A hydrogen-encapsulated inorganic clathrate, which is stable at ambient temperature and pressure, has been prepared in high yield. Na5.5(H2)2.15Si46 is a sodium-deficient, hydrogen-encapsulated, type I silicon clathrate. It was prepared by the reaction between NaSi and NH4Br under dynamic vacuum at 300 degrees C. The Rietveld refinement of the powder X-ray diffraction data is consistent with the clathrate type I structure. The type I clathrate structure has two types of cages where the guest species, in this case Na and H2, can reside: a large cage composed of 24 Si, in which the guest resides in the 6d crystallographic position, and a smaller one composed of 20 Si, in which the guest occupies the 2a position. Solid-state 23Na, 1H, and 29Si MAS NMR confirmed the presence of both sodium and hydrogen in the clathrate cages. 23Na NMR shows that sodium completely fills the small cage and is deficient in the larger cage. The 1H NMR spectrum shows a pattern consistent with mobile hydrogen in the large cage. 29Si NMR spectrum is consistent with phase pure type I clathrate framework. Elemental analysis is consistent with the stoichiometry Na5.5(H2.15)2Si46. The sodium occupancy was also examined using spherical aberration (Cs) corrected scanning transmission electron microscopy (STEM). The high-angle annular dark-field (HAADF) STEM experimental and simulated images indicated that the Na occupancy of the large cage, 6d sites, is less than 2/3, consistent with the NMR and elemental analysis.     

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.     

Molecular dynamics simulation of NMR powder lineshapes of linear guests in structure I clathrate hydrates

We perform molecular dynamics simulations (up to 6 ns) for the structure I clathrate hydrates of linear molecules CS, CS(2), OCS, and C(2)H(2) in large cages at different temperatures in the stability range to determine the angular distribution and dynamics of the guests in the large cages. The long axes of linear guest molecules in the oblate large structure I clathrate hydrate cages are primarily confined near the equatorial plane of the cage rather than axial regions. This non-uniform spatial distribution leads to well-known anisotropic lineshapes in the solid-state NMR spectra of the guest species. We use the dynamic distribution of guest orientations in the cages during the MD simulations at different temperatures to predict the (13)C NMR powder lineshapes of the guests in the large cages. The length of the guests and intermolecular interactions of the guests in the water cages determine the angular distribution and the mobility of the guests in the sI large cages at different temperatures. At low temperatures the range of motion of the guests in the cages are limited and this is reflected in the skew of the predicted (13)C lineshapes. As the guest molecules reach the fast motion limit at higher temperatures, the lineshapes for CS, OCS, and C(2)H(2) are predicted to have the "standard" powder lineshapes of guest molecules.