Heterogeneous crystal growth of methane hydrate on its sII [001] crystallographic face

This paper presents a systematic molecular simulation study of the heterogeneous crystal growth of methane hydrate sII from supersaturated aqueous methane solutions. The growth of sII hydrate on the [001] crystallographic face is achieved through utilization of a recently proposed methodology, and rates of crystal growth of 1 A/ns were sustained for the molecular models and specific conditions employed in this work. Characteristics of the crystals grown as well as properties and structure of the interface are examined. Water cages with a 5(12)6(3) arrangement, which are improper to both sI and sII structures, are identified during the heterogeneous growth of sII methane hydrate. We show that the growth of a [001] face of sII hydrate can produce an sI crystalline structure, confirming that cross-nucleation of methane hydrate structures is possible. Defects consisting of two methane molecules trapped in large 5(12)6(4) cages and water molecules trapped in small and large cages are observed, where in one instance we have found a large 5(12)6(4) cage containing three water molecules.     

Molecular dynamics methodology to investigate steady-state heterogeneous crystal growth

In this paper a new molecular dynamics simulation methodology to investigate steady-state heterogeneous crystal growth from a supercooled liquid is presented. The method is tested on pure component systems such as Lennard-Jonesium and water/ice, as well as multicomponent systems such as methane hydrate crystals. The setup uses periodicity in all three directions and two interfaces; at one interface, crystallization occurs, while at the other, melting is enforced by locally heating the crystal only near that interface. Steady-state conditions are achieved when the crystal is melted at the same rate as the growth occurs. A self-adaptive scheme that automatically modifies the rate of melting to match the rate of growth, crucial for establishing steady-state conditions, is described. In contrast with the recently developed method of Razul et al. [Mol. Phys. 103, 1929 (2005)], where the rates of growth (melting) were constant and the temperatures determined, the present approach fixes the supercooling temperature at the growing interface and identifies the corresponding steady-state crystal growth rate that corresponds to the thermodynamic force provided. The static properties of the interface (e.g., the interfacial widths) and the kinetics of the crystal growth are found to reproduce well previous findings. The importance of establishing steady-state conditions in such investigations is also briefly discussed.     

Crystal growth simulations of methane hydrates in the presence of silica surfaces

We present a molecular dynamics simulation study of the crystal growth of methane hydrates in the presence of model silica (SiO(2)) surfaces. The crystal growth under apparent steady-state conditions shows a clear preference for bulk solution. We observe rather disordered water arrangements very close to the silica surface within about 5 ? in both liquid and crystalline regions of the system. These disordered structures have dynamic and structural properties intermediate between those exhibited by molecules in bulk liquid and crystalline phases. The presence of methane molecules appears to help stabilize these structures. We observe that under appropriate conditions, the hydroxylated silica surfaces can serve as a source of methane molecules which can help promote hydrate growth near the surfaces.     

Analysis of Decomposition for Structure I Methane Hydrate by Molecular Dynamics Simulation

Under multi-nodes of temperatures and pressures, microscopic decomposition mechanisms of structure I methane hydrate in contact with bulk water molecules have been studied through LAMMPS software by molecular dynamics simulation. Simulation system consists of 482 methane molecules in hydrate and 3027 randomly distributed bulk water molecules. Through analyses of simulation results, decomposition number of hydrate cages, density of methane molecules, radial distribution function for oxygen atoms, mean square displacement and coefficient of diffusion of methane molecules have been studied. A significant result shows that structure I methane hydrate decomposes from hydrate-bulk water interface to hydrate interior. As temperature rises and pressure drops, the stabilization of hydrate will weaken, decomposition extent will go deep, and mean square displacement and coefficient of diffusion of methane molecules will increase. The studies can provide important meanings for the microscopic decomposition mechanisms analyses of methane hydrate.     

Gas hydrate nucleation and cage formation at a water/methane interface

Nucleation of gas hydrates remains a poorly understood phenomenon, despite its importance as a critical step in understanding the performance and mode of action of low dosage hydrate inhibitors. We present here a detailed analysis of the structural and mechanistic processes by which gas hydrates nucleate in a molecular dynamics simulation of dissolved methane at a methane/water interface. It was found that hydrate initially nucleates into a phase consistent with none of the common bulk crystal structures, but containing structural units of all of them. The process of water cage formation has been found to correlate strongly with the collective arrangement of methane molecules.     

甲烷水合物导热机理的分子动力学模拟

甲烷水合物导热系数是甲烷水合物勘探、开采、储运以及其他应用过程中一个十分重要的物理参数.我们采用平衡分子动力学（EMD）方法Green-Kubo理论计算温度203.15～263.15K、压力范围3～100MPa、晶穴占有率为0～1的sI甲烷水合物的导热系数,采用的水分子模型包括TIP4P、TIP4P-Ew、TIP4P-FQ、TIP4P/2005、TIP4P/Ice.研究了主客体分子、外界温压条件等对甲烷水合物导热性能的影响.研究结果显示甲烷水合物的低导热性能由主体分子构建的sI笼型结构决定,而客体分子进入笼型结构后,使得笼型结构导热性能增强,同时进入笼型结构的客体分子越多,甲烷水合物导热性能越强.研究结果还显示在高温区域（T〉TDebye/3）内不同温度作用下,所有sI水合物具有相似的导热规律.压力对导热系数有一定影响,尤其是在较高压力条件下,压力越高,导热系数越大.而在不同温度和不同压力作用过程中,密度的改变对导热系数的增大或减小几乎没有影响.     

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.     

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.     

Molecular simulation of non-equilibrium methane hydrate decomposition process

We recently performed constant energy molecular dynamics simulations of the endothermic decomposition of methane hydrate in contact with water to study phenomenologically the role of mass and heat transfer in the decomposition rate [S. Alavi, J.A. Ripmeester, J. Chem. Phys. 132 (2010) 144703]. We observed that with the progress of the decomposition front temperature gradients are established between the remaining solid hydrate and the solution phases. In this work, we provide further quantitative macroscopic and molecular level analysis of the methane hydrate decomposition process with an emphasis on elucidating microscopic details and how they affect the predicted rate of methane hydrate decomposition in natural methane hydrate reservoirs. A quantitative criterion is used to characterize the decomposition of the hydrate phase at different times. Hydrate dissociation occurs in a stepwise fashion with rows of sI cages parallel to the interface decomposing simultaneously. The correlations between decomposition times of subsequent layers of the hydrate phase are discussed.     

Calculation of free energies and chemical potentials for gas hydrates using Monte Carlo simulations

We describe a method for calculating free energies and chemical potentials for molecular models of gas hydrate systems using Monte Carlo simulations. The method has two components: (i) thermodynamic integration to obtain the water and guest molecule chemical potentials as functions of the hydrate occupancy; (ii) calculation of the free energy of the zero-occupancy hydrate system using thermodynamic integration from an Einstein crystal reference state. The approach is applicable to any classical molecular model of a hydrate. We illustrate the methodology with an application to the structure-I methane hydrate using two molecular models. Results from the method are also used to assess approximations in the van der Waals-Platteeuw theory and some of its extensions. It is shown that the success of the van der Waals-Platteeuw theory is in part due to a cancellation of the error arising from the assumption of a fixed configuration of water molecules in the hydrate framework with that arising from the neglect of methane-methane interactions.     

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.     

常见客体分子对笼型水合物晶格常数的影响

Natural gas hydrates are considered as ideal alternative energy resources for the future, and the relevant basic and applied research has become more attractive in recent years. The influence of guest molecules on the hydrate crystal lattice parameters is of great significances to the understanding of hydrate structural characteristics, hydrate formation/decomposition mechanisms, and phase stability behaviors. In this study, we test a series of artificial hydrate samples containing different guest molecules (e.g. methane, ethane, propane, iso-butane, carbon dioxide, tetrahydrofuran, methane + 2, 2-dimethylbutane, and methane + methyl cyclohexane) by a low-temperature powder X-ray diffraction (PXRD). Results show that PXRD effectively elucidates structural characteristics of the natural gas hydrate samples, including crystal lattice parameters and structure types. The relationships between guest molecule sizes and crystal lattice parameters reveal that different guest molecules have different controlling behaviors on the hydrate types and crystal lattice constants. First, a positive correlation between the lattice constants and the van der Waals diameters of homologous hydrocarbon gases was observed in the single-guest-component hydrates. Small hydrocarbon homologous gases, such as methane and ethane, tended to form sI hydrates, whereas relatively larger molecules, such as propane and iso-butane, generated sⅡ hydrates. The hydrate crystal lattice constants increased with increasing guest molecule size. The types of hydrates composed of oxygen-containing guest molecules (such as CO₂ and THF) were also controlled by the van der Waals diameters. However, no positive correlation between the lattice constants and the van der Waals diameters of guest molecules in hydrocarbon hydrates was observed for CO₂ hydrate and THF hydrate, probably due to the special interactions between the guest oxygen atoms and hydrate "cages". Furthermore, the influences of the macromolecules and auxiliary small molecules on the lengths of the different crystal axes of the sH hydrates showed inverse trends. Compared to the methane + 2, 2-dimethylbutane hydrate sample, the length of the a-axis direction of the methane + methyl cyclohexane hydrate sample was slightly smaller, whereas the length of the c-axis direction was slightly longer. The crystal a-axis length of the sH hydrate sample formed with nitrogen molecules was slightly longer, whereas the c-axis was shorter than that of the methane + 2, 2-dimethylbutane hydrate sample at the same temperature.     

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.     

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.     

An atomistic level description of guest molecule effect on the formation of hydrate crystal nuclei by <Emphasis Type="Italic">ab initio</Emphasis> calculations

In the present study, we report the results of a systematic investigation of cage-like water structures using the first-principles calculations. These results show that, in the case of methane hydrate, the following nucleation mechanism can be revealed. The formation of small water cavities filled with methane is the first step of the formation of methane hydrate. It is not necessary to occupy all dodecahedral cages by guest molecules. After that small cavities start to form the H-bonding network with surrounding water molecules and a small number of water molecules is enough for the formation of a stable hydrogen-bonding network. The structural information contained in such nuclei is conserved in the forming crystal. Moreover, the presence of a methane molecule between small cages is also important to prevent the adhesion of cavities. It found that the ozone molecule can also stabilize the small cage since the value of the interaction energy between the ozone guest and the water host framework is very close to that obtained for the methane case. However, ozone affects the structure of large cavities and hence, the second guest is necessary to stabilize the hydrate structure.     

Molecular Insights into the Effect of Nitrogen Bubbles on the Formation of Tetrahydrofuran Hydrates

In this work, a molecular dynamics simulation was conducted to study the microscopic mechanism of how nitrogen bubbles affect the formation of THF hydrates at the molecular level. The results obtained reveal that the nitrogen bubble can promote the formation of THF hydrates. In the system with a nitrogen bubble, more THF-filled cages were generated, and the crystal structure was more orderly. The promotion of nitrogen bubbles on hydrate crystallization comes from the dissolution of nitrogen molecules. Some of dissolved nitrogen molecules can be enclosed in small hydrate cages near the nitrogen bubble, which can serve as stable sites for hydrate crystal growth, resulting in the fact that THF-filled cages connected with N₂-filled cages are much more stable and have a long lifetime. The results in this work can help to understand the promotion effect of micro- and nano-air bubbles on the crystallization of THF hydrates.     

Molecular vibrations of methane molecules in the structure I clathrate hydrate from ab initio molecular dynamics simulation

Vibrational frequencies of guest molecules in clathrate hydrates reflect the molecular environment and dynamical behavior of molecules. A detailed understanding of the mechanism for the vibrational frequency changes of the guest molecules in the clathrate hydrate cages is still incomplete. In this study, molecular vibrations of methane molecules in a structure I clathrate hydrate are calculated from ab initio molecular dynamics simulation. The vibrational spectra of methane are computed by Fourier transform of autocorrelation functions, which reveal distinct separation of each vibrational mode. Calculated symmetric and asymmetric stretching vibrational frequencies of methane molecules are lower in the large cages than in the small cages (8 and 16 cm(-1) for symmetric and asymmetric stretching, respectively). These changes are closely linked with the C-H bond length. The vibrational frequencies for the bending and rocking vibrational modes nearly overlap in each of the cages.     

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.     

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.     

Molecular dynamics study of gas hydrate formation

Long-time-scale molecular dynamics simulations are presented of the spontaneous formation of methane hydrate at a methane/liquid water interface. The water film was prepared at 300 K, 30 bar and showed no significant hydrate order. On crash cooling to 250 K, 300 bar (about 20 K subcooling), the system showed a rapid growth of hydrate clusters. Contrary to popular models for hydrate nucleation, the clusters formed first as two-dimensional arrangements and only later into three-dimensional cage structures; the results are, however, consistent with the local order model proposed recently. The hydrate clusters showed clear signatures of the type II hydrate structure even though the type I structure is the thermodynamically stable form for methane hydrate; this is in accord with the results of recent diffraction experiments.