Molecular dynamics simulation of oxygen diffusion in dry and water-containing brown coal

Moisture content of brown coal has a pronounced impact on O₂ diffusion during coal oxidation. The transport behaviour of O₂ in dry and hydrated brown coal matrix with 10, 20, and 30?wt% moisture contents have been studied by molecular dynamics (MD) simulations at 298.15?K and 0.1?MPa. By analysing the diffusion characteristics of O₂, it is found that the diffusion process results from jumps of O₂ molecules between adjacent cavities in the coal matrix. With the increase of moisture content, the swelling extent of coal matrix increases in the range of 10–30?wt%, resulting in increased cavity volume, cavity interconnectivity, and mobility of coal macromolecule. These factors are all beneficial to the transport process of O₂ molecule. In addition, O₂ molecule prefers to reside in the available water-poor coal cavity and unoccupied region of water-rich cavity according to the estimated binding energies.     

Relaxation dynamics and power spectra of liquid water: a molecular dynamics simulation study

We present molecular dynamics simulations of liquid water at normal and supercooled conditions. Autocorrelation functions (ACFs) of several structural quantities and their fourier transforms are obtained and analysed. Structural correlations and relaxation times increase linearly with degree of supercooling. Power spectra of ACFs show increase in librational motion of liquid water with cooling. These modes intensify with supercooling because of structuring and ordering of water molecules. Overall, liquid water structure is homogenous over the temperatures and pressures studied and undergoes fluctuation–dissipation in its local-density variations [English and Tse, Phys. Rev. Lett. 106, 037801 (2011)].     

界面氢键对受限水结构和动态特性的影响

应用反应力场分子动力学方法, 模拟了水限制在全羟基化二氧化硅晶体表面间的弛豫过程, 研究了基底表面与水形成的界面氢键, 及其对受限水结构和动态特性行为的影响. 当基底表面硅醇固定时, 靠近基底表面水分子中的氧原子与基底表面的氢原子形成强氢键, 这使得靠近表面水分子中的氧原子比对应的氢原子更靠近基底表面, 从而水分子的偶极矩远离表面. 当基底表面硅醇可动时, 靠近基底表面水分子与基底表面原子形成两种强氢键, 一种是水分子中的氧原子与表面的氢原子形成的强氢键, 数量较少, 另一种是水分子中的氢原子与表面的氧原子形成的强氢键, 数量较多, 这使得靠近表面水分子中的氢原子比对应的氧原子更靠近表面, 从而水分子的偶极矩指向表面. 在相同几何间距下, 当基底表面硅醇可动时, 表面的活动性使得几何限制作用减弱, 导致了受限水分层现象没有固定表面限制下的明显. 此外, 固定表面比可动表面与水形成的界面氢键作用较强, 数量较多, 导致了可动表面限制下水的运动更为剧烈.     

改性活性炭内丙酮吸附的动力学模拟

本文通过在活性炭孔表面植入不同摩尔比例的羰基、羟基和羧基实现孔的改性,并采用分子动力学模拟研究了丙酮分子在改性活性炭孔内的动力学特性.模拟结果表明:改性后,孔内丙酮分子自扩散系数显著降低,其中羰基和羧基改性孔的丙酮分子扩散系数随着官能团的增多逐渐降低,但羟基改性孔的分子扩散系数并未呈现单调递减的趋势.植入官能团摩尔比例相同时,羧基对丙酮分子扩散的影响最高,羟基次之,羰基最低.总结得出改性孔内丙酮分子扩散主要受势能平滑度,自由体积分数以及分子与孔表面相互作用能三大因素的影响.     

水滴撞击结冰过程的分子动力学模拟

利用三维分子动力学模拟方法,研究了纳米尺度水滴撞击冷壁面的结冰过程.数值模拟中,统计系统采用微正则系综,势能函数选用TIP4P/ice模型,温度校正使用速度定标法,牛顿运动方程的求解采用文莱特算法,水滴内部结冰过程则通过统计垂直方向水分子温度分布来判定.研究发现,当冷壁面温度降低时,水滴完全结冰的时间减小,但水滴降至壁面温度的时间却增大;同时随着壁面亲水性降低,水滴内部热传递速度减慢(尤其是冷壁面与水滴底端分子层间),水滴内部温度趋于均匀,但水滴完全结冰时间延长.     

Pt(111)表面低能溅射现象的分子动力学模拟

利用嵌入原子方法的原子间相互作用势，通过分子动力学模拟，详细研究了贵金属原子在Pt (111)表面的低能溅射现象.模拟结果显示：对于垂直入射情况，入射原子的质量对Pt (11 1)表面的溅射阈值影响不大.当入射原子的能量小于溅射阈值时，入射原子基本以沉积为主 ；当入射原子的能量大于溅射阈值时，溅射产额随入射原子能量的增加而线性增大；当入射 原子能量达到200 eV时，各种入射原子的溅射产额都达到或接近1，此时入射原子主要起溅 射作用.溅射原子发射的角分布概率和溅射花样与高能溅射相类似.研究表明：与基于二体碰 撞近似的线性级联溅射理论不同，当入射原子能量大于溅射阈值时，低能入射原子的溅射产 额正比于入射原子的约化能量和入射原子与基体原子的质量比.通过对低能入射原子的钉扎 能力分析，提出了支配低能溅射的入射原子反射物理机理. 关键词： 分子动力学模拟、低能溅射     

Influence of water on HFO-1234yf oxidation pyrolysis via ReaxFF molecular dynamics simulation

ABSTRACT

The effect of water molecules on HFO-1234yf oxidation pyrolysis was investigated by ReaxFF-molecular dynamics simulation from 1900 to 4200?K. The initial pyrolysis of HFO-1234yf starts around 2500?K and the water molecules participate in chemical reactions at 2800?K when the reactants pyrolysis reached the highest reaction rate. The primary products including HF, COF₂ and CO₂ are observed at 2600, 2700 and 2900?K, respectively. The influence of water molecules on products is mainly reflected in the promotion activity on the conversion from COF₂ to CO₂ and the generation of HF molecules. Four formation pathways are observed and calculated to further elucidate the procedure of pyrolysis. The main conversion process from H₂O to HF is the ?F?+?H₂O?=?HF+?OH reaction, and the paths from H₂O to ?OH radical and COF₂ to ?CFO radical which are promoted by ?F and ?H radical, respectively, have relatively low energy barriers of 10.44 and 40.29?kJ/mol, and both reaction processes released HF molecules.     

石墨烯表面的特征水分子排布及其湿润透明特性的分子动力学模拟

石墨烯因其独特的分子构型、卓越的物理化学性能而受到广泛关注.本文首先利用分子动力学模拟比较了单层石墨烯、铜、二氧化硅三者表面的浸润性,除了接触角的比较,还分析了基底表面的水分子排布,得到石墨烯表面的特征水分子排布为:表面有两层密集的水分子层,其中靠近基底的密集水分子层中O—H键与垂直基底方向夹角集中在90°附近,并且基底表面的氢键几乎都垂直于基底.另一方面,本文研究了石墨烯浸润透明特性,发现在铜和二氧化硅上添加一层石墨烯,对铜的浸润性影响较小,对二氧化硅的浸润性影响很大,不仅使其上接触角显著增大,还使得基底表面的水分子排布呈现出类似单层石墨烯上的规律.本文使用分子动力学模拟方法从微观尺度验证了文献的实验结果,从基底表面水分子排布角度分析了石墨烯独特的浸润透明特性,为进一步开发石墨烯在微结构设计上的应用提供了理论指导.     

分子动力学模拟压水反应堆中联氨对水的影响

本文采用分子动力学方法模拟在常温常压下(1 atm,298 K)和在压水堆环境下(155 atm,626 K),水分子数为256,联氨(N₂H₄)分子数为0,25,50,75等不同数目时,水和联氨粒子系统的动力性质和微观结构.同时探讨了联氨分子的引入对水中溶解氧的影响.从模拟结果可知,在常温常压下,当联氨的分子数为0,25,50,75时,粒子系统的均方位移会随联氨分子数的增加而增加;联氨分子数为0与为25,50,75比较时会少一个数量级;压水堆环境下,联氨分子数 关键词： 分子动力学 压水堆 联氨     

Dynamical and structural properties of adsorbed water molecules at the TiO2 rutile-(110) surface: interfacial hydrogen bonding probed by ab-initio molecular dynamics

ABSTRACT

Ab-initio molecular dynamics (AIMD) simulations have been carried out to study a range of different and energetically-accessible adsorbed-water configurations and motifs for their vibrational and structural characteristics, in contact with rutile-(110) interfaces at 100?K. The radial pair distribution function between the titanium atoms at the interface and the hydrogen and oxygen atoms in the water monolayer show an orientation of the water molecules parallel to the surface of titania, and with hydrogen atoms pointed in the opposite direction to the surface. In some cases, a distinctive vibrational frequency region between 2500 and 3000?cm^?1 has also been observed, due to a strong dispersion interaction between water molecules. This behaviour is also seen in experimental studies of thin-film water coverage on TiO₂ surfaces.     

A study of cavitation nucleation in pure water using molecular dynamics simulation

To discover the microscopic mechanism responsible for cavitation nucleation in pure water, nucleation processes in pure water are simulated using the molecular dynamics method. Cavitation nucleation is generated by uniformly stretching the system under isothermal conditions, and the formation and development of cavitation nuclei are simulated and discussed at the molecular level. The processes of energy, pressure, and density are analyzed, and the tensile strength of the pure water and the critical volume of the bubble nuclei are investigated. The results show that critical states exist in the process of cavitation nucleation. In the critical state, the energy, density, and pressure of the system change abruptly, and a stable cavitation nucleus is produced if the energy barrier is broken and the critical volume is exceeded. System pressure and water density are the key factors in the generation of cavitation nuclei. When the critical state is surpassed, the liquid is completely ruptured, and the volume of the cavitation nucleus rapidly increases to larger than 100 nm³; at this point, the surface tension of the bubble dominates the cavitation nucleus, instead of intermolecular forces. The negative critical pressure for bubble nucleation is -198.6 MPa, the corresponding critical volume is 13.84 nm³, and the nucleation rate is 2.42×10³² m^-3·^-1 in pure water at 300 K. Temperature has a significant effect on nucleation: as the temperature rises, nucleation thresholds decrease, and cavitation nucleation occurs earlier.     

石墨烯纳米孔道内受限水介电常数影响因素的分子动力学模拟

受限水的介电特性对于石墨烯电容器的储能效率有着重大的影响.因此本文使用分子动力模方法,研究了不同石墨烯纳米孔道宽度(0.812～10 nm)下内受限水介电常数的分布及其影响因素.结果表明,在石墨烯纳米孔道内部受限水双电层结构可分为空隙层、界面层和体相层三个部分.其整体的介电常数随孔道大小的降低而线性减少.这是由内部体相层的宽度的降低而导致的.此外电极电压大小和石墨烯-水相互作用参数ε_CO的强弱也会显著改变电极表面的双电层(Electric Double Layers, EDLs)结构.其中电压的增大使得介电常数分布的震荡的程度也随之增加,最终导致了整体介电常数的减小.与之类似,亲水态的石墨烯表面(高ε_CO)下受限水分布的震荡程度也显著增加,这导致了整体介电常数的降低.     

分子动力学模拟研究方解石表面润湿性反转机理

利用分子动力学模拟技术从分子尺度探究方解石表面润湿性反转机理.首先,研究方解石表面润湿性反转过程;而后,从原油分子-方解石表面与原油分子-原油分子/水分子相互作用两个方面系统揭示方解石表面润湿性反转机理.结果:(1)水分子能够驱离方解石表面弱吸附的非极性分子造成润湿性的改变,但不能驱离强吸附的极性分子使润湿性反转难以实现;(2)原油分子极性越强与方解石表面相互作用越强,极性分子与方解石表面之间主要为静电力,非极性分子与方解石表面之间主要为范德华力;(3)原油分子极性越相近分子之间的相互作用越强,分子极性相差越大分子之间的相互作用越弱.非极性分子之间主要是范德华力,极性分子之间主要是静电力;(4)原油分子在方解石表面和水分子的共同作用下形成乙酸-吡啶-水-甲苯-己烷的稳定吸附序列.本研究为靶向提高采收率技术的设计与应用提供理论基础.     

低密度Ar熔化及结晶的分子动力学模拟

采用有位移力的ＬｅｎｎａｒｄＪｏｎｅｓ（１２－６）势对微正则系综下低密度Ａｒ系统（简约数密度为ρ＝０．８５）的一级相变过程进行了细致的分子动力学模拟,发现Ａｒ系统的熔化过程是原子的崩塌过程,结晶过程和理想的完整晶体不同,是一活化过程：形核长大过程随温度的降低进行,原来均匀分布在系统中的自由体积呈集中分布,由此系统达到了更稳定的结晶状态。     

接触角与液固界面热阻关系的分子动力学模拟

本文利用分子动力学方法模拟了液体在固体表面的 接触角及液固界面热阻, 并探讨了二者之间的关系. 通过分别改变液固结合强度和固体的原子性质来分析接触角和界面热阻的关系及变化趋势. 模拟结果显示增强液固间相互作用时, 接触角减小的同时界面热阻也随之单调减小; 而改变固体原子间结合强度和原子质量时, 接触角几乎保持不变, 但界面热阻显著改变. 固体原子间结合强度和原子质量影响界面热阻的原因是其改变了固体的振动频率分布, 导致液固原子间的振动耦合程度发生变化. 本文的结果表明界面热阻不仅与由接触角所表征的液固结合强度有关, 还与液固原子间的振动耦合程度有关. 接触角与界面热阻间不存在单值的对应关系, 不能单一地将接触角作为液固界面热阻的评价标准. 关键词： 液固界面 接触角 界面热阻 分子动力学模拟     

Investigation of the dissociative adsorption for cyclopropane on the copper surface by density functional theory and quantum chemical molecular dynamics method

The dissociative adsorption of cyclopropane on the copper surface was studied using quantum chemical molecular dynamics method with “Colors-Excite” code and density functional theory by Amsterdam Density Functional program (ADF2000). The excited state of cyclopropane was used as adsorbate to simulate the dissociated adsorption under an irradiation energy of ca. 10 eV. One of the C-C bonds in cyclopropane was broken and the two new bonds between cyclopropane and copper surface were formed. The electrons transferred from the copper atoms to cyclopropane with a value of about 0.2e. The shorter distances between the carbons and surface copper atoms showed the existence of strong interaction. Consistently, the results indicated metallacyclopentane was the most possible intermediate species in dissociative adsorption by ADF2000 and “Colors-Excite” method.     

High-temperature dynamic behavior in bulk liquid water: A molecular dynamics simulation study using the OPC and TIP4P-Ew potentials

Classical molecular dynamics simulations were performed to study the high-temperature (above 300 K) dynamic behavior of bulk water, specifically the behavior of the diffusion coefficient, hydrogen bond, and nearest-neighbor lifetimes. Two water potentials were compared: the recently proposed “globally optimal” point charge (OPC) model and the well-known TIP4P-Ew model. By considering the Arrhenius plots of the computed inverse diffusion coefficient and rotational relaxation constants, a crossover from Vogel–Fulcher–Tammann behavior to a linear trend with increasing temperature was detected at T* ≈ 309 and T* ≈ 285 K for the OPC and TIP4P-Ew models, respectively. Experimentally, the crossover point was previously observed at T* ± 315–5 K. We also verified that for the coefficient of thermal expansion α _P (T, P), the isobaric α _P (T) curves cross at about the same T* as in the experiment. The lifetimes of water hydrogen bonds and of the nearest neighbors were evaluated and were found to cross near T*, where the lifetimes are about 1 ps. For T < T*, hydrogen bonds persist longer than nearest neighbors, suggesting that the hydrogen bonding network dominates the water structure at T < T*, whereas for T > T*, water behaves more like a simple liquid. The fact that T* falls within the biologically relevant temperature range is a strong motivation for further analysis of the phenomenon and its possible consequences for biomolecular systems.     

分子动力学模拟压水反应堆中氢气对水的影响

通过分子动力学方法模拟了在常温常压下(1 atm, 298 K)和在压水堆环境下(155 atm, 626 K), 水分子数为256, 氢分子数为0, 25, 50, 75和100等不同数目时, 粒子系统的动力学性质和微观结构, 分析了不同氢气对水中溶解氧的影响. 从模拟结果可知, 在常温常压和压水堆环境下, 当氢粒子数分别为0, 25, 50, 75和100时, 粒子系统的均方位移会随氢分子数增加而增加, 并且常温常压下的增长幅度远小于压水堆环境下的增长幅度, 如压水堆环境下氢分子数为75时系统的均方位移约是常温常压下氢分子数为75时系统的均方位移的6.02倍, 比压水堆环境下氢分子数0时系统的均方位移增加了131.88%. 此外, 粒子系统的微观结构, 从径向分布函数看, 在常温常压下随着氢分子数目的增加而小幅度增加, 这与常温常压下因氢气溶解在水中增大了氧离子周围的粒子密度相符合. 而在压水堆环境下, 氢分子数为75, 50, 25与为0时的水比较, 其径向分布均不会有太大的变化, 而分子数为100时会出现明显增加, 与为0时的水比较其径向分布增加了22.00%. 模拟结果表明, 往压水堆中的水加入氢气能明显地抑制水中的溶解氧.     

Plasma-based introduction of monosort functional groups of different type and density onto polymer surfaces. Part 1: Behaviour of polymers exposed to oxygen plasma

Several approaches were investigated to produce monosort functionalized polymer surfaces with a high density and homogeneity of functional groups: (i) Plasma oxidation followed by wet-chemical reduction, (ii) formation of radicals and grafting on of functional group carrying molecules, (iii) plasma bromination followed by (iv) Williamson or Gabriel-like synthesis of spacer molecules, and (v) a pulsed plasma polymerization of functional groups bearing monomers or (vi) their copolymerisation with other comonomers. The formation of hydroxyl (OH), primary amino (NH₂), and carboxyl (COOH) groups was studied in detail. The oxygen plasma treatment (i) in a low-pressure non-isothermal glow discharge results in the formation of a wide variety of O functional groups, polymer degradation and crosslinking. Low power densities and short exposure times (0.1 to 2 s) are required to functionalize a surface while preserving the original polymer structure. Carbonate, ester, and aromatic groups are rapidly degraded by an oxygen plasma treatment leading to scissions of polymer backbones and loss in molecular weight. Also the formation of macrocycles and C=C bonds was observed in a region of around 4 nm in depth. The investigated polymers could be classified by their degradation behaviour on exposure to the oxygen plasma. In order to maximize the process selectivity for OH groups, the variety of oxygen functionalities formed by the oxygen plasma was wet-chemically reduced by diborane, vitride? (Na complex), and LiAlH₄. Typical yields were 9 to 14 OH groups per 100 carbon atoms. Plasma bromination (iii) (40 Br per 100 C atoms) of polymers, followed by grafting of spacer molecules (iv), has been proved to be a highly selective reaction. Another way to produce high densities of monosort functionalities was the pulsed plasma polymerization of functional group bearing monomers such as allylamine, allylalcohol or acrylic acid (v). The retention of chemical structure and functional groups during plasma polymerization was achieved by using low power densities and the pulsed plasma technique. The maximum yields were 30 OH, 18 NH₂, and 24 COOH groups per 100 C atoms. To vary the density of functional groups a chemical copolymerization with 'chain-extending' comonomers such as butadiene and ethylene was initiated in the pulsed plasma (vi). Additionally, the often-observed post-plasma oxidations of such layers initiated by reaction of trapped radicals with oxygen from the air were successfully suppressed by using NO gas as radical quencher.