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.     

First determination of volume changes and enthalpies of the high-pressure decomposition reaction of the structure H methane hydrate to the cubic structure I methane hydrate and fluid methane

Pressure-temperature (P-T) conditions of the decomposition reaction of the structure H high-pressure methane hydrate to the cubic structure I methane hydrate and fluid methane were studied with a piston-cylinder apparatus at room temperature. For the first time, volume changes accompanying this reaction were determined. With the use of the Clausius-Clapeyron equation the enthalpies of the decomposition reaction of the structure H high-pressure methane hydrate to the cubic structure I methane hydrate and fluid methane have been calculated.     

Molecular insights into the heterogeneous crystal growth of si methane hydrate

In this paper we report a successful molecular simulation study exploring the heterogeneous crystal growth of sI methane hydrate along its [001] crystallographic face. The molecular modeling of the crystal growth of methane hydrate has proven in the past to be very challenging, and a reasonable framework to overcome the difficulties related to the simulation of such systems is presented. Both the microscopic mechanisms of heterogeneous crystal growth as well as interfacial properties of methane hydrate are probed. In the presence of the appropriate crystal template, a strong tendency for water molecules to organize into cages around methane at the growing interface is observed; the interface also demonstrates a strong affinity for methane molecules. The maximum growth rate measured for a hydrate crystal is about 4 times higher than the value previously determined for ice I in a similar framework (Gulam Razul, M. S.; Hendry, J. G.; Kusalik, P. G. J. Chem. Phys. 2005, 123, 204722).     

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.     

Influence of pyrogenic silica hydration on water diffusion in surface layers

The surface hydration of pyrogenic silica (aerosil) has been studied by slow neutron scattering. It was shown that the rate of diffusion changes with the thickness of the layer monotonically. The mean square displacement of water molecules from the equilibrium position in aerosil hydration shell is smaller than in bulk water, but twice larger than in ice.     

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.     

Search for memory effects in methane hydrate: structure of water before hydrate formation and after hydrate decomposition

Neutron diffraction with HD isotope substitution has been used to study the formation and decomposition of the methane clathrate hydrate. Using this atomistic technique coupled with simultaneous gas consumption measurements, we have successfully tracked the formation of the sI methane hydrate from a water/gas mixture and then the subsequent decomposition of the hydrate from initiation to completion. These studies demonstrate that the application of neutron diffraction with simultaneous gas consumption measurements provides a powerful method for studying the clathrate hydrate crystal growth and decomposition. We have also used neutron diffraction to examine the water structure before the hydrate growth and after the hydrate decomposition. From the neutron-scattering curves and the empirical potential structure refinement analysis of the data, we find that there is no significant difference between the structure of water before the hydrate formation and the structure of water after the hydrate decomposition. Nor is there any significant change to the methane hydration shell. These results are discussed in the context of widely held views on the existence of memory effects after the hydrate decomposition.     

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.     

Comparison of Ethene Diffusion Characteristics in H[Al]ZSM‐5 at 300 K and 400 K by Molecular Dynamics

The ethene diffusion characteristics in the framework of H[Al]ZSM‐5 at 300 K and 400 K have been studied by molecular dynamics (MD) simulation. The data have been obtained for the molecular kinetic, potential and total energies, mean square displacement and self‐diffusion coefficient, interaction and beat of adsorption. The dependence of molecular diffusion on temperature has been explored.     

聚乙烯唑啉作用下甲烷水合物分解的分子动力学模拟

利用分子动力学模拟系统研究了不同质量浓度下(1.25%、2.50%、6.06%)聚乙烯唑啉(PEtO)对甲烷水合物的分解作用. 模拟体系为甲烷水合物2′2′2的超胞和聚合物对接体系. 模拟发现水分子间氢键构架的水合物笼型结构在PEtO的作用下出现扭曲, 最终导致水合物笼型结构完全坍塌. 通过氧原子径向分布函数、均方位移以及扩散系数比较不同浓度PEtO的作用, 证实在一定浓度范围内, PEtO的浓度越高, 其水合物分解作用越强. 此外, PEtO 具有一定的可生物降解性. PEtO 对水合物的作用为: PEtO 吸附在水合物表面, 其中的酰胺基(N―C＝O)与成笼的水分子形成氢键, 破坏邻近的笼形结构, 令水合物分解; PEtO不断分解表面的水合物, 直到水合物笼完全分解.     

Experimental and modeling study on decomposition kinetics of methane hydrates in different media

The decomposition kinetic behaviors of methane hydrates formed in 5 cm3 porous wet activated carbon were studied experimentally in a closed system in the temperature range of 275.8-264.4 K. The decomposition rates of methane hydrates formed from 5 cm3 of pure free water and an aqueous solution of 650 g x m(-3) sodium dodecyl sulfate (SDS) were also measured for comparison. The decomposition rates of methane hydrates in seven different cases were compared. The results showed that the methane hydrates dissociate more rapidly in porous activated carbon than in free systems. A mathematical model was developed for describing the decomposition kinetic behavior of methane hydrates below ice point based on an ice-shielding mechanism in which a porous ice layer was assumed to be formed during the decomposition of hydrate, and the diffusion of methane molecules through it was assumed to be one of the control steps. The parameters of the model were determined by correlating the decomposition rate data, and the activation energies were further determined with respect to three different media. The model was found to well describe the decomposition kinetic behavior of methane hydrate in different media.     

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.     

水合物CH4-CO2置换的双重机理:空位辅助与竞争

过去的CO2置换甲烷水合物的微观机理研究,主要集中在客体分子(CH4、CO2)之间的交换、占据状态,孤立地研究分解过程或生成过程,忽视主体-客体之间的作用、主体分子(H2O)的空位辅助和客体分子的多重竞争通道。本文基于水合物分解的过冷水及其水空位辅助,以及水合物生成的串滴链及其客体分子竞争的研究,进一步评论水合物CH4-CO2置换的双重机理。然后,对微观机理的动态性和未来研究的相关问题进行讨论。通过综述和评论,文章得出以下初步结果:过冷水通过水空位推动客体分子的跳跃、扩散,实现置换过程的自组装;CO2分子在分解前沿形成一个有序结构的CO2串滴链,其动态性伴随水的组织到获取包合物笼的结构,以及非晶形包合物转变成晶形包合物的生长过程;CO2和CH4在中晶穴中必然产生竞争,并且存在多种竞争类型;成核过程中,不稳定簇导致竞争结构,且有多重竞争通道。最后,结果表明水合物的CH4-CO2置换机理具有双重性,即主体分子的空位辅助和客体分子的竞争,是分解过程和生成过程的自然统一。     

NMR study of methane + ethane structure I hydrate decomposition

The thermally activated decomposition of methane + ethane structure I hydrate was studied with use of 13C magic-angle spinning (MAS) NMR as a function of composition and temperature. The observed higher decomposition rate of large sI cages initially filled with ethane gas can be described in terms of a model where a distribution of sI unit cells exists such that a particular unit cell contains zero, one, or two methane molecules in the unit cell; this distribution of unit cells is combined to form the observed equilibrium composition. In this model, unit cells with zero methane molecules are the least stable and decompose more rapidly than those populated with one or two methane molecules leading to the observed overall faster decomposition rate of the large cages containing ethane molecules.     

石墨烯狭缝受限孔道中水分子的分子动力学模拟

石墨烯是一种具有广泛应用前景的纳米材料,特别是由石墨烯片层自组装形成的二维纳米通道能够应用于物质的过滤分离.本文采用分子动力学模拟方法研究了原态石墨烯/羟基改性石墨烯狭缝孔道中水分子的微观行为,模拟计算了水的界面结构性质和扩散动力学性质,所研究的石墨烯孔宽为0.6-1.5 nm.模拟结果表明,在石墨烯狭缝孔道中,水分子受限结构呈现层状分布,在超微石墨烯孔道(0.6-0.8 nm)中水分子可形成特殊的环状有序结构,石墨烯表面可诱导产生特殊的水分子界面取向.在石墨烯孔道中,水分子的扩散运动低于主体相水分子的扩散运动,羟基化石墨烯孔道可以促使水分子的扩散能力降低.对于改性石墨烯狭缝孔道,由于羟基的作用,水分子可以自发渗入0.6 nm的石墨烯孔道内.模拟所得到的受限水分子的动力学性质与水分子在石墨烯孔道内的有序结构有关.本文研究结果将有助于分析理解水分子通过石墨烯纳米通道的渗透机理,为设计基于石墨烯的纳米膜提供理论指导.     

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.     

Thermodynamic properties of methane hydrate in quartz powder

Using the experimental method of precision adiabatic calorimetry, the thermodynamic (equilibrium) properties of methane hydrate in quartz sand with a grain size of 90-100 microm have been studied in the temperature range of 260-290 K and at pressures up to 10 MPa. The equilibrium curves for the water-methane hydrate-gas and ice-methane hydrate-gas transitions, hydration number, latent heat of hydrate decomposition along the equilibrium three-phase curves, and the specific heat capacity of the hydrate have been obtained. It has been experimentally shown that the equilibrium three-phase curves of the methane hydrate in porous media are shifted to the lower temperature and high pressure with respect to the equilibrium curves of the bulk hydrate. In these experiments, we have found that the specific heat capacity of the hydrate, within the accuracy of our measurements, coincides with the heat capacity of ice. The latent heat of the hydrate dissociation for the ice-hydrate-gas transition is equal to 143 +/- 10 J/g, whereas, for the transition from hydrate to water and gas, the latent heat is 415 +/- 15 J/g. The hydration number has been evaluated in the different hydrate conditions and has been found to be equal to n = 6.16 +/- 0.06. In addition, the influence of the water saturation of the porous media and its distribution over the porous space on the measured parameters has been experimentally studied.     

New effective method for quantitative analysis of diffusion jumps, applied in molecular dynamics simulations of small molecules dispersed in short chain systems

Diffusion jumps of small molecules dispersed in chain molecules or other kinds of slow-moving matrices have already been observed in many previous simulations of such systems, and their treatment led to important qualitative conclusions. In the present work, a new, very simple yet effective method is described, allowing for both identification of individual penetrant jump events and their quantitative treatment in a statistical sense. The method is applied in equilibrium Molecular Dynamics simulations for systems of gaseous alkanes, methane through n-butane, including also a mixture of methane and n-butane, dispersed in n-decane or n-eicosane. Equilibration and attainment of a linear diffusion regime is confirmed by means of various criteria, and the jumps detection method is applied to all systems studied. The results obtained clearly show the existence of distinct jump events in all cases, although the average jump length is reduced with penetrant or liquid alkane molecular weight. The method allows one to determine the average jump length and the corresponding jumps frequency. On the basis of these results, it was possible to estimate a random walk type diffusion coefficient, D(s,jumps), of the penetrants, which was found to be substantially lower compared with the overall diffusion coefficient D(s,MSD) obtained by the mean square displacement method. This finding led us to assume that the overall penetrants' diffusion in the studied systems is a combination of longer jumps with a smoother and more gradual displacement, a result that confirms assumptions suggested in previous studies.     

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.