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.     

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.     

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

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.     

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.     

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.     

Guest-Host Interaction Study in Clathrate Hydrates Using Lattice Dynamics Simulation

Lattice dynamics simulation of several gas hydrates （helium, argon, and methane） with different occupancy rates has been performed using TIP3P potential model. Results show that the coupling between the guest and host is not simple as depicted by the conventional viewpoints. For clathrate hydrate enclosing small guest, the small cages are dominantly responsible for the thermodynamic stability of clathrate hydrates. And the spectrum of methane hydrate is studied compared with argon hydrate, then as a result, shrink effect from positive hydrogen shell is proposed.     

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.     

Hydrogen storage in sH hydrates: a Monte Carlo study

Grand canonical Monte Carlo simulations are performed to evaluate the hydrogen-storage capacity of the recently discovered hydrogen hydrates of the sH type, at 274 K and up to 500 MPa. First, the pure H2 hydrate is investigated in order to determine the upper limit of H 2 content in sH hydrates. It is found that the storage capacity of the hypothetical pure H2 hydrate could reach 3.6 wt % at 500 MPa. Depending on pressure, the large cavity of this hydrate can accommodate up to eight H2 molecules, while the small and medium ones are singly occupied even at pressures as high as 500 MPa. Next, the binary H2-methylcyclohexane sH hydrate is examined. In this case, the small and medium cavities are again singly occupied, resulting in a maximum H2 uptake of 1.4 wt %. Finally, the results from simulations on pure H2 and binary hydrates are utilized to investigate the potential of H2 storage in sH hydrates where the promoter molecules occupy the medium instead of the large cavities.     

Phase equilibrium measurements and the tuning behavior of new sII clathrate hydrates

We suggest two types of new amine-type sII formers: pyrrolidine and piperidine. These guest compounds fail to form clathrate hydrate structures with host water, but instead have to combine with light gaseous guest molecules (methane) for enclathration. First, two binary clathrate hydrates of (pyrrolidine + methane) and (piperidine + methane) were synthesized at various amine concentrations. ¹³C NMR and Raman analysis were done to identify the clathrate hydrate structure and guest distribution over sII-S and sII-L cages. XRD was also used to find the exact structure and corresponding cell parameters. At a dilute pyrrolidine concentration of less than 5.56 mol%, the tuning phenomenon is observed such that methane molecules surprisingly occupy sII-L cages. At the critical guest concentration of about 0.1 mol%, the cage occupancy ratio reaches the maximum of approximately 0.5. At very dilute guest concentration below 0.1 mol%, the methane molecules fail to occupy large cages on account of their rarefied distribution in the network. Direct-release experiments were performed to determine the actual guest compositions in the clathrate hydrate phases. Finally, we measured the clathrate hydrate phase equilibria of (pyrrolidine + methane) and (piperidine + methane).     

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.     

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.     

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.     

The cages, dynamics, and structuring of incipient methane clathrate hydrates

Interest in describing clathrate hydrate formation mechanisms spans multiple fields of science and technical applications. Here, we report findings from multiple molecular dynamics simulations of spontaneous methane clathrate hydrate nucleation and growth from fully demixed and disordered two-phase fluid systems of methane and water. Across a range of thermodynamic conditions and simulation geometries and sizes, a set of seven cage types comprises approximately 95% of all cages formed in the nucleated solids. This set includes the ubiquitous 5(12) cage, the 5(12)6(n) subset (where n ranges from 2-4), and the 4(1)5(10)6(n) subset (where n also ranges from 2-4). Transformations among these cages occur via water pair insertions/removals and rotations, and may elucidate the mechanisms of solid-solid structural rearrangements observed experimentally. Some consistency is observed in the relative abundance of cages among all nucleation trajectories. 5(12) cages are always among the two most abundant cage types in the nucleated solids and are usually the most abundant cage type. In all simulations, the 5(12)6(n) cages outnumber their 4(1)5(10)6(n) counterparts with the same number of water molecules. Within these consistent features, some stochasticity is observed in certain cage ratios and in the long-range ordering of the nucleated solids. Even when comparing simulations performed at the same conditions, some trajectories yield swaths of multiple adjacent sI unit cells and long-range order over 5 nm, while others yield only isolated sI unit cells and little long-range order. The nucleated solids containing long-range order have higher 5(12)6(2)/5(12) and 5(12)6(3)/4(1)5(10)6(2) cage ratios when compared to systems that nucleate with little long-range order. The formation of multiple adjacent unit cells of sI hydrate at high driving forces suggests an alternative or addition to the prevailing hydrate nucleation hypotheses which involve formation through amorphous intermediates.     

Molecular dynamics study of carbon dioxide hydrate dissociation

We present results from a molecular dynamics study of the dissociation behavior of carbon dioxide (CO(2)) hydrates. We explore the effects of hydrate occupancy and temperature on the rate of hydrate dissociation. We quantify the rate of dissociation by tracking CO(2) release into the liquid water phase as well as the velocity of the hydrate-liquid water interface. Our results show that the rate of dissociation is dependent on the fractional occupancy of each cage type and cannot be described simply in terms of overall hydrate occupancy. Specifically, we find that hydrates with similar overall occupancy differ in their dissociation behavior depending on whether the small or large cages are empty. In addition, individual cages behave differently depending on their surrounding environment. For the same overall occupancy, filled small and large cages dissociate faster in the presence of empty large cages than when empty small cages are present. Therefore, hydrate dissociation is a collective phenomenon that cannot be described by focusing solely on individual cage behavior.     

The effect of the water/methane interface on methane hydrate cages: the potential of mean force and cage lifetimes

Molecular dynamics simulations were used to determine the influence of a methane-water interface on the position and stability of methane hydrate cages. A potential of mean force was calculated as a function of the separation of a methane hydrate cage and a methane-water interface. The hydrate cages are found to be strongly repelled from the methane gas into the water phase. At low enough temperatures, however, the most favorable location for the hydrate cage is at the interface on the water side. Cage lifetime simulations were performed in bulk water and near a methane-water interface. The methane-water interface increases the cage lifetime by almost a factor of 2 compared to cage lifetimes of cages in bulk water. The potential of mean force and the cage lifetime results give additional explanations for the proposed nucleation of gas hydrates at gas-water interfaces.     

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.     

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.