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

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

不同浓度聚乙烯醇对甲烷水合物分解作用的分子动力学模拟

采用分子动力学方法,模拟了273.15K下聚乙烯醇(PVA)对甲烷水合物的分解作用.研究发现,PVA浓度为2.7wt%时,水分子始终在其平衡位置附近波动,扩散系数仅为1.04×10-11m2/s;PVA浓度为5.2wt%时,水合物笼型结构坍塌,水分子以液态水的形式存在,甲烷分子从孔穴中逸出,聚乙烯醇的羟基在分子内部形成氢键,形成团簇,产生空间位阻,阻止了水分子再生成水合物,水分子的扩散系数1.61×10-9m2/s;PVA浓度为7.6wt%时,甲烷水合物周围有部分笼型结构被破坏,部分甲烷分子从孔穴中逸出,水分子扩散系数为3.55×10-10m2/s.得出聚乙烯醇对甲烷水合物的分解作用大小为:5.2wt%7.6wt%2.7wt%,为PVA溶液促进甲烷水合物分解实验研究提供参考.     

过氧化氢水溶液作用下甲烷水合物分解特性的分子动力学模拟研究

采用分子动力学模拟方法研究过氧化氢水(HP)溶液作用下结构I型(SI)甲烷水合物晶体分解特性. 系统分析甲烷水合物在过氧化氢水溶液作用下由晶态向液态转变过程的机理, 对比相同摩尔浓度乙二醇(EG)溶液作用下甲烷水合物分解变化规律, 得出HP与水合物热力学抑制剂EG一样对甲烷水合物分解具有促进作用, 为HP溶液促进甲烷水合物分解实验研究提供参考.     

客体分子所处晶穴结构不同对甲烷水合物稳定性的影响

采用分子动力学方法模拟了SⅠ型甲烷水合物受热分解微观过程,并对水合物分解过程中不同晶穴结构内客体分子对甲烷水合物稳定性的作用进行了研究.通过最终构象、均方位移和势能等性质的变化规律对分别缺失大晶穴和小晶穴中客体分子的2种水合物体系随模拟温度升高稳定性的变化进行了分析.模拟结果显示,随温度的上升,水合物稳定性逐渐下降直至彻底分解;而水合物分解速度与2种晶穴各自部分晶穴占有率相关,不能简单的通过整体晶穴占有率表示.对比相同注热过程中2种水合物体系分解状况,发现位于大晶穴内的客体分子对水合物稳定性影响更大,缺失大晶穴内客体分子的水合物更容易随温度升高而分解.     

多孔介质中气体水合物降压分解的分形理论研究

通过将水合物的分解过程看作是无固态产物层生成的气固反应过程, 结合粒径缩小的收缩核反应模型和分形理论, 建立了多孔介质中水合物降压分解的分数维动力学模型, 提出了基于水合物分解实验数据计算多孔介质分形维数的方法. 分别利用前人的甲烷水合物和CO2水合物降压分解实验数据, 对上述分数维动力学模型进行了验证. 计算结果表明, 用提出的方法所计算得到的多孔介质分形维数与前人的测定结果基本符合; 对甲烷水合物和CO2水合物的降压分解过程, 提出的分数维动力学分解模型得出了和实验结果基本一致的预测, 绝对平均误差(AAD)小于10%.     

第45届国际化学奥林匹克理论试题1

正第1题笼中枪(Clathrate gun)在海洋的底部,蕴藏着大量甲烷,这些甲烷以笼形包合物形式存在,被称为甲烷水合物。这些宝藏被开采后不仅可以提供能源,也可以用作有机合成的原料。但是,科学家们也非常担心由于海洋温度的升高而引发这些水合物的自发分解。科学家认为,如果一定量的甲烷被释放到大气中,其产生的温室效应会引     

甲烷水合物分解及自保护效应的分子动力学模拟

采用分子动力学(MD)方法, 在温度T = 240, 260, 280和300 K的条件下模拟了Ⅰ型甲烷水合物晶体的分解过程. 研究发现,水合物分解后将在相界面上形成一层“准液膜”,准液膜中水分子的结构性质、空间取向和动力学性质均出现由“似晶”到“似液”的渐变过程. 在水合物分解过程中, 准液膜的存在对水分子和甲烷分子的扩散形成传质阻力. 由于甲烷分子必须穿过准液膜才能进入气相, 准液膜的传质阻力抑制了甲烷分子向气相的扩散过程, 致使水合物的分解速率随之降低, 从而产生自保护效应. 当温度低于水的冰点时, 准液膜中水分子的“似晶”程度较高, 准液膜的传质阻力较大, 自保护效应较明显. 当温度高于水的冰点时, 准液膜中水分子的“似液”程度较高, 准液膜的传质阻力显著下降, 水合物的自保护效应明显减弱.     

南海神狐海域及祁连山冻土区天然气水合物的拉曼光谱特征

采用显微激光拉曼光谱对我国在南海神狐海域及祁连山冻土区首次钻获的天然气水合物实物样品进行了详细的研究, 探讨了其笼型结构特征及其气体组成. 结果表明, 南海神狐海域天然气水合物样品是典型的I型结构(sI)水合物, 气体组分主要是甲烷, 占99%以上; 水合物大笼的甲烷占有率大于99%, 小笼为86%, 水合指数为5.99. 祁连山冻土区天然气水合物气体组分相对复杂, 主要成分除甲烷外(70%左右), 还有相当数量的乙烷、丙烷及丁烷等烃类气体, 从拉曼谱图上可初步判断其为II型结构(sII)水合物; 水合物的小、大笼的甲烷占有率的比值(θ_S/θ_L)为26.38, 远远大于南海神弧海域水合物的0.87, 这主要是由于祁连山水合物气体组分中的大分子(乙烷、丙烷及丁烷等)优先占据水合物的大笼, 大大减少了大笼中甲烷分子的数量.     

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

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

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

甲烷水合物导热系数是甲烷水合物勘探、开采、储运以及其他应用过程中一个十分重要的物理参数.我们采用平衡分子动力学（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水合物具有相似的导热规律.压力对导热系数有一定影响,尤其是在较高压力条件下,压力越高,导热系数越大.而在不同温度和不同压力作用过程中,密度的改变对导热系数的增大或减小几乎没有影响.     

Time-resolved in situ neutron diffraction studies of gas hydrate: transformation of structure II (sII) to structure I (sI).

We report the in situ observation from diffraction data of the conversion of a gas hydrate with the structure II (sII) lattice to one with the structure I (sI) lattice. Initially, the in situ formation, dissociation, and reactivity of argon gas clathrate hydrate was investigated by time-of-flight neutron powder diffraction at temperatures ranging from 230 to 263 K and pressures up to 5000 psi (34.5 MPa). These samples were prepared from deuterated ice crystals and transformed to hydrate by pressurizing the system with argon gas. Complete transformation from D(2)O ice to sII Ar hydrate was observed as the sample temperature was slowly increased through the D(2)O ice melting point. The transformation of sII argon hydrate to sI hydrate was achieved by removing excess Ar gas and exposing the hydrate to liquid CO(2) by pressurizing the Ar hydrate with CO(2). Results suggest the sI hydrate formed from CO(2) exchange in argon sII hydrate is a mixed Ar/CO(2) hydrate. The proposed exchange mechanism is consistent with clathrate hydrate being an equilibrium system in which guest molecules are exchanging between encapsulated molecules in the solid hydrate and free molecules in the surrounding gas or liquid phase.     

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.     

Argon–water system at low temperatures

The molecular dynamics method is used to simulate argon solutions in water and a thin water film–argon system at low temperatures. The correlation in motions of two closely spaced argon atoms is of another nature than the correlation of two neon atoms in a neon solid solution in ice II. The structure of hydrate shells of argon atoms contains five-membered rings composed of water molecules. The solubility of argon in a water film at low temperatures is noticeably higher than at room temperature. If a water film is first cooled to the glassy state and then argon atoms are added to it, then approximately as many argon atoms are absorbed on the film surface as they are present in a cooled film in equilibrium with the argon atmosphere. Argon atoms migrate from one pit to another on the rough surface of a solid water film.     

Influences of hydration force and elastic strain energy on the stability of solid film in a very thin solid-on-liquid structure

Since hydration forces become very strong at short range and are particularly important for determining the magnitude of the adhesion between two surfaces or interaction energy, the influences of the hydration force and elastic strain energy due to hydration-induced layering of liquid molecules close to a solid film surface on the stability of a solid film in a solid-on-liquid (SOL) nanostructure are studied in this paper. The liquid of this thin SOL structure is a kind of water solution. Since the surface forces play an important role in the structure, the total free energy change of SOL structures consists of the changes in the bulk elastic energy within the solid film, the surface energy at the solid-liquid interface and the solid-air interface, and highly nonlinear volumetric component associated with interfacial forces. The critical wavelength of one-dimensional undulation, the critical thickness of the solid film, and the critical thickness of the liquid layer are studied, and the stability regions of the solid film have been determined. Emphasis is placed on calculation of critical values, which are the basis of analyzing the stability of the very thin solid film.     

NMR/MRI study of clathrate hydrate mechanisms

Clathrate hydrates are of great importance in many aspects. However, hydrate formation and dissociation mechanisms, essential to all hydrate applications, are still not well understood due to the limitations of experimental techniques capable of providing dynamic and structural information on a molecular level. NMR has been shown to be a powerful tool to noninvasively measure molecular level dynamic information. In this work, we measured nuclear magnetic resonance (NMR) spin lattice relaxation times (T1's) of tetrahydrofuran (THF) in liquid deuterium oxide (D2O) during THF hydrate formation and dissociation. At the same time, we also used magnetic resonance imaging (MRI) to monitor hydrate formation and dissociation patterns. The results showed that solid hydrate significantly influences coexisting fluid structure. Molecular evidence of residual structure was identified. Hydrate formation and dissociation mechanisms were proposed based on the NMR/MRI observations.     

Stabilization of methane hydrate by pressurization with He or N2 gas

The behavior of methane hydrate was investigated after it was pressurized with helium or nitrogen gas in a test system by monitoring the gas compositions. The results obtained indicate that even when the partial pressure of methane gas in such a system is lower than the equilibrium pressure at a certain temperature, the dissociation rate of methane hydrate is greatly depressed by pressurization with helium or nitrogen gas. This phenomenon is only observed when the total pressure of methane and helium (or nitrogen) gas in the system is greater than the equilibrium pressure required to stabilize methane hydrate with just methane gas. The following model has been proposed to explain the observed phenomenon: (1) Gas bubbles develop at the hydrate surface during hydrate dissociation, and there is a pressure balance between the methane gas inside the gas bubbles and the external pressurizing gas (methane and helium or nitrogen), as transmitted through the water film; as a result the methane gas in the gas bubbles stabilizes the hydrate surface covered with bubbles when the total gas pressure is greater than the equilibrium pressure of the methane hydrate at that temperature; this situation persists until the gas in the bubbles becomes sufficiently dilute in methane or until the surface becomes bubble-free. (2) In case of direct contact of methane hydrate with water, the water surrounding the hydrate is supersaturated with methane released upon hydrate dissociation; consequently, methane hydrate is stabilized when the hydrostatic pressure is above the equilibrium pressure of methane hydrate at a certain temperature, again until the dissolved gas at the surface becomes sufficiently dilute in methane. In essence, the phenomenon is due to the presence of a nonequilibrium state where there is a chemical potential gradient from the solid hydrate particles to the bulk solution that exists as long as solid hydrate remains.     

十二烷基苯磺酸钠在SiO2表面聚集的分子动力学模拟

采用分子动力学方法研究了阴离子表面活性剂十二烷基苯磺酸钠(SDBS)在无定形SiO2固体表面的吸附. 设置不同的水层厚度, 观察固液界面和气液界面吸附的差异. 模拟发现表面活性剂分子能够在短时间内吸附到SiO2表面, 受碳链和固体表面之间相互作用的影响形成表面活性剂分子层, 并依据吸附量的大小形成不同的聚集结构; 在水层足够厚的情况下, 由于有较多的表面活性剂分子吸附在固体表面,从而形成带有疏水核心的半胶束结构; 计算得到的成对势表明极性头与钠离子或水分子之间的结合或解离与二者之间的能垒有关, 解离能垒远大于结合能垒, 引起更多Na+聚集在极性头周围而只有少数Na+存在于溶液中; 无论气液还是固液界面, 极性头均伸向水相, 与水分子形成不同类型的氢键. 模拟表明, 分子动力学方法可以作为实验的一种补充, 为实验提供必要的微观结构信息.     

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.     

In CH3NH3PbI3 Perovskite Film,the Surface Termination Layer Dominates the Moisture Degradation Pathway

Theoretical studies have shown that surface terminations, such as MAI or PbI layers, greatly affect the environmental stability of organic–inorganic perovskite. However, until now, there has been little effort to experimentally detect the existence of MAI or PbI terminations on MAPbI₃ grains, let alone disclose their effects on the humidity degradation pathway of perovskite solar cell. Here, we successfully modified and detected the surface terminations of MAI and PbI species on polycrystalline MAPbI₃ films. MAI-terminated perovskite film followed the moisture degradation process from MAPbI₃ to hydrate MAPbI₃⋅H₂O and then into PbI₂, with penetration of water molecules being the main driving force leading to the degradation of MAPbI₃ layer by layer. In contrast, for the PbI-terminated perovskite film in a humid atmosphere, a deprotonation degradation pathway was confirmed, in which the film preferentially degraded directly from MAPbI₃ into PbI₂, here the iodine defects played a key role in promoting the dissociation of water molecules into OH⁻ and further catalyzing the decomposition of perovskite.