First-principles molecular dynamics simulation of biased electrode/solution interface

First-principles molecular dynamics simulations have been carried out for water in contact with Pt(1 1 1) surface. To apply negative bias potential to the water/Pt interface, excess electrons were added to our slab model using the recently developed computational scheme called “effective screening medium (ESM)”. Water molecules located away from the surface reoriented themselves to screen the electric field, but they responded differently near the surface. Water molecules nearest to the surface, forming a distinct layered structure with the hydrogen atom directed to the surface, increased the density with increasing field. On these bases, we discuss microscopic aspects of the electric double layer.     

A modified equation of state for Xe at high pressures by molecular dynamics simulation

The exact equation of state （EOS） for the fission gas Xe is necessary for the accurate prediction of the fission gas behavior in uranium dioxide nuclear fuel, However, the comparison with the experimental data indicates that the applicable pressure ranges of existing EOS for Xe published in the literature cannot cover the overpressure of the rim fission gas bubble at the typical UO2 fuel pellet rim structure. Based on the interatomic potential of Xe, the pressure-volume-temperature data are calculated by the molecular dynamics （MD） simulation. The results indicate that the data of MD simulation with Ross and McMahan＇s potential [M. Ross and A. K. McMahan 1980 Phys. Rev. B 21 1658] are in good agreement with the experimental data. A preferable EOS for Xe is proposed based on the MD simulation. The comparison with the MD simulation data shows that the proposed EOS can be applied at pressures up to 550 MPa and 3 GPa and temperatures 900 K and 1373 K respectively. The applicable pressure range of this EOS is wider than those of the other existing EOS for Xe published in the literature.     

Molecular dynamics simulation of structural change at metal/semiconductor interface induced by nanoindenter

The structures of the Si/Cu heterogenous interface impacted by a nanoindenter with different incident angles and depths are investigated in detail using molecular dynamics simulation.The simulation results suggest that for certain incident angles,the nanoindenter with increasing depth can firstly increase the stress of each atom at the interface and it then introduces more serious structural deformation of the Si/Cu heterogenous interface.A nanoindenter with increasing incident angle(absolute value) can increase the length of the Si or Cu extended atom layer.It is worth mentioning that when the incident angle of the nanoindenter is between-45° and 45°,these Si or Cu atoms near the nanoindenter reach a stable state,which has a lower stress and a shorter length of the Si or Cu extended atom layer than those of the other incident angles.This may give a direction to the planarizing process of very large scale integration circuits manufacture.     

纳米铜团簇凝结规律的分子动力学研究

采用分子动力学方法对包含147，309和561个原子数的液态纳米铜团簇凝结过程进行模拟研 究，结果表明降温速率及团簇原子数对凝结得到常温下的固态团簇结构有重要影响：在模拟 时间内，降温速度越慢，团簇原子数越少，凝结得到铜团簇越倾向生成二十面体结构，反之 则倾向生成面心立方结构.同时探讨了该现象的物理机理. 关键词： 铜团簇 凝结 结构 分子动力学     

SiC/C界面辐照性能的分子动力学研究

本文通过分子动力学模拟的方法,研究了5种含不同空间结构的SiC/C界面的材料受辐照后的缺陷分布随时间以及PKA位置的变化关系,并与单质SiC中缺陷分布情况进行对比.利用径向分布函数分析了辐照对界面原子排列情况的影响.研究结果表明,SiC/C界面的抗辐照能力明显低于SiC内部,不同的空间结构对界面缺陷数量存在一定影响.由径向分布函数推得界面区域石墨原子密度高则界面原子排列情况受辐照影响越大.     

Surface diffusion of Si,Ge and C adatoms on Si （001） substrate studied the molecular dynamics simulation

Depositions of Si, Ge and C atoms onto a preliminary Si (001) substrate at different temperatures are investigated by using the molecular dynamics method. The mechanism of atomic self-assembling occurring locally on the flat terraces between steps is suggested. Diffusion and arrangement patterns of adatoms at different temperatures are observed. At 900 K, the deposited atoms are more likely to form dimers in the perpendicular [110] direction due to the more favourable movement along the perpendicular [110] direction. C adatoms are more likely to break or reconstruct the dimers on the substrate surface and have larger diffusion distances than Ge and Si adatoms. Exchange between C adatoms and substrate atoms are obvious and the epitaxial thickness is small. Total potential energies of adatoms and substrate atoms involved in the simulation cell are computed. When a newly arrived adatom reaches the stable position, the potential energy of the system will decrease and the curves turns into a ladder-like shape. It is found that C adatoms can lead to more reduction of the system energy and the potential energy of the system will increase as temperature increases.     

高能量短脉冲激光作用水分子的动力学模拟

基于高能量短脉冲激光作用水分子实验的复杂性,需用采取一种比较准确且较为方便的方法来预测实验结果是必要的.本文采用分子动力学方法(Molecular Dynamics)对高能量密度短脉冲激光加载下的超临界水分子进行热力学分析及结构研究.结果表明,随着能量快速加入加载区域,水分子系统的温度迅速提高.同时,由于分子迅速向四周扩散,在加载区域产生空泡.其中,只有25.5%的能量用以提高水分子系统的动能(温度),其余的能量都用于增大水分子系统的势能(水分子间距).伴随着温度的提高,水分子热运动加快,有序程度逐渐减弱,O—H间距增加,水分子间的氢键作用减弱,分子极性降低.分子动力学模拟方法研究短脉冲高能量激光对水分子的作用,对水下激光微加工研究有一定指导意义.     

利用分子动力学方法铁素体-渗碳体相界面效应探究

珠光体是十分重要的组织结构,因此本文构建了含铁素体-渗碳体相界面的模型,并采用分子动力学模拟方法模拟纳米压入的过程。通过对模拟结果的力学性能和组织结构分析,探究了铁素体-渗碳体相界面效应。研究发现,距铁素体-渗碳体晶界不同距离（位置压入）,在压入最初阶段,压头载荷随着压头与晶界距离的增大而增大,当压入深度达到一定深度后,载荷随着距离的增大而减小。杨氏模量和最大剪切模量受压头尖端下方原子结构的直接影响,硬度受到结构完整性和类型的共同影响。铁素体-渗碳体相界面影响了纳米压入过程中位错形核、增殖和扩展,宏观表现为在相同压入深度下,不同压入位置压头载荷的差异。     

镍基单晶高温合金γ/γ′（001）相界面上原子构型的分子动力学研究

用分子动力学方法研究了镍基单晶高温合金γ/γ′（001）相界面上三种各具特征的原子堆垛结构. 能量学计算发现，存在最优构型，动力学模拟显示不同构型的界面弛豫后，在相界面上都“成对”出现刃型错配位错. 相关计算表明体系能量、界面形成能及弛豫能都依赖于界面原子堆垛特征，而几何特征则具共性，即不同原子构型的界面具有同一的应力释放模式.     

单壁碳纳米管在石墨基底上运动的分子动力学模拟

采用分子动力学模拟方法研究单壁碳纳米管在石墨基底上的运动.首先碳纳米管在基底弛豫至平衡状态，然后对其施加一固定外力，撤去外力后，碳纳米管在基底上逐渐减速至停止.为了研究管径、手性角对运动方式的影响，本文选择了C(10,10)，C(10,9)，C(10,8)，C(10,5)，C(10,0)，C(8,8)六种单壁碳纳米管进行模拟.结果表明，碳纳米管在石墨基底上的运动方式由手性角决定，与管径无关.手性角等于30°时，碳纳米管与石墨基底之间为公度结构，碳纳米管的运动出现周期性的滑动和翻滚现象；手性角大于28.3°小于30°时，碳纳米管一边向前滑动一边滚动；手性角小于26.3°时，碳纳米管在基底上滑动.碳纳米管的手性角决定了它与石墨基底接触界面的微观构型，从而决定了碳纳米管的运动方式. 关键词： 分子动力学模拟 碳纳米管 动能 结构公度性     

镍基单晶高温合金γ/γ′（001）相界面上原子构型的分子动力学研究

用分子动力学方法研究了镍基单晶高温合金γ/γ′（001）相界面上三种各具特征的原子堆垛结构. 能量学计算发现，存在最优构型，动力学模拟显示不同构型的界面弛豫后，在相界面上都“成对”出现刃型错配位错. 相关计算表明体系能量、界面形成能及弛豫能都依赖于界面原子堆垛特征，而几何特征则具共性，即不同原子构型的界面具有同一的应力释放模式.     

广义等温等压系综-分子动力学模拟全原子水的气液共存形貌

探测相变过程中瞬时共存相的形貌等特征对理解其微观机制十分重要.本文应用广义等温等压系综-分子动力学模拟方法,研究全原子水模型的气液两相平衡及相变的中间过程.研究发现,此广义系综方法能够通过持续降温,连续地历经从气态、气液共存到液态的整个相变过程,通过持续升温历经其相反过程,而不会发生标准正则系综中的过饱和热滞现象.该方法不需要使用副本交换等增强抽样方法,因此可以用于较大体系的研究,多个独立的模拟即能获得整个气态液态区的平衡性质及共存相特征.本文还提出了计算气液共存界面面积的新方法,给出了水的气液共存界面形状随温度、压强的变化规律.结果表明,低压时水的气液共存界面因其较大的表面张力而接近球面,符合经典成核理论的描述,但随着压强升高接近其临界压强时,气态和液态的差别减小,界面的表面张力变小,界面形状变为无规则的枝杈结构,表现为二阶相变特征.     

纳米尺度下Si/Ge界面应力释放机制的分子动力学研究

采用分子动力学方法研究了纳米尺度下硅(Si)基锗(Ge)结构的Si/Ge界面应力分布特征,以及点缺陷层在应力释放过程中的作用机制.结果表明:在纳米尺度下, Si/Ge界面应力分布曲线与Ge尺寸密切相关,界面应力下降速度与Ge尺寸存在近似的线性递减关系;同时,在Si/Ge界面处增加一个富含空位缺陷的缓冲层,可显著改变Si/Ge界面应力分布,在此基础上对比分析了点缺陷在纯Ge结构内部引起应力变化与缺陷密度的关系,缺陷层的引入和缺陷密度的增加可加速界面应力的释放.参考对Si/Ge界面结构的研究结果,可在Si基纯Ge薄膜生长过程中引入缺陷层,并对其结构进行设计,降低界面应力水平,进而降低界面处产生位错缺陷的概率,提高Si基Ge薄膜质量,这一思想在研究报道的Si基Ge膜低温缓冲层生长方法中初步得到了证实.     

磁盘润滑膜全氟聚醚的分子动力学模拟研究

纳米润滑膜全氟聚醚(perfluoropolyether，简称PFPE)在固体表面的结构和迁移特性对于计 算机磁盘磁头驱动系统的可靠性具有重要的作用. 采用基于粗粒珠簧模型的分子动力学模拟 方法，考察了不同壁面和端基极性对于PFPE膜静态特性（如分子构型、单体密度分布、端基 密度分布）以及动态特性（如自扩散系数）的影响. 静态特性的模拟结果表明，非极性PFPE 膜在壁面附近具有单层厚度为一个单体直径的层状结构，而极性PFPE膜则具有复杂的层状结 构. 动态特性的模拟结果表明，PFPE膜的扩散能力因壁面作用而增强并因端基极性作用而减 弱. 关键词： 全氟聚醚膜 分子动力学模拟 薄膜润滑 固液界面结构     

Coarse-grained molecular dynamics simulations of capillary evaporation of water confined in hydrophilic mesopores

Non-equilibrium molecular dynamics (MD) simulations were performed to investigate the capillary evaporation of water confined in hydrophilic mesopores. The electrostatics-based (ELBA) coarse-grained water model was employed to calculate the duration of the time-consuming capillary evaporation process. To evaluate the effect of hydrophilicity of mesopores on the capillary evaporation of water, three types of thin films with a cylindrical mesopore were modelled by tuning the interactions between water and wall atoms. Initially, the cylindrical mesopore was filled with water, and evaporation of the water into vacuum was simulated. The calculation results showed that when capillary evaporation occurred, the desorption rate of water was almost constant in a highly hydrophilic mesopore where a stable water layer was formed on the pore surface, whereas the rate decreased with time in a weakly hydrophilic mesopore where the water layer did not remain stable. As time progressed, the water column shortened and then broke up. The number of water molecules in the mesopores decreased exponentially with time. The difference in the hydrophilicity of the mesopores resulted in different relaxation curves of water desorption from the mesopores.     

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

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

一种描述金属界面原子扩散的加速分子动力学方法

传统分子动力学(MD)的纳秒级时间尺度限制了对固体界面原子的深层扩散、渗透以及相形成等长时动力学性质的模拟研究.在Voter的超动力学框架内,提出了一种更为简单的偏移势的构建方法.该方法通过在偏移势中引入一个加速因子,抬高了原势阱,从而加速了原子的跃迁,将MD模拟的时间尺度提高了若干个数量级.更为重要的是,该方法不需要预知体系势能的势阱及鞍点分布,还能够将原势能曲面的特性完全保留.以Mg/Zn界面扩散为例,选取简单的Lennard-Jones双体势,考察了不同加速因子对界面原子扩散速度的影响.结果显示,该 关键词： 超动力学 加速因子 原子扩散 金属界面     

水的密度极值和相关性质的分子动力学研究

本文采用分子动力学方法对NVT系综中水的密度极值问题以及分布函数和扩散系数等相关性质进行了模拟.运用极值思想,通过引入新的能量-压力关联判据和线性插值的方法,可以将水的密度极值的温度值点上移到2℃,与实验值更加接近,相对于前人用TIP4P水分子模型的模拟结果-13℃或-23℃具有较好的改进,并且避开了临界点附近统计压力值的困难.模拟获得的径向分布函数和扩散系数的值与实验基本吻合,分子间作用能分布也与文献报道的趋势一致.     

Molecular dynamics simulation of gas diffusion behavior in polyethylene terephthalate/aluminium/polyethylene interface

Molecular dynamics (MD) simulations were performed to estimate the diffusion coefficients of O₂ and H₂O molecules in polyethylene terephthalate/aluminum/polyethylene interface at the temperature of 298 K. It came out that the diffusion coefficient of gasses in the interface is smaller than that of a single polymer, and the diffusion coefficients compare well with experimental data as well as previously published work. Furthermore, the diffusion coefficients of H₂O molecules in the interface are preferable to that of O₂ molecules. Interestingly, the largest diffusion coefficient was detected in the polyethylene terephthalate/aluminum(1 0 0)/polyethylene interface, while the smallest value of the diffusion coefficients was found in the polyethylene terephthalate/aluminum(1 1 1)/polyethylene interface. Calculation and analysis of the interaction between aluminum and polymers indicated that the interaction of polymer/aluminum(1 1 0) has the most interface strength, and crystal density of the metal surface has a definite effect on the planar interface energy. What’s more, the figure of gas molecule concentration is further resulted that the interface make contribution to adsorption of gas molecules. Moreover, the diffusion is belonging to the Einstein diffusion in the multilayer materials, and this work provides some key clues to improve the performance of polymer materials.     

Alternative ensemble averages in molecular dynamics simulation of hard spheres

ABSTRACT

We consider application to the hard sphere (HS) model of the mapped-averaging framework for generating alternative ensemble averages for thermodynamic properties. Specifically, we develop and examine new formulas for the pressure, the singlet and pair densities, and the cavity-correlation function inside the HS core; the pressure formula in particular is constructed such that it gives an ensemble average that exactly corrects the second-order virial equation of state. The force plays a central role in mapped-averaging expressions, and we write them in a way that accounts for the impulsive, event-driven nature of the HS dynamics. Comparison between results obtained conventionally versus mapped averaging finds that the latter has some advantage at low density, while both perform equally well (in terms of uncertainties for a given amount of sampling) at higher densities.