Molecular dynamics study of methane hydrate formation at a water/methane interface

We present molecular dynamics simulation results of a liquid water/methane interface, with and without an oligomer of poly(methylaminoethylmethacrylate), PMAEMA. PMAEMA is an active component of a commercial low dosage hydrate inhibitor (LDHI). Simulations were performed in the constant NPT ensemble at temperatures of 220, 235, 240, 245, and 250 K and a pressure of 300 bar. The simulations show the onset of methane hydrate growth within 30 ns for temperatures below 245 K in the methane/water systems; at 240 K there is an induction period of ca. 20 ns, but at lower temperatures growth commences immediately. The simulations were analyzed to calculate hydrate content, the propensity for hydrogen bond formation, and how these were affected by both temperature and the presence of the LDHI. As expected, both the hydrogen bond number and hydrate content decreased with increasing temperature, though little difference was observed between the lowest two temperatures considered. In the presence of PMAEMA, the temperature below which sustained hydrate growth occurred was observed to decrease. Some of the implications for the role of PMAEMA in LDHIs are discussed.     

Phase transformation of methane hydrate under high pressure

In situ Raman spectroscopy is employed to study the phase behavior of methane hydrate at high pressure. The structure 1 of methane hydrate can be maintained up to 950 MPa at 299 K. The transformation of structure I<-->structure H+water+CH4 occurs at 880 MPa and 323 K. The structure H of methane hydrate, however, decomposes to methane and water at 960 MPa and 348 K. The initiation mechanism of methane hydrate sI is also discussed.     

Ordering and Clathrate Hydrate Formation in Co‐deposits of Xenon and Water at Low Temperatures

Local ordering in co‐deposits of water and xenon atoms produced at low temperatures can be followed uniquely by ¹²⁹Xe NMR spectroscopy. In water‐rich samples deposited at 10 K and observed at 77 K, xenon NMR results show that there is a wide distribution of arrangements of water molecules around xenon atoms. This starts to order into the definite coordination for the structure I, large and small cages, when samples are annealed at ～140 K, although the process is not complete until a temperature of 180 K is reached, as shown by powder Xray diffraction. There is evidence that Xe ? 20 H₂O clusters are prominent in the early stages of crystallization. In xenon‐rich deposits at 77 K there is evidence of xenon atoms trapped in Xe ? 20 H₂O clusters, which are similar to the small hydration shells or cages observed in hydrate structures, but not in the larger water clusters consisting of 24 or 28 water molecules. These observations are in agreement with results obtained on the formation of Xe hydrate on the surface of ice surfaces by using hyperpolarized Xe NMR spectroscopy. The results indicate that for the various different modes of hydrate formation, both from Xe reacting with amorphous water and with crystalline ice surfaces, versions of the small cage are important structures in the early stages of crystallization.     

Molecular-dynamics simulations of methane hydrate dissociation

Nonequilibrium molecular-dynamics simulations have been carried out at 276.65 K and 68 bar for the dissolution of spherical methane hydrate crystallites surrounded by a liquid phase. The liquid was composed of pure water or a water-methane mixture ranging in methane composition from 50% to 100% of the corresponding theoretical maximum for the hydrate and ranged in size from about 1600 to 2200 water molecules. Four different crystallites ranging in size from 115 to 230 water molecules were used in the two-phase systems; the nanocrystals were either empty or had a methane occupation from 80% to 100% of the theoretical maximum. The crystal-liquid systems were prepared in two distinct ways, involving constrained melting of a bulk hydrate system or implantation of the crystallite into a separate liquid phase. The breakup rates were very similar for the four different crystal sizes investigated. The method of system preparation was not found to affect the eventual dissociation rates, despite a lag time of approximately 70 ps associated with relaxation of the liquid interfacial layer in the constrained melting approach. The dissolution rates were not affected substantially by methane occupation of the hydrate phase in the 80%-100% range. In contrast, empty hydrate clusters were found to break up significantly more quickly. Our simulations indicate that the diffusion of methane molecules to the surrounding liquid layer from the crystal surface appears to be the rate-controlling step in hydrate breakup. Increasing the size of the liquid phase was found to reduce the initial delay in breakup. We have compared breakup rates computed using different long-range electrostatic methods. Use of the Ewald, minimum image, and spherical cut-off techniques led to more rapid dissociation relative to the Lekner method.     

The structure of methane hydrate under geological conditions a combined Rietveld and maximum entropy analysis

We present a study of the structure of a fully deuterated methane hydrate under the geological conditions found in the world's oceans. In situ high-resolution neutron diffraction experiments have been performed at temperatures of 220, 275, and 280 K and a pressure of 100 bar, corresponding to the conditions at 1000 m water depth. The data were analyzed with a combination of Rietveld refinement and maximum entropy methods. From the Rietveld refinement, precise atomic parameters of the host lattice could be determined, indicating increasing distortions of the structure of the cages at elevated temperatures and pressures. Debye-Waller factors of the encaged CD(4) molecules have been found to exceed the values of the Debye-Waller factors of the D(2)O molecules considerably. In the large cage of structure type I the thermal center-of-mass displacements of the guests are 5-10 times larger than those of the water molecules. From the maximum entropy analysis maps of the scattering length density have been obtained, showing details of the vibrational amplitudes of the atoms in methane hydrate. The Debye-Waller factors of all molecules have been found to deviate considerably from a simple spherical geometry.     

Methane hydrate formation and decomposition: structural studies via neutron diffraction and empirical potential structure refinement

Neutron diffraction studies with hydrogen/deuterium isotope substitution measurements are performed to investigate the water structure at the early, medium, and late periods of methane clathrate hydrate formation and decomposition. These measurements are coupled with simultaneous gas consumption measurements to track the formation of methane hydrate from a gas/water mixture, and then the complete decomposition of hydrate. Empirical potential structure refinement computer simulations are used to analyze the neutron diffraction data and extract from the data the water structure in the bulk methane hydrate solution. The results highlight the significant changes in the water structure of the remaining liquid at various stages of hydrate formation and decomposition, and give further insight into the way in which hydrates form. The results also have important implications on the memory effect, suggesting that the water structure in the presence of hydrate crystallites is significantly different at equivalent stages of forming compared to decomposing. These results are in sharp contrast to the previously reported cases when all remaining hydrate crystallites are absent from the solution. For these systems there is no detectable change in the water structure or the methane hydration shell before hydrate formation and after decomposition. Based on the new results presented in this paper, it is clear that the local water structure is affected by the presence of hydrate crystallites, which may in turn be responsible for the "history" or "memory" effect where the production of hydrate from a solution of formed and then subsequently melted hydrate is reportedly much quicker than producing hydrate from a fresh water/gas mixture.     

Mechanisms for thermal conduction in hydrogen hydrate

Extensive equilibrium molecular dynamics simulations have been performed to investigate thermal conduction mechanisms via the Green-Kubo approach for (type II) hydrogen hydrate, at 0.05 kbar and between 30 and 250 K, for both lightly filled H(2) hydrates (1s4l) and for more densely filled H(2) systems (2s4l), in which four H(2) molecules are present in the large cavities, with respective single- and double-occupation of the small cages. The TIP4P water model was used in conjunction with a fully atomistic hydrogen potential along with long-range Ewald electrostatics. It was found that substantially less damping in guest-host energy transfer is present in hydrogen hydrate as is observed in common type I clathrates (e.g., methane hydrate), but more akin in to previous results for type II and H methane hydrate polymorphs. This gives rise to larger thermal conductivities relative to common type I hydrates, and also larger than type II and H methane hydrate polymorphs, and a more crystal-like temperature dependence of the thermal conductivity.     

A computer simulation of the mechanism of self-conservation of gas hydrates

Local density profiles and local component pressure profiles were obtained for two model systems containing methane hydrate and ice by molecular dynamics simulation. The ice matrix with methane hydrate clusters inserted into it was shown to be stable at normal pressure and even at a temperature higher than the temperature of methane hydrate dissociation. Calculations showed that the pressure in such a methane hydrate cluster inserted into ice was higher than in the ice phase. There were, however, no strong structure distortions because of the formation of a network of strong hydrogen bonds between the hydrate and ice phases.     

Preservation phenomena of methane hydrate in pore spaces

Dissociation processes of methane hydrate synthesized with glass beads were investigated using powder X-ray diffraction and calorimetry. Methane hydrate formed with coarse glass beads dissociated quickly at 150-200 K; in this temperature range methane hydrate dissociates at atmospheric pressure. In contrast, methane hydrate formed with glass beads less than a few microns in size showed very high stability up to just below the melting point of ice, even though this temperature is well outside the zone of thermodynamic stability of the hydrate. The rate-determining steps for methane hydrate dissociation within pores are also discussed. The experimental results suggest that methane hydrate existing naturally within the pores of fine particles such as mud at low temperatures would be significantly more stable than expected thermodynamically.     

Molecular dynamics simulation for surface melting and self-preservation effect of methane hydrate

The surface melting process of structure sI methane hydrate is simulated at T = 240, 260, 280, and 300 K using NVT molecular dynamics method. The simulation results show that a quasi-liquid layer will be formed during the melting process. The density distribution, translation, orientation, and dynamic properties of water molecules in the quasi-liquid layer are calculated as a function of the distance normal to the interface, which indicates the performance of quasi-liquid layer exhibits a continuous change from crystal-like to liquid-like. The quasi-liquid layer plays as a resistance of mass transfer restraining the diffusion of water and methane molecules during the melting process. The resistance of quasi-liquid layer will restrain methane molecules diffuse from hydrate phase to gas phase and slow the melting process, which can be considered as a possible mechanism of self-preservation effect. The performance of quasi-liquid layer is more crystal-like when the temperature is lower than the melt- ing-point of water, which will exhibit an obvious self-preservation. The self-preservation will weaken while the temperature is higher than the melting-point of water because of the liquid-like performance of the quasi-liquid layer.     

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.     

Gas hydrates of argon and methane synthesized at high pressures: composition, thermal expansion, and self-preservation

For the first time, the compositions of argon and methane high-pressure gas hydrates have been directly determined. The studied samples of the gas hydrates were prepared under high-pressure conditions and quenched at 77 K. The composition of the argon hydrate (structure H, stable at 460-770 MPa) was found to be Ar.(3.27 +/- 0.17)H(2)O. This result shows a good agreement with the refinement of the argon hydrate structure using neutron powder diffraction data and helps to rationalize the evolution of hydrate structures in the Ar-H(2)O system at high pressures. The quenched argon hydrate was found to dissociate in two steps. The first step (170-190 K) corresponds to a partial dissociation of the hydrate and the self-preservation of a residual part of the hydrate with an ice cover. Presumably, significant amounts of ice Ic form at this stage. The second step (210-230 K) corresponds to the dissociation of the residual part of the hydrate. The composition of the methane hydrate (cubic structure I, stable up to 620 MPa) was found to be CH(4).5.76H(2)O. Temperature dependence of the unit cell parameters for both hydrates has been also studied. Calculated from these results, the thermal expansivities for the structure H argon hydrate are alpha(a) = 76.6 K(-1) and alpha(c) = 77.4 K(-1) (in the 100-250 K temperature range) and for the cubic structure I methane hydrate are alpha(a) = 32.2 K(-1), alpha(a) = 53.0 K(-1), and alpha(a) = 73.5 K(-1) at 100, 150, and 200 K, respectively.     

水合物生成条件下甲烷在水中的吸附

采用悬滴法系统地测定了温度274.2 ~ 282.2 K、压力0.1 ~ 10.1 MPa下甲烷/纯水间界面张力。实验结果表明在恒定温度下界面张力随压力的增加而增大。在高压条件下,压力对界面张力有很大的影响。不同温度和压力下计算出的甲烷在水中的表面过剩浓度结果表明,压力越高,温度越低,甲烷在水溶液中的吸附浓度越高。同时,计算出的甲烷在水溶液中的表面吸附自由能结果表明,在水合物生成条件下,甲烷在水中的吸附比298.2 K更容易。     

Molecular Dynamics Simulation of Methane Hydrate Decomposition in the Presence of Alcohol Additives

The decomposition process of methane hydrate in pure water and methanol aqueous solution was studied by molecular dynamics simulation. The effects of temperature and pressure on hydrate structure and decomposition rate are discussed. The results show that decreasing pressure and increasing temperature can significantly enhance the decomposition rate of hydrate. After adding a small amount of methanol molecules, bubbles with a diameter of about 2 nm are formed, and the methanol molecules are mainly distributed at the gas-liquid interface, which greatly accelerates the decomposition rate and gas-liquid separation efficiency. The radial distribution function and sequence parameter analysis show that the water molecules of the undecomposed hydrate with ordered ice-like configuration at a temperature of 275 K evolve gradually into a long-range disordered liquid structure in the dynamic relaxation process. It was found that at temperatures above 280 K and pressures between 10 atm and 100 atm, the pressure has no significant effect on hydrate decomposition rate, but when the pressure is reduced to 1 atm, the decomposition rate increases sharply. These findings provided a theoretical insight for the industrial exploitation of hydrates.     

Molecular dynamics simulations of the mechanisms of thermal conduction in methane hydrates

The thermal conductivity of methane hydrate is an important physical parameter affecting the processes of methane hydrate exploration,mining,gas hydrate storage and transportation as well as other applications.Equilibrium molecular dynamics simulations and the Green-Kubo method have been employed for systems from fully occupied to vacant occupied sI methane hydrate in order to estimate their thermal conductivity.The estimations were carried out at temperatures from 203.15 to 263.15 K and at pressures from 3 to 100 MPa.Potential models selected for water were TIP4P,TIP4P-Ew,TIP4P/2005,TIP4P-FQ and TIP4P/Ice.The effects of varying the ratio of the host and guest molecules and the external thermobaric conditions on the thermal conductivity of methane hydrate were studied.The results indicated that the thermal conductivity of methane hydrate is essentially determined by the cage framework which constitutes the hydrate lattice and the cage framework has only slightly higher thermal conductivity in the presence of the guest molecules.Inclusion of more guest molecules in the cage improves the thermal conductivity of methane hydrate.It is also revealed that the thermal conductivity of the sI hydrate shows a similar variation with temperature.Pressure also has an effect on the thermal conductivity,particularly at higher pressures.As the pressure increases,slightly higher thermal conductivities result.Changes in density have little impact on the thermal conductivity of methane hydrate.     

Antifreezes Act as Catalysts for Methane Hydrate Formation from Ice

Contrary to the thermodynamic inhibiting effect of methanol on methane hydrate formation from aqueous phases, hydrate forms quickly at high yield by exposing frozen water–methanol mixtures with methanol concentrations ranging from 0.6–10 wt % to methane gas at pressures from 125 bars at 253 K. Formation rates are some two orders of magnitude greater than those obtained for samples without methanol and conversion of ice is essentially complete. Ammonia has a similar catalytic effect when used in concentrations of 0.3–2.7 wt %. The structure I methane hydrate formed in this manner was characterized by powder X‐ray diffraction and Raman spectroscopy. Steps in the possible mechanism of action of methanol were studied with molecular dynamics simulations of the Ih (0001) basal plane exposed to methanol and methane gas. Simulations show that methanol from a surface aqueous layer slowly migrates into the ice lattice. Methane gas is preferentially adsorbed into the aqueous methanol surface layer. Possible consequences of the catalytic methane hydrate formation on hydrate plug formation in gas pipelines, on large scale energy‐efficient gas hydrate formation, and in planetary science are discussed.     

Lifetimes of cagelike water clusters immersed in bulk liquid water: a molecular dynamics study on gas hydrate nucleation mechanisms

Molecular dynamics simulations were performed to observe the evolution of cagelike water clusters immersed in bulk liquid water at 250 and 230 K. Totally, we considered four types of clusters--dodecahedral (5(12)) and tetrakaidecahedral (5(12)6(2)) cagelike water clusters filled with or without a methane molecule, respectively. The lifetimes of these clusters were calculated according to their Lindemann index (delta) using the criterion of delta> or =0.07. The lifetimes of the clusters at 230 K are longer than that at 250 K, and their ratios are the same as the ratio of structure relaxation times of bulk water at these temperatures. For both the filled and empty clusters, the lifetimes of 5(12)6(2) cagelike clusters are similar to that of 5(12) cagelike clusters. Although the methane molecules indeed make the filled cagelike water clusters live longer than the empty ones, the empty cagelike water clusters still have the chance of being long lived. These observations support the cluster nucleation hypothesis for the formation mechanisms of gas hydrates.     

Methane hydrate stability in the presence of water-soluble hydroxyalkyl cellulose

The effect of low-dosage water-soluble hydroxyethyl cellulose (approximate MW～90,000 and 250,000) as a member of hydroxyalkyl cellulosic polymer group on methane hydrate stability was investigated by monitoring hydrate dissociation at pressures greater than atmospheric pressure in a closed vessel. In particular, the influence of molecular weight and mass concentration of hydroxyethyl cellulose (HEC) was studied with respect to hydrate formation and dissociation. Methane hydrate formation was performed at 2℃ and at a pressure greater than 100 bar. Afterwards, hydrate dissociation was initiated by step heating from -10℃ at a mild pressure of 13 bar to 3℃, 0℃ and 2℃. With respect to the results obtained for methane hydrate formation/dissociation and the amount of gas uptake, we concluded that HEC 90,000 at 5000 ppm is suitable for long-term gas storage and transportation under a mild pressure of 13 bar and at temperatures below the freezing point.     

Application of the cell potential method to predict phase equilibria of multicomponent gas hydrate systems

We present the application of a mathematical method reported earlier by which the van der Waals-Platteeuw statistical mechanical model with the Lennard-Jones and Devonshire approximation can be posed as an integral equation with the unknown function being the intermolecular potential between the guest molecules and the host molecules. This method allows us to solve for the potential directly for hydrates for which the Langmuir constants are computed, either from experimental data or from ab initio data. Given the assumptions made in the van der Waals-Platteeuw model with the spherical-cell approximation, there are an infinite number of solutions; however, the only solution without cusps is a unique central-well solution in which the potential is at a finite minimum at the center to the cage. From this central-well solution, we have found the potential well depths and volumes of negative energy for 16 single-component hydrate systems: ethane (C2H6), cyclopropane (C3H6), methane (CH4), argon (Ar), and chlorodifluoromethane (R-22) in structure I; and ethane (C2H6), cyclopropane (C3H6), propane (C3H8), isobutane (C4H10), methane (CH4), argon (Ar), trichlorofluoromethane (R-11), dichlorodifluoromethane (R-12), bromotrifluoromethane (R-13B1), chloroform (CHCl3), and 1,1,1,2-tetrafluoroethane (R-134a) in structure II. This method and the calculated cell potentials were validated by predicting existing mixed hydrate phase equilibrium data without any fitting parameters and calculating mixture phase diagrams for methane, ethane, isobutane, and cyclopropane mixtures. Several structural transitions that have been determined experimentally as well as some structural transitions that have not been examined experimentally were also predicted. In the methane-cyclopropane hydrate system, a structural transition from structure I to structure II and back to structure I is predicted to occur outside of the known structure II range for the cyclopropane hydrate. Quintuple (L(w)-sI-sII-L(hc)-V) points have been predicted for the ethane-propane-water (277.3 K, 12.28 bar, and x(eth,waterfree) = 0.676) and ethane-isobutane-water (274.7 K, 7.18 bar, and x(eth,waterfree) = 0.81) systems.