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 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.     

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.     

Molecular-dynamics simulations of binary structure II hydrogen and tetrahydrofurane clathrates

The binary structure II hydrogen and tetrahydrofurane (THF) clathrates are studied with molecular-dynamics simulations. Simulations are done at pressures of 120 and 1.013 bars for temperatures ranging from 100 to 273 K. For the small cages of the structure II unit cell, H2 guest molecule occupancies of 0, 16 (single occupancy), and 32 (double occupancy) are considered. THF occupancies of 0-8 in the large cages are studied. For cases in which THF does not occupy all large cages in a unit cell, the remaining large cages can be occupied with sets of four H2 guest molecules. The unit-cell volumes and configurational energies are compared in the different occupancy cases. Increasing the small cage occupancy leads to an increase in the unit-cell volume and thermal-expansion coefficient. Among simulations with the same small cage occupancy, those with the large cages containing 4H2 guests have the largest volumes. The THF guest molecules have a stabilizing effect on the clathrate and the configurational energy of the unit cell decreases linearly as the THF content increases. For binary THF + H2 clathrates, the substitution of the THF molecules in the large cages with sets of 4H2 molecules increases the configurational energy. For the binary clathrates, various combinations of THF and H2 occupancies have similar configurational energies.     

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.     

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.     

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.     

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.     

Determination of NMR Lineshape Anisotropy of Guest Molecules within Inclusion Complexes from Molecular Dynamics Simulations

Nonspherical cages in inclusion compounds can result in non‐uniform motion of guest species in these cages and anisotropic lineshapes in NMR spectra of the guest. Herein, we develop a methodology to calculate lineshape anisotropy of guest species in cages based on molecular dynamics simulations of the inclusion compound. The methodology is valid for guest atoms with spin 1/2 nuclei and does not depend on the temperature and type of inclusion compound or guest species studied. As an example, the nonspherical shape of the structure I (sI) clathrate hydrate large cages leads to preferential alignment of linear CO₂ molecules in directions parallel to the two hexagonal faces of the cages. The angular distribution of the CO₂ guests in terms of a polar angle θ and azimuth angle ? and small amplitude vibrational motions in the large cage are characterized by molecular dynamics simulations at different temperatures in the stability range of the CO₂ sI clathrate. The experimental ¹³C NMR lineshapes of CO₂ guests in the large cages show a reversal of the skew between the low temperature (77 K) and the high temperature (238 K) limits of the stability of the clathrate. We determine the angular distributions of the guests in the cages by classical MD simulations of the sI clathrate and calculate the ¹³C NMR lineshapes over a range of temperatures. Good agreement between experimental lineshapes and calculated lineshapes is obtained. No assumptions regarding the nature of the guest motions in the cages are required.     

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.     

Thermal expansivity of tetrahydrofuran clathrate hydrate with diatomic guest molecules

The guest dynamics and thermal behavior occurring in the cages of clathrate hydrates appear to be too complex to be clearly understood through various structural and spectroscopic approaches, even for the well-known structures of sI, sII, and sH. Neutron diffraction studies have recently been carried out to clarify the special role of guests in expanding the host water lattices and have contributed to revealing the influence factors on thermal expansivity. Through this letter we attempt to address three noteworthy features occurring in guest inclusion: (1) the effect of guest dimension on host water lattice expansion; (2) the effect of thermal history on host water lattice expansion; and (3) the effect of coherent/incoherent scattering cross sections on guest thermal patterns. The diatomic guests of H 2, D 2, N 2, and O 2 have been selected for study, and their size and mass dependence on the degree of lattice expansion have been examined, and four sII clathrate hydrates with tetrahydrofuran (THF) have been synthesized in order to determine their neutron powder diffraction patterns. After thermal cycling, the THF + H 2 clathrate hydrate is observed to exhibit an irreversible plastic deformation-like pattern, implying that the expanded lattices fail to recover the original state by contraction. The host-water cage dimension after degassing the guest molecules remains as it was expanded, and thus host-guest as well as guest-guest interactions will be altered if guest uptake reoccurs.     

13C NMR studies of hydrocarbon guests in synthetic structure H gas hydrates: experiment and computation

(13)C NMR chemical shifts were measured for pure (neat) liquids and synthetic binary hydrate samples (with methane help gas) for 2-methylbutane, 2,2-dimethylbutane, 2,3-dimethylbutane, 2-methylpentane, 3-methylpentane, methylcyclopentane, and methylcyclohexane and ternary structure H (sH) clathrate hydrates of n-pentane and n-hexane with methane and 2,2-dimethylbutane, all of which form sH hydrates. The (13)C chemical shifts of the guest atoms in the hydrate are different from those in the free form, with some carbon atoms shifting specifically upfield. Such changes can be attributed to conformational changes upon fitting the large guest molecules in hydrate cages and/or interactions between the guests and the water molecules of the hydrate cages. In addition, powder X-ray diffraction measurements revealed that for the hexagonal unit cell, the lattice parameter along the a-axis changes with guest hydrate former molecule size and shape (in the range of 0.1 ?) but a much smaller change in the c-axis (in the range of 0.01 ?) is observed. The (13)C NMR chemical shifts for the pure hydrocarbons and all conformers were calculated using the gauge invariant atomic orbital method at the MP2/6-311+G(2d,p) level of theory to quantify the variation of the chemical shifts with the dihedral angles of the guest molecules. Calculated and measured chemical shifts are compared to determine the relative contribution of changes in the conformation and guest-water interactions to the change in chemical shift of the guest upon clathrate hydrate formation. Understanding factors that affect experimental chemical shifts for the enclathrated hydrocarbons will help in assigning spectra for complex hydrates recovered from natural sites.     

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.     

Effect of Guest–Host Hydrogen Bonding on the Structures and Properties of Clathrate Hydrates

To provide improved understanding of guest–host interactions in clathrate hydrates, we present some correlations between guest chemical structures and observations on the corresponding hydrate properties. From these correlations it is clear that directional interactions such as hydrogen bonding between guest and host are likely, although these have been ignored to greater or lesser degrees because there has been no direct structural evidence for such interactions. For the first time, single‐crystal X‐ray crystallography has been used to detect guest–host hydrogen bonding in structure II (sII) and structure H (sH) clathrate hydrates. The clathrates studied are the tert‐butylamine (tBA) sII clathrate with H₂S/Xe help gases and the pinacolone + H₂S binary sH clathrate. X‐ray structural analysis shows that the tBA nitrogen atom lies at a distance of 2.64 Å from the closest clathrate hydrate water oxygen atom, whereas the pinacolone oxygen atom is determined to lie at a distance of 2.96 Å from the closest water oxygen atom. These distances are compatible with guest–water hydrogen bonding. Results of molecular dynamics simulations on these systems are consistent with the X‐ray crystallographic observations. The tBA guest shows long‐lived guest–host hydrogen bonding with the nitrogen atom tethered to a water HO group that rotates towards the cage center to face the guest nitrogen atom. Pinacolone forms thermally activated guest–host hydrogen bonds with the lattice water molecules; these have been studied for temperatures in the range of 100–250 K. Guest–host hydrogen bonding leads to the formation of Bjerrum L‐defects in the clathrate water lattice between two adjacent water molecules, and these are implicated in the stabilities of the hydrate lattices, the water dynamics, and the dielectric properties. The reported stable hydrogen‐bonded guest–host structures also tend to blur the longstanding distinction between true clathrates and semiclathrates.     

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.     

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.     

Intermolecular hydrogen transfer between guest species in small and large cages of methane + propane mixed gas hydrates

To investigate the molecular interaction between guest species inside of the small and large cages of methane + propane mixed gas hydrates, thermal stabilities of the methyl radical (possibly induced in small cages) and the normal propyl and isopropyl radicals (induced in large cages) were investigated by means of electron spin resonance measurements. The increase of the total amount of the normal propyl and isopropyl radicals reveals that the methyl radical in the small cage withdraws one hydrogen atom from the propane molecule enclathrated in the adjacent large cage of the structure-II hydrate. A guest species in a hydrate cage has the ability to interact closely with the other one in the adjacent cages. The clathrate hydrate may be utilized as a possible nanoscale reaction field.     

13C chemical shifts of propane molecules encaged in structure II clathrate hydrate

Experimental NMR measurements for (13)C chemical shifts of propane molecules encaged in 16-hedral cages of structure II clathrate hydrate were conducted to investigate the effects of guest-host interaction of pure propane clathrate on the (13)C chemical shifts of propane guests. Experimental (13)C NMR measurements revealed that the clathrate hydration of propane reverses the (13)C chemical shifts of methyl and methylene carbons in propane guests to gaseous propane at room temperature and atmospheric pressure or isolated propane, suggesting a change in magnetic environment around the propane guest by the clathrate hydration. Inversion of the (13)C chemical shifts of propane clathrate suggests that the deshielding effect of the water cage on the methyl carbons of the propane molecule encaged in the 16-hedral cage is greater than that on its methylene carbon.     

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.     

The thermodynamic stability of clathrate hydrate. Encaging non-spherical propane molecules

The thermodynamic stability of a clathrate hydrate encaging non-spherical molecules has been investigated by examining the free energy of cage occupancy. In the present study, a generalized van der Waals and Platteeuw theory is extended to treat the rotational motion of guest molecules in clathrate hydrate cages. The vibrational free energy of both guest and host molecules is divided into harmonic and anharmonic contributions. The anharmonic free energy associated with the non-spherical nature of the guest molecules is evaluated as a perturbation from the spherical guest. Predicted thermodynamic properties are compared with measured values. It is shown that this anharmonic contribution is important in the free energy of the hindered rotation of the guests.