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

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

低温下Cu13团簇负载于Cu(001)表面上结构变化的分子动力学研究

应用分子动力学方法研究温度为10和50 K时具有二十面体结构的Cu13团簇以不同接触条件与Cu(001)表面结合后的结构变化,原子间的相互作用势采用Johnson的嵌入原子方法模型.通过基于原子密度分布函数的分析表明,负载团簇与表面的结合能主要受团簇与载体相接触的最低层原子数及这些原子所具有的不同几何构型影响,同时更高层的原子呈现出不同的几何结构.温度为10 K时,负载团簇的初始位置对团簇几何结构和结合能影响较大.     

低温下Cu13团簇负载于Cu(001)表面上结构变化的分子动力学研究

应用分子动力学方法研究温度为10和50 K时具有二十面体结构的Cu₁₃团簇以不同接触条件与Cu（001）表面结合后的结构变化,原子间的相互作用势采用Johnson的嵌入原子方法模型.通过基于原子密度分布函数的分析表明,负载团簇与表面的结合能主要受团簇与载体相接触的最低层原子数及这些原子所具有的不同几何构型影响,同时更高层的原子呈现出不同的几何结构.温度为10 K时,负载团簇的初始位置对团簇几何结构和结合能影响较大. 关键词： 分子动力学 团簇 表面 计算机模拟     

Cu偏析诱导Co团簇结构及性质异常的动力学模拟

采用分子动力学结合镶嵌原子势方法,模拟研究了Cu原子分别分布于基体Co团簇内层和表面构成Cu-Co合金团簇的结构和热力学性质,研究表明,相同数目的 Cu原子掺杂到基体中因掺杂层的不同,会诱导内层Co团簇和外层Co团簇结构、能量及熔点表现出巨大差异;Cu原子在团簇各层掺杂位置的差异,会导致原子向低能态位置偏移,但相对移动后后续原子的补位,使团簇结构随温度呈相对无扩散度相变;Cu原子由内层向表面偏析是内层Co团簇与相同原子数比例的外层Co团簇熔点产生巨大差异的主要原因.     

利用Gupta势研究ConCu13-n (n=1～12)混合团簇的基态结构和熔化行为

利用Gupta多体相互作用势结合遗传算法和分子动力学方法模拟研究了ConCu13-n(n==0～13)单质及混合团簇的基态结构和熔化行为,结果表明:ConCu13-n(n=1～12)混合团簇的基态结构均是在单质Co13、Cu13基态二十面体基础之上的畸变,Co原子先占据中心后占据表面,表面上的Co原子总连接在一起,抱团分布;分析二级差分能和混合能发现Co1Cu12、Co7C6具有相对高的稳定性,可视为幻数结构团簇;ConCu13-n(n=1～12)混合团簇的熔点均位于Co13、Cu12单质团簇的熔点之间,且随着Co原子数目的增多总体呈上升趋势,但在n=3时出现反常(熔点降低),这可归因于Co3Cu10的基态与第一激发态之间的能量差远小于Co2Cu11的相应值.     

铂纳米粒子热稳定性的原子级模拟研究

本文采用分子动力学结合嵌入原子多体势,模拟了铂纳米粒子在升温过程中的热稳定性和熔化机制,并利用共近邻分析方法分析了它的微结构演化过程。模拟的结果表明：铂纳米粒子的熔点明显低于体材料的熔点;由于表面层原子的结合力较弱,在升温过程中表面会首先出现预熔;纳米粒子的熔化是从表面层开始的,并随着温度的升高,熔化的表面层会逐渐向内部扩展,最终导致纳米粒子整体转变为液态结构;当温度低于表面预熔温度时,纳米粒子保持良好的晶态结构。     

利用Gupta势研究ConCu13-n (n=1～12)混合团簇的基态结构和熔化行为

利用Gupta多体相互作用势结合遗传算法和分子动力学方法模拟研究了ConCu13-n (n=0～13)单质及混合团簇的基态结构和熔化行为,结果表明：ConCu13-n (n=1~12)混合团簇的基态结构均是在单质Co13、Cu13基态二十面体基础之上的畸变,Co原子先占据中心后占据表面,表面上的Co原子总连接在一起,抱团分布;分析二级差分能和混合能发现Co1Cu12、Co7Cu6具有相对高的稳定性,可视为幻数结构团簇;ConCu13-n (n=1~12)混合团簇的熔点均位于Co13、Cu13单质团簇的熔点之间,且随着Co原子数目的增多总体呈上升趋势,但在n=3时出现反常（熔点降低）,这可归因于Co3Cu10的基态与第一激发态之间的能量差远小于Co2Cu11的相应值。     

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

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

CH2和CH3自由基吸附在Cun (n=1-6)团簇上C-H对称振动模式的软化

利用密度泛函理论中杂化泛函理论方法计算了CH₂和CH₃自由基吸附在Cu_n(n=1～6)团簇上时C?H对称伸缩振动模式的软化性质,结果表明,CH₂在Cun团簇上的吸附要比CH₃的吸附强. 计算得到的C-H键的振动频率与实验上测量的这两个自由基吸附在Cu(111)表面的结果符合得很好,随着团簇尺寸的增加,C-H对称伸缩振动频率的软化(红移)越来越大.     

Cu-Co双金属团簇结构演化及其性质的分子动力学模拟

采用分子动力学结合嵌入原子方法对比研究了Co分布于Cu-Co团簇不同层的结构和性质. 研究表明：Co原子分层掺杂可对团簇的结构转变点和熔点进行诱导控制;分层掺杂的Cu-Co团簇第一相变是一种扩散度较小的由立方八面体转变为二十面体的相变;Co原子易于向低能态团簇的亚表层（111）面偏析, 从而诱导团簇结构紊乱, 造成其熔点差异.     

低能Cu原子簇沉积薄膜的分子动力学模拟

用分子动力学计算机模拟研究了能量为5—20eV/atorn,结构为正二十面体的(Cu)₁₃原子簇在Cu（001）表面的沉积过程.采用紧束缚势同Moliers势的结合描述Cu原子间相互作用通过原子簇-衬底相互作用的“快照”研究沉积的动态过程．结果表明,当入射能量较低时,轰击弛豫后,入射原子簇在衬底表面发生重构,生成很好的外延层,靶没有任何损伤．随着轰击能量的增加,原子簇原子穿入靶的深度增加.当入射能量达到20eV/atom时,原子簇完全穿入靶并开始造成辐照损伤,表面出现空位,靶内产生间     

Ternary alloying effect on the melting of metal clusters

Canonical ensemble Monte Carlo simulations are applied to investigate the melting of the icosahedral 55-atom Ag-Cu-Au clusters. The clusters are modeled by the second-moment approximation of the tight-binding (TB-SMA) many-body potentials. Results show that the introduction of the only Cu atom of the third alloying metal in the bimetallic Ag₄₃Au₁₂ cluster, forming the Ag₄₂Cu₁Au₁₂ cluster, can greatly increase the melting point of the cluster by about 100 K. It is also found that the substitution of the only Cu atom of the third alloying metal in the Ag₁Au₅₄ clusters, forming the Ag₁Cu₁Au₅₃ cluster, can result in an increase of 40 K in the melting point. It can be concluded that the melting points of the bimetallic clusters can be tuned by the third metal impurities doping. In addition, the surface segregation of Ag atoms in the Ag-Cu-Au trimetallic clusters occurs even after melting.     

Molecular Dynamics Simulation of Icosahedral Transformations in Solid Cu-Co Clusters

We study the icosahedral transformations of solid Cu Co clusters with different initial configurations by using molecular dynamics with the embedded atom method. It is found that the formation of symmetric icosahedral cluster is strongly related to the atomic number and initial configuration. The transformation originates from the surface into the interior of the cluster and is a structural change which is rapid and diffusionless. The icosahedral clusters with any composition and configuration, such as core-shell or three-shell duster, can be prepared by the means of solid-solid phase transition in bimetallic clusters.     

Melting behavior of large disordered sodium clusters

The melting-like transition in disordered sodium clusters Na₉₂ and Na₁₄₂ is studied by performing density functional constant-energy molecular dynamics simulations. The orbital-free version of the density functional formalism is used. In Na₁₄₂ the atoms are distributed in two distinct shells (surface and inner shells) and this cluster melts in two steps: the first one, at ≈130 K, is characterized by the development of a high intrashell atomic mobility, and the second, homogeneous melting at ≈270 K, involves diffusive motion of all the atoms across the whole cluster volume. On the contrary, the melting of Na₉₂ proceeds smoothly over a very broad temperature interval, without any abrupt change in the thermal or structural indicators. The occurrence of two steps in the melting transition is suggested to be related to the existence of a grouping of the atoms in radial shells, even if those shells present a large degree of internal disorder. It then appears that a cluster can be considered fully amorphous (totally disordered) only when there are no radial regions of low atomic density separating shells. The isomer of Na₉₂ studied here fulfills this criterion and its thermal behavior can be considered as representative of that expected for fully amorphous clusters. Disordered Na₁₄₂, on the other hand, that has a discernible structure of an inner and a surface shell, should be considered as not fully disordered. The thermal behavior of these two clusters is also compared to that of icosahedral (totally ordered) sodium clusters of the same sizes. Received 5 February 2001 and Received in final form 21 May 2001     

Temperature-accelerated molecular dynamics simulation of the evolution of a low-energy incident Cu<Subscript>3</Subscript> cluster on the Cu(100) surface with a monoatomic step

The molecular dynamics simulation of the normal incidence of a Cu₃ cluster with an energy of 0.2–1 eV/atom on a Cu(100) substrate at an equilibrium temperature of 500–700 K is performed. This substrate contains a rectilinear step, whose height is one atomic layer. After 20-ps relaxation of the atomic cluster on the substrate, its further thermally activated motion on the surface is simulated using the method of temperature- accelerated molecular dynamics. The Cu₃ cluster mainly demonstrates rotational-translational motion. The time interval between two successive atomic transitions of the cluster is found to be reduced as the distance to the step is decreased.     

低温下Au959团簇负载于MgO(100)表面后结构变化的分子动力学研究

采用基于Chen-Mbius反演方法,从金属/金属氧化物界面第一性原理计算的粘结能结果中推导出的Au/MgO原子间相互作用势的正则系综（NVT）分子动力学,模拟了在10 K条件下,Au₉₅₉团簇负载于MgO（100）表面后团簇结构的变化.根据原子对分析技术和对分布函数的分析表明,由于团簇界面处原子间距与载体原子间距相匹配,置于载体上的Au团簇经过一个变形过程后,较其孤立自由表面时的团簇体积变大. 关键词： 团簇 分子动力学 计算机模拟 表面     

熔融Cu55团簇在Cu(010)表面上凝固过程的分子动力学模拟

采用基于嵌入原子方法的分子动力学,模拟了熔融Cu₅₅团簇在Cu衬底（010）表面上以两个不同降温速率降温过程中结构的变化.模拟结果表明,降温速率对团簇结构的变化有很大影响.较快的降温速率使得降温过程中团簇原子具有较低的能量;较慢的降温速率有助于高温时位于衬底内的原子向衬底表面扩散,排列形成面心立方结构. 关键词： 团簇 凝固 分子动力学 表面     

Irregular phenomenon induced by Ag atomic segregation during the heating process of Ag-Pd cluster

This paper studies the melting of icosahedral Ag-Pd bimetallic clusters by using molecular dynamics with the embedded atom method.It finds that the mixed Ag-Pd cluster shows an irregular phenomenon before melting,i.e.,the atomic energy decreases with the increase of temperature.It indicates that the segregation of Ag atoms results in this phenomenon by analysing atomic radius distribution.Since the surface energy of Ag is lower than that of Pd,this leads to the result that the decreased energy by the Ag atomic segregation is larger than the increased energy by the heating.This provides a new method to obtain irregular thermodynamic properties.     

The initial oxidation of molybdenum: IV. Epitaxial growth of oxide on molybdenum

The oxide which grows in low oxygen pressure and at temperatures between 700 and 1000 K on molybdenum is shown to be MoO₂. The epitaxial relationships between the oxide and the metal (100), (110) and (111) surfaces are given. The epitaxial relationships of oxide on the molybdenum (100) and (110) surfaces are geometrically equivalent. The oxide grows on the (111) molybdenum surface with no major oxide plane parallel to the substrate. It is suggested that the epitaxy of MoO₂ on the (111) surface is a consequence of growth on {211} molybdenum facets. The atomic positions in the pairs of interfacial planes found are given. There is little agreement between the positions of ions in the oxide and substrate lattice sites. Only in the postulated case of MoO₂ on {211} Mo facets is a small misfit found.