金属Al表面熔化各向异性的分子动力学模拟

采用嵌入原子势,使用分子动力学方法对金属Al不同低指数晶面的表面熔化现象分别进行了模拟.分析了熔化过程中样品结构组态的变化.模拟结果表明对于不同的自由表面,表面熔化呈现出明显的各向异性行为.Al（110）面在低于熔点的温度之下发生预熔化;（111）与（001）面都出现过热现象.与（111）面不同,（001）面发生过热现象时表面原子层为类液层,而（111）面仍然保持很好的晶格结构,即预熔化的Al（001）面在高于熔点的温度下,仍可以在很长的时间内处于相对稳定的亚稳态.由模拟得到Al的热力学熔点为950 K左右,与实验值基本符合. 关键词： 分子动力学 表面熔化 过热     

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

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

Molecular dynamics of MgSiO3 perovskite melting

The melting curve of MgSiO分子动力学 MgSiO3钙钛矿 熔化温度 高压melting temperature, molecular dynamics, high pressureProject supported by the National Natural Science Foundation of China (Grant Nos 10274055 and 10376021),the Natural Science Foundation of Gansu Province, China (Grant No 3ZS051-A25-027) and the Scientific Research Foundation of Education Bureau of Gansu Province, China (Grant No 0410-01).2005-01-125/8/2005 12:00:00 AMThe melting curve of MgSiO3 perovskite is simulated using molecular dynamics simulations method at high pressure. It is shown that the simulated equation of state of MgSiO3 perovskite is very successful in reproducing accurately the experimental data. The pressure dependence of the simulated melting temperature of MgSiO3 perovskite reproduces the stability of the orthorhombic perovskite phase up to high pressure of 130GPa at ambient temperature, consistent with the theoretical data of the other calculations. It is shown that its transformation to the cubic phase and melting at high pressure and high temperature are in agreement with recent experiments.     

支撑Co团簇熔化行为的分子动力学模拟

运用分子动力学结合退火方法,模拟研究了所含原子数介于400～20000之间的自由及支撑Co团簇的熔化行为.模拟中团簇Co原子间相互作用采用半经验的Gupta多体相互作用势.讨论了不同衬底势(Lennard-Jones势与Morse势)对支撑团簇熔化行为的影响.结果表明:对于给定原子数目的支撑Co团簇,在衬底势强度相同时,两种衬底势下对支撑团簇熔点及预熔化区间给出了一致的描述;随衬底势强度增加,支撑团簇熔点升高,且两种衬底势所描述的支撑Co团簇也类同于自由Co团簇都表现出较好的熔点与尺寸依赖的线性关系.     

储氢笼型水合物生成促进机理的分子动力学模拟研究

用分子动力学(MD)模拟方法研究水合物法储氢的促进机理,系统研究纯H₂水合物、H₂+四氢呋喃(THF)水合物、H₂+四丁基溴化铵(TBAB)半笼型水合物和H₂+四异戊基溴化铵(TiAAB)半笼型水合物的微观结构及性质.模拟分析客体与笼子之间的稳定能ΔE_GH,得出水合物中大笼子对稳定水合物起到主要作用.THF进入大笼子能促进H₂水合物稳定,降低H₂水合物形成压力,模拟结果与实验一致.模拟对比不同客体在大笼子中的ΔE_GH值,得出从小到大的顺序依次为TiAAB,TBAB,THF,H₂.模拟结果表明半笼型水合物的稳定性比结构Ⅱ型水合物强,同时得出H₂+TiAAB半笼型水合物的结构最稳定.MD模拟为TiAAB成为一种水合物新型促进剂和新型储氢材料提供了理论依据. 关键词： 2笼型水合物')" href="#">H₂笼型水合物 分子动力学模拟 储氢 半笼型水合物     

纳米团簇熔化过程的分子动力学模拟

本文采用分子动力学结合嵌入原子多体势,模拟了不同半径的Ni纳米团簇的升温熔化过程,研究团簇尺寸对熔点和表面能的影响.模拟结果表明:团簇的熔点显著低于体材料的熔点.团簇熔化的过程首先是在团簇的表面出现预熔,然后向团簇内部扩展,直到整个团簇完全熔为液态.在模拟的纳米尺度范围内,团簇的熔点与团簇尺寸基本成线性关系.团簇的表面能随着团簇尺寸的增大而减小,而且表面能均高于体材料的表面能.     

Molecular dynamics simulations on the melting of gold nanoparticles

Molecular dynamics is employed to study the melting of bulk gold and gold nanoparticles. PCFF, Sutton-Chen and COMPASS force fields are adopted to study the melting point of bulk gold and we find out that the Sutton-Chen force field is the most accurate model in predicting the melting point of bulk gold. Consequently, the Sutton-Chen force field is applied to study the melting points of spherical gold nanoparticles with different diameters. Variations of diffusion coefficient, potential energy and translational order parameter with temperature are analyzed. The simulated melting points of gold nanoparticles are between 615～1115 K, which are much lower than that of bulk gold (1336 K). As the diameter of gold nanoparticle drops, the melting point also descends. The melting mechanism is also analyzed for gold nanoparticles.     

载能氢同位素原子与石墨(001)面碰撞的分子动力学研究

采用分子动力学方法研究了载能H同位素原子与石墨晶体碰撞的同位素效应. 碳氢系统的强共价键作用和石墨层间的弱van der Waals力分别用REBO和Ito半经验势函数来描述. 研究发现: 随着入射原子质量的增加, 上表面吸附几率和反射几率的峰值都会向高能区移动; 相比于H, ²H入射原子, ³H入射原子具有较高的吸附几率——包括上表面吸附和内部吸附; 穿透石墨晶体, ²H, ³H原子所需的能量较高; 原子质量和原子入射能量都会影响入射粒子与不同石墨层之间的能量传递过程. 这些结果对理解碳基材料的³H滞留机制有重要意义.     

分子动力学模拟Cu(010)基体对负载Co-Cu双金属团簇熔化过程的影响

纳米团簇负载到基体上的结构演化和热稳定性是其走向技术应用的关键. 本文用分子动力学结合嵌入原子方法模拟了具有二十面体初始结构的Co₂₈₁Cu₂₈₀ 混合双金属团簇在Cu(010)基体上的熔化过程, 考察了基体的Cu原子可以自由移动(自由基体)和固定(固定基体)两种条件对负载团簇熔化的影响. 发现基体条件对团簇的熔化有明显的影响. 在自由基体上团簇原子的温度-能量曲线存在明显的团簇熔化时的能量突变点, 熔点为1320 K, 低于固定基体上团簇的熔点1630 K. 在升温过程中团簇的二十面体结构会在基体表面发生外延生长. 外延团簇随着温度增加发生表面预熔, 预熔原子会逐渐向基体表面扩散形成薄层, 直至完全熔化. 自由基体上团簇原子的嵌入行为会使原子的分布状态产生不同于固定基体上的演变.     

铜-铝扩散焊及拉伸的分子动力学模拟

采用分子动力学方法模拟了铜-铝扩散焊过程，分析了理想平面铜-铝试件(001)晶面间扩散焊的过渡层厚度，并利用径向分布、键对分析方法分析了在不同的降温速率下过渡层的结构变化.降温速率大时，过渡层保持原有无序结构，降温速率小时，过渡层从无序结构向面心立方结构转变.还对扩散焊后的铜-铝试件进行了拉伸模拟，并与尺寸大小相近的单晶铜和单晶铝的拉伸模拟结果进行比较.结果发现焊接后的强度比单晶铝和单晶铜的强度都要小，最大应变值也小.     

Molecular dynamics study for the melting curve of MgO at high pressure

Shell-model molecular dynamics method is used to study the melting temperatures of MgO at elevated temperatures and high pressures using interaction potentials. Equations of state for MgO simulated by molecular dynamics are in good agreement with available experimental data. The pressure dependence of the melting curve of MgO has been calculated. The surface melting and superheating are considered in the correction of experimental data and the calculated values, respectively. The results of corrections are compared with those of previous work. The corrected melting temperature of MgO is consistent with corrected experimental measurements. The melting temperature of MgO up to 140GPa is calculated.     

应变效应对金属Cu表面熔化影响的分子动力学模拟

采用Mishin镶嵌原子势，通过分子动力学方法模拟了金属Cu 的(110)表面在不同应变下的熔 化行为，分析了表面熔化过程中系统结构组态和能量的变化以及固液界面迁移情况.金属Cu 的(110)表面在低于热力学熔点的温度下发生预熔化，准液体层的厚度随温度升高而增加.当 温度高于热力学熔点时，固液界面的移动速度与温度成正比，外推得到热力学熔点为1380K ，与实验结果1358K吻合良好.应变效应（包括拉伸和压缩）导致热力学熔点降低，并促进表 面预熔化进程.在相同温度条件下，准液体层的厚度随应变绝对值的增加而增大.应变效应导 致的固相自由能增加是金属Cu(110)表面热稳定性下降的主要因素，且表面应力和应变方向 的异同也会影响表面预熔化的进程. 关键词： 表面预熔化 热力学熔点 表面应力 分子动力学     

载能氢原子与石墨(001)面碰撞过程中的能量传递行为的分子动力学研究

本文采用分子动力学方法研究了单一载能氢原子与石墨碰撞时氢原子被石墨反射、 吸附和石墨被氢原子穿透的发生系数以及碰撞中的能量传递机理. 研究发现: 与单层石墨相比, 多层石墨之间的长程相互作用增加了氢原子发生反射的能量范围, 尤其当入射能量大于20.0 eV时, 对反射过程的影响很明显; 当氢原子的入射能量大于25.0 eV时, 有一定的概率穿透四层石墨; 当氢原子入射能量高于28.0 eV时, 载能氢原子的能量传递给第二层石墨烯的比传递给第一层石墨烯的多. 这些结果对理解聚变反应中, 碳基材料的化学腐蚀及氚滞留有重要意义.     

不同势下铱团簇结构和熔化行为的分子动力学模拟

采用分子动力学方法及淬火技术,结合Gupta势和Sutton-Chen势,模拟研究了包含13,14,55,56,147及148个原子数的铱团簇的熔化行为. 结果表明,Gupta势和Sutton-Chen势对所研究的Ir团簇的基态几何结构和熔化行为给出了基本一致的描述：两种不同势给出了完全相同的基态几何结构;两种势给出的Ir团簇熔点及预熔化区间随团簇尺寸的变化关系基本一致;对于小Ir_n团簇（n=13,14）两类势均表现出比热峰值相对于均方根键长涨落饱和值滞后的现象. 但是 关键词： 铱团簇 分子动力学 Gupta势和Sutton-Chen势 熔化     

Mechanism of microweld formation and breakage during Cu-Cu wire bonding investigated by molecular dynamics simulation

Currently,wire bonding is the most popular first-level interconnection technology used between the die and package terminals,but even with its long-term and excessive usage,the mechanism of wire bonding has not been completely evaluated.Therefore,fundamental research is still needed.In this study,the mechanism of microweld formation and breakage during Cu-Cu wire bonding was investigated by using molecular dynamics simulation.The contact model for the nanoindentation process between the wire and substrate was developed to simulate the contact process of the Cu wire and Cu substrate.Elastic contact and plastic instability were investigated through the loading and unloading processes.Moreover,the evolution of the indentation morphology and distributions of the atomic stress were also investigated.It was shown that the loading and unloading curves do not coincide,and the unloading curve exhibited hysteresis.For the substrate,in the loading process,the main force changed from attractive to repulsive.The maximum von Mises stress increased and shifted from the center toward the edge of the contact area.During the unloading process,the main force changed from repulsive to attractive.The Mises stress reduced first and then increased.Stress concentration occurs around dislocations in the middle area of the Cu wire.     

Molecular dynamics simulation model of hydrogen recycling on carbon divertor for neutral transport analysis in large helical device

To perform the neutral-transport simulation with processes in which hydrogen molecules contribute to the reaction such as molecular assisted recombination, the parameters of emitted neutral particles at the wall such as the energy distributions and the form (atom or molecule) of emitted neutral particles are necessary as a boundary condition of the calculation. Therefore, in order to provide information of recycled hydrogen on the divertor to neutral-transport code, molecular dynamics simulation of a hydrogen atom injection into a carbon material is performed to obtain the distributions of emission angle and translational energy of emitted hydrogen atoms and molecules. The distributions of rotational and vibrational energies are also investigated in the case of molecular hydrogen emission. Moreover, the quantum rotational state J, and vibrational state v are estimated from the classical value obtained by the simulation.     

Molecular dynamics simulation of surface melting behaviours of the V（110） plane

The modified analytic embedded-atom method and molecular dynamics simulations are applied to the investigation of the surface premelting and melting behaviours of the V（110） plane by calculating the interlayer relaxation, the layer structure factor and atomic snapshots in this paper. The results obtained indicate that the premelting phenomenon occurs on the V（110） surface at about 1800K and then a liquid-like layer, which approximately keeps the same thickness up to 2020K, emerges on it. We discover that the temperature 2020K the V（110） surface starts to melt and is in a completely disordered state at the temperature of 2140K under the melting point for the bulk vanadium.     

银纳米杆高温熔化断裂弛豫性能的原子级模拟研究

采用分子动力学模拟方法,研究了不同长度银纳米杆在不同温度弛豫过程中的结构演变过程.结果表明:银纳米杆存在一与杆长相关的临界熔断温度,该临界熔断温度随杆长增加而显著降低.当温度大于熔点而小于临界熔断温度时,体系形成一个高度无序的球形团簇,而温度大于临界熔断温度时,体系则熔断成两个球形团簇.并给出了银纳米杆的产生该熔断现象的机理. 关键词： 纳米杆 分子动力学 弛豫 熔化     

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

通过分子动力学方法模拟了在常温常压下(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%. 模拟结果表明, 往压水堆中的水加入氢气能明显地抑制水中的溶解氧.     

从常温常压到超临界乙醇的分子动力学模拟

采用分子动力学方法系统地研究了从常温常压到超临界状态乙醇的热力学性质、结构性质和动力学性质.模拟发现随着温度的升高,体系焓值增大,乙醇分子间的氢键作用减弱,自扩散系数增大;随着压强的增大,乙醇分子间的氢键作用增强,自扩散系数减小;乙醇自扩散系数在液相区随温度变化不明显,在气相区随压强增大很快减小,超临界区乙醇的自扩散系数比液相区大十几倍.温度和压强对乙醇自扩散系数的影响可通过密度来体现.与常温常压相比,超临界条件下的乙醇体系因密度涨落存在分子聚集现象,且在低密度区域更显著;乙醇分子间的氢键作用明显减弱,结