用分子动力学模拟甲烷水合物热激法结合化学试剂法分解

用分子动力学模拟方法研究甲烷水合物的热激法，化学试剂法，以及热激法结合化学试剂法分解，系统研究温度为277K和340K时添加液态水(WTR)和30wt％乙二醇(EG)溶液对水合物分解的影响.模拟显示WTR与水合物表面水分子形成氢键，破坏水合物原有的氢键平衡，造成笼状结构坍塌，水合物分解.EG分子中的羟基与水合物表面水分子形成氢键，从而破坏原有的稳定结构，造成水合物笼状结构被破坏，达到促进水合物分解，释放甲烷气体的效果.比较温度为277K和340K时添加WTR和30wt％EG溶液对水合物分解效果得出EG(340K)> WTR(340K)>EG(277K)>WTR(277K)，热激法结合化学试剂法能更好促进水合物分解.     

醇类抑制剂对甲烷水合物形成的影响

在本工作中,甲烷水合物的生长动力学是通过甲醇、乙醇、乙二醇三种不同醇类抑制剂存在下的分子动力学模拟研究的.模拟结果发现,三种醇类都可作为甲烷水合物的抑制剂,醇类分子中的亲水性羟基极大地破坏了水合物笼的结构,并且羟基可以与局部的液态水分子形成氢键,从而增加了形成水合物笼型结构的难度,导致甲烷水合物的生长速率降低.对于甲醇分子,甲醇分子的亲水性羟基与水分子形成氢键从而破坏了水分子结构,而亲油性甲基对周围的水分子具有簇效应,两者都会降低水合物生长速率;对于乙二醇和乙醇分子,它们只含有羟基,特别是乙二醇分子含有两个羟基,其对H₂O分子有很强的吸附作用,导致水合物生长速率降低.在抑制效果方面,甲醇分子最优,乙二醇稍微优于乙醇.     

高压下Ⅰ型甲烷水合物导热机理的分子动力学模拟

本文采用平衡态分子动力学方法研究了多种Ⅰ型甲烷水合物结构在高压下的导热性能。结果显示各水合物结构中水分子的排布构型差别不大,但各水合物结构的热导率存在一定差异。其中空穴水合物结构具有较好的导热性能,而含水分子晶格缺陷水合物结构的热导率较低。高压可以促进水合物中的声子传热,高温也能促进甲烷分子平动相关声子的传热,但将削弱水分子运动相关声子的传热。此外水合物结构中的水分子晶格缺陷会导致传热过程中声子大量散射。     

显微激光拉曼光谱测定甲烷水合物的水合指数

甲烷水合物是由甲烷气体分子与水分子在低温高压下形成的一种笼型结构化合物,广泛存在于海底陆架区和陆地冻土区,被认为是一种潜在的能源资源。在水合物的晶格中,水分子在氢键的作用下形成大小不同的笼子,甲烷分子可分别进入大笼(51262)和小笼(512)中。在自行研制的实验装置上,分别合成了一系列不同体系下的甲烷水合物,包括十二烷基硫酸钠(SDS)水溶液-甲烷体系、冰粉-甲烷体系以及冰粉-不同粒度砂-甲烷体系。对这些甲烷水合物样品进行了激光拉曼光谱分析,测定了其水合指数,笼占有率等结构参数。结果表明,这些甲烷水合物都为Ⅰ型结构,其水合指数和笼占有率基本不受沉积物粒径大小的影响。在3种体系中生成的水合物,大笼中甲烷分子基本占满,占有率大于97%;小笼中甲烷分子占有率为80%～86%,测得的水合指数为6.05～6.15。     

CH4水合物在SDS和CTAB溶液中的生成特性研究

为探究不同促进剂在甲烷水合物生成过程的微观作用机理,选取动力学促进剂十二烷基硫酸钠(SDS)和热力学促进剂十六烷基三甲基溴化铵(CTAB)作为添加剂,采用分子动力学方法研究其对甲烷水合物生成速率的影响.通过分析势能变化、均方位移、径向分布函数、分子簇生长速率,发现质量分数为0.9%SDS、1.2%SDS、1.2%CTAB、1.6%CTAB的溶液均可促进水合物生成.质量分数为1.2%的SDS溶液水合物生长速率最快,且SDS促进效果优于CTAB.通过分析甲烷分子密度分布云图,发现呈阴性的SDS分子头部基团吸附了大量甲烷分子,水分子受挤压向中间聚集;CTAB含氮的头部基团朝向均相溶液,包含在不稳定的水合物笼中,形成半笼型水合物.相比之下,CTAB溶液中水合物含气率更高.     

SⅡ型R134a水合物稳定性的分子动力学模拟北大核心

SⅡ型R134a制冷剂笼状水合物是理想的蓄冷介质之一。基于分子动力学(MD)方法,建立晶体结构模型,在NPT系综下运算,通过观察最终构象,分析计算体系中粒子的均方位移(MSD)和径向分布函数(RDF),从微观上揭示了小分子和温度对SⅡ型水合物晶体结构稳定性的影响。结果表明:小分子对小晶穴(512)占有率从0.25变化到1.00,占有率为0.25和0.50时,水合物晶体结构分解严重,体系中分子混乱,晶体结构难以稳定存在,占有率在0.75以上时,水合物晶体结构清晰有序,能够稳定存在;当体系温度由260K增大到290K时,水合物结构的稳定性逐渐降低,在280K时晶体结构出现扭曲变形,此时有部分水合物开始分解。     

SⅡ型R134a水合物稳定性的分子动力学模拟

SⅡ型R134a制冷剂笼状水合物是理想的蓄冷介质之一。基于分子动力学(MD)方法,建立晶体结构模型,在NPT系综下运算,通过观察最终构象,分析计算体系中粒子的均方位移(MSD)和径向分布函数(RDF),从微观上揭示了小分子和温度对SⅡ型水合物晶体结构稳定性的影响。结果表明:小分子对小晶穴(512)占有率从0.25变化到1.00,占有率为0.25和0.50时,水合物晶体结构分解严重,体系中分子混乱,晶体结构难以稳定存在,占有率在0.75以上时,水合物晶体结构清晰有序,能够稳定存在;当体系温度由260K增大到290K时,水合物结构的稳定性逐渐降低,在280K时晶体结构出现扭曲变形,此时有部分水合物开始分解。     

甲烷水合物成核过程的实验研究进展(特邀)

甲烷水合物(可燃冰)作为一种储量巨大、分布广泛的清洁能源而备受关注。围绕甲烷水合物的开采、以固态水合物形式储运天然气和氢气等问题开展基础科学研究,具有重要的科学意义和应用价值。成核过程是甲烷水合物形成的关键第一步,由甲烷、水分子形成团簇并逐渐演化形成水合物的微观过程。然而,由于缺乏在高压环境下研究成核微观过程的有效实验方法,针对成核过程的实验研究进展缓慢。本文首先对气体水合物的结构及其性质进行了回顾;然后,以水合物演化的分子动力学模拟为基础,梳理现有关于水合物演化路径的初步认知;最后,以超快非线性光谱学方法为主,讨论了水合物实验研究的进展和展望。     

光解水的原子尺度机理和量子动力学

利用阳光直接将水分解为不含碳的氢气燃料和氧气是面向全球能源危机环保且低成本的解决方案.得益于电子结构理论和量子模拟方法的进步,人们已经能够直接研究在纳米颗粒上等离激元诱导光解水过程在原子尺度上的反应机理和超快动力学.本文简述近年来的相关工作进展.吸附在氧化物薄膜上的金纳米颗粒很有希望成为水分解的高效新型光催化剂.在光激发条件下,水分解反应速率和光强、热电子转移之间有强相关性.水分解速率不仅取决于光吸收强度,还受到等离激元量子振动模式的调控.这对于太阳能光解水器件中纳米颗粒的设计有借鉴意义.我们发现液态水在金团簇等离激元催化下100 fs内就能产生氢气.超快量子动力学模拟表明,该过程中场增强起主导作用,从金属到水反键态的超快电荷转移也扮演着重要角色.综合这些原子尺度上的量子动力学研究,我们提出受激水分子中氢原子高速碰撞(速度远远超出其热速度)合成氢分子的"链式反应"机理.     

Simulation Study on the Structure and Dynamics of Water in Sodium Tetrafluoroborate/Water

应用分子动力学模拟方法研究了室温条件下四氟硼酸钠(NaBF₄)/水混合体系中水分子的微观结构、IR光谱以及转动动力学. 考察了混合物体系中水分子的摩尔分数浓度分别为6.25%、25.0%、50.0%、75.0%、90.0% 和99.6%时体系的结构和动力学性质. 研究显示在不同水分子含量的混合物中水分子以自由分子存在,随着混合物中水分子摩尔分数的增加,水分子的转动和弯曲振动带红移,而O-H伸缩振动蓝移,混合物中水分子内和分子间的氢键和非谐性相互作用增强,分子平动和转动变得困难和缓慢,研究结果与实验观测一致.     

显微激光拉曼光谱原位观测甲烷水合物生成与分解的微观过程

采用显微激光拉曼光谱技术对高压透明毛细管中甲烷水合物的生成与分解的微观过程进行了原位观测,初步探讨了甲烷水合物笼型结构的变化规律.结果表明,甲烷水合物在生成过程中,甲烷分子的拉曼峰(2 917 cm-1)逐渐分裂为两个峰(2 905和2 915 cm-1),表明溶解态甲烷分子从单一的化学环境进入了两个有差异的化学环境中...     

Methane combustion in hydrate systems: Water-methane and water-methane-isopropanol

Kinetics of dissociation of synthetic and natural methane gas hydrates, and also double isopropanol-methane hydrate is investigated. Thermal fields of the sample surfaces are measured by means of thermal imaging in combustion of released methane with clathrate dissociation. The dissociation rates of natural hydrate and double hydrate with isopropanol are many times lower than those of synthetic methane hydrate. Methane combustion is accompanied by formation of a thin water film on the powder surface, which has a strong effect on the heat and mass transfer mechanisms. The experiments demonstrated partial self-preservation for methane hydrate and the absence of self-preservation for double isopropanol-methane hydrate. The experimentally observed dissociation rate of double isopropanol-methane hydrate is considerably lower than that of methane hydrate.     

Molecular dynamics simulations of PVP kinetic inhibitor in liquid water and hydrate/liquid water systems

The possible effects of PVP (poly(N-vinylpyrrolidone)) on the properties of liquid and water in clathrate hydrate has been investigated using NVT molecular dynamics simulations. A model for a monomer of the PVP polymer is immersed in three systems, liquid water, a unit cell of a hydrate in liquid water with a hydrate former and a third system where some of the liquid water molecules of this last system are replaced by a PVP monomer. Both molecular dynamics simulation and integral equation theory predict hydrogen bonding between the double bonded oxygen in the PVP ring and hydrogens in water. For the composite system, the PVP monomer has a preference for hydrogen bonding to hydrogens from the water molecules at the surface of the hydrate lattice. The simulations indicate that the PVP monomer tends to orient perpendicular to the hydrate surface. For the model systems in this study PVP may form hydrogen bonds with liquid water through the double bonded oxygen in the ring. When a hydrate crystal is immersed in the liquid water phase this hydrogen bonding is shifted towards the hydrate due to a more favourable Coulomb interaction involving hydrogens from more than one water molecule at the hydrate surface. The PVP monomer has a preference for perpendicular orientation with respect to the hydrate surface. A scheme is suggested for the characterization of kinetic hydrate inhibitors based on molecular dynamics simulations and on three basic properties. In addition to the energy between the active groups of the inhibitor and hydrate water another point of focus is the free energy changes in the interactions between the inhibitor and water as the charges are changed from zero to the original model charges. In particular the difference between this integral for the (hydrate water)–(PVP monomer) interaction and the (liquid water)–(PVP inhibitor) interaction should reflect the driving forces in freezing the inhibitor out from the liquid water phase and onto the hydrate surface. The third property in the characterization scheme is the diffusivities of groups connecting to the hydrate crystal, relative to the diffusivities of the hydrate crystal. Results are presented from simulations where a small cavity with a methane model as a guest is immersed in liquid water with free methane molecules at a temperature of 150K. Changes in structure, diffusivities and energy indicate a tendency towards a more solid–like structurde around the cavity.     

Methane hydrate crystal growth in a porous medium filled with methane-saturated liquid water

The nucleation, growth and ageing of methane hydrate crystals were observed visually in a porous medium filled with liquid water presaturated with methane. The pore space dimensions of the porous medium were 1.0?×?10²?µm. The pressure?temperature conditions at which hydrate formation was initiated corresponded to system subcoolings of 3.4?K, 6.7?K, 12.3?K and 14.1?K, respectively, where the system subcooling denotes the difference of the system temperature from the triple methane?hydrate?water equilibrium temperature under a given pressure. Faceted (skeletal) hydrate crystals grew and bridged the pore spaces without intervention of a liquid water layer when the subcoolings were equal or smaller than 6.7?K. The faceted crystals may form a physical bonding with the walls of the porous medium. At the higher subcoolings, the dispersive formation of dendritic crystals and subsequent morphological change into particulate crystals were observed. The bridging of the dendritic crystals is unlikely in the absence of a large amount of additional methane supply due to the dispersive spatial distribution of the dendritic crystals that have dimensions smaller than those of the pore spaces. As a result of the interpretation of the observed variation in the crystal morphology of the methane hydrate formed in liquid water, the dependence of the crystal morphology on the magnitude of the mass transfer of methane molecules in liquid water observed in the porous medium was consistent with that previously observed in a bulk methane–water system.     

Theoretical investigation of structures and compositions of double neon-methane clathrate hydrates,depending on gas phase composition and pressure

The region of existence of neon clathrate hydrates is an actual problem of hydrate chemistry. The current work presents theoretical study of the equilibrium formation conditions of pure neon clathrate hydrates and double clathrate hydrates of neon-methane mixture. The structures and properties of double clathrate hydrates were described within the scope of the previously developed molecular clathrate hydrate model that takes into account the influence of guest molecules on the host lattice, interaction of guest molecules between themselves, and the possibility of multiple filling of host lattice cages by guest molecules. The model makes it possible to find an equilibrium state and thermodynamic properties of clathrate hydrates at given values of p and T. In the present work, we considered the properties of double clathrate hydrates in the range of pressures from 0 to 4 kbar at 250 K. The results of modeling have shown that the mass fraction of neon in double clathrate hydrate of Ne and CH₄ mixture of cubic structure I (sI) can reach 26%, and 22.5% in double hydrate of cubic structure II (sII) even at a low methane concentration (1%) in gas phase, at high pressure. It is shown that in double clathrate hydrates of the Ne and CH₄ mixture at high pressures, phase transition sII-sI can occur.     

Raman Spectroscopy of Nitrogen Clathrate Hydrates

在253 K和16 MPa的压力下,于实验室内合成了氮气水合物,用显微共焦拉曼光谱对其N-N和O-H键伸缩振动的光谱特征进行了研究．结果表明,氮气水合物中的N-N和O-H键的拉曼峰分别为2322.4和3092.1 cm^-1,与天然的空气水合物中的数据十分接近．另外,还测定了液氮和溶解于水中的氮分子中N-N键的拉曼峰值,分别为2326.6和2325.0 cm^-1．氮气笼型水合物分解的拉曼谱图表明,氮分子同时进入水合物的大笼和小笼中,但由于氮分子在大、小笼中的环境氛围十分接近,其拉曼位移相差不大,故拉曼谱图只能显示N-N键伸缩振动一个峰．     

Density functional theory study of structural stability for gas hydrate

Using the first-principles method based on the density functional theory(DFT),the structures and electronic properties of different gas hydrates(CO_2,CO,CH_4,and H_2) are investigated within the generalized gradient approximation.The structural stability of methane hydrate is studied in this paper.The results show that the carbon dioxide hydrate is more stable than the other three gas hydrates and its binding energy is-2.36 e V,and that the hydrogen hydrate is less stable and the binding energy is-0.36 e V.Water cages experience repulsion from inner gas molecules,which makes the hydrate structure more stable.Comparing the electronic properties of two kinds of water cages,the energy region of the hydrate with methane is low and the peak is close to the left,indicating that the existence of methane increases the stability of the hydrate structure.Comparing the methane molecule in water cages and a single methane molecule,the energy of electron distribution area of the former is low,showing that the filling of methane enhances the stability of hydrate structure.     

Application of low field NMR T2 measurements to clathrate hydrates

Low field (2 MHz) Nuclear Magnetic Resonance (NMR) proton spin–spin relaxation time (T₂) distribution measurements were employed to investigate tetrahydrofuran (THF)—deuterium oxide (D₂O) clathrate hydrate formation and dissociation processes. In particular, T₂ distributions were obtained at the point of hydrate phase transition as a function of the co-existing solid/liquid ratios. Because T₂ of the target molecules reflect the structural arrangements of the molecules surrounding them, T₂ changes of THF in D₂O during hydrate formation and dissociation should yield insights into the hydrate mechanisms on a molecular level. This work demonstrated that such T₂ measurements could easily distinguish THF in the solid hydrate phase from THF in the coexisting liquid phase. To our knowledge, this is the first time that T₂ of guest molecules in hydrate cages has been measured using this low frequency NMR T₂ distribution technique. At this low frequency, results also proved that the technique can accurately capture the percentages of THF molecules residing in the solid and liquid phases and quantify the hydrate conversion progress. Therefore, an extension of this technique can be applied to measure hydrate kinetics. It was found that T₂ of THF in the liquid phase changed as hydrate formation/dissociation progressed, implying that the presence of solid hydrate influenced the coexisting fluid structure. The rotational activation measured from the proton response of THF in the hydrate phase was 31 KJ/mole, which agreed with values reported in the literature.