Aln(n=13–32)团簇熔化行为的分子动力学模拟研究

基于半经验的Gupta原子间多体相互作用势, 采用分子动力学方法并结合模拟退火及淬火技术, 系统研究了小尺寸铝团簇Al_n (n=13–32)的熔化行为. 模拟结果表明: 除个别尺寸(Al₁₃ 和Al₁₉) 外, 团簇熔化过程比热曲线普遍呈现杂乱无规(无明显单峰)现象, 这与实验观测小Al团簇比热普遍无规的结果完全一致. 通过分析不同温度点上团簇淬火构型的势能分布图给出了小Al团簇比热呈现无规或有规现象的成因. 对于比热无规团簇, 可以利用原子等价指数判断给出团簇熔点, 所得团簇熔点随团簇尺寸的变化趋势与实验观测结果完全一致. 关键词： Gupta势 团簇 分子动力学 熔化     

Structures, stabilities and ordering in Ni-Al nanoalloy clusters

A detailed study is made of the structures and relative stabilities of nickel and aluminium clusters and nickel-aluminium "nanoalloy" clusters, with up to 55 atoms, modelled by the many-body Gupta potential. Random search and genetic algorithms are used to find the lowest energy isomers (both geometrical and permutational). For the pure Al and Ni clusters, the lowest energy structures are identical for most nuclearities but different structures are found for clusters with 15-18 atoms -- the Al clusters having non-icosahedral structures. For these nuclearities, we investigate the effect of doping Al atoms into pure Ni clusters and vice versa, finding that the replacement of a single atom by a dopant atom is often sufficient to change the cluster geometry. The lowest energy isomers of nanoalloy clusters with the approximate composition "Ni₃Al" generally have structures based on icosahedral packing, though truncated octahedral (fcc packing) motifs are also observed. In agreement with previous studies, the atom ordering in the mixed Ni-Al clusters is found to depend on the maximization of the number of Ni-Al interactions, the minimization of the cluster surface energy, and atom size effects.Received: 17 January 2003, Published online: 17 July 2003PACS: 36.40.Mr Spectroscopy and geometrical structure of clusters - 61.46.+w Nanoscale materials: clusters, nanoparticles, nanotubes, and nanocrystals - 61.66.Dk AlloysM.S. Bailey: Present address: Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, NY 14853-1301, USA.     

Density functional study of Al-doped Au clusters

The equilibrium geometries and electronic properties of Au_nAl, up to n=13, have been systematically investigated using the density functional theory. The results show that, for the Au_nAl clusters, two patterns are identified. Pattern one (n=2, 3, 8), the lowest-energy geometries prefer two-dimensional structures. Pattern two (n=4-7, 9--13), the lowest-energy geometries prefer three-dimensional structures. According to the analysis of the binding energy and the fragmentation energy, Au_nAl clusters with odd n are found to be more stable than those with even n. The same trend of alternation can be illuminated according to the computational results in the HOMO--LUMO gap, the ionization potential, and the electron affinities. The Al atom not only changes the structures of pure gold clusters, but also enhances their stabilities. NBO analysis indicates 6s orbital of Au atom hybridizes with 3p orbital of Al atom.     

Mag netic properties of $$\hbox{Ni}_{13- n}\hbox{Al}_ n$$ clusters with n = 0–13

We report the magnetic properties of small Ni_13-nAl_n\hbox{Ni}_{13-n}\hbox{Al}_n clusters with n = 0–13 calculated in the framework of density functional theory. The cluster magnetic moment decreases with the sequential substitution of Ni by Al atoms, which can be attributed to a greater degree of hybridization that forces the pairing of the electrons in the molecular orbitals of Ni and Al. For Ni₇Al₆, the complete quenching of the cluster magnetic moment appears to be due to the antiferromagnetic alignment of atomic spins as revealed by the spin density plots.     

Molecular dynamics simulations of the melting of Al–Ni nanowires

Molecular dynamics (MD) simulations were performed to investigate the influence of nickel (Ni) composition and nanowire thickness on the thermal properties of Al-x%Ni (at%) nanowires using the embedded atom model (EAM) potential. The melting of the nanowire was characterised by studying the temperature dependence of the cohesive energy and mean square displacement. The effect of the nanowire thickness on the cohesive energy, melting temperature, heat capacity as well as latent heat was studied in canonical ensemble. Moreover, the crystal stability of Al, Al-20%Ni, Al-40%Ni, Al-60%Ni, Al-80%Ni, Al₃Ni, Ni₃Al and Ni nanowires was studied at different temperatures using mean square displacement and cohesive energy.     

密度泛函方法研究NiSin(n=1~13)团簇

基于第一性原理,利用密度泛函理论中的广义梯度近似(GGA)系统研究了NiSin(n=1～13)团簇,在充分考虑自旋多重度的基础上讨论了这些团簇的生长行为,电子性质及其磁性,结果表明:NiSin 1的基态结构是在NiSin的基态结构上带帽一个Si原子而得到;随着团簇尺寸的增大,Ni原子逐渐从吸附在sin团簇的表面位置移动到Sin团簇笼内;掺杂Ni原子提高了硅团簇的稳定性;NiSi10团簇的稳定性在所有团簇中是最高的;电子总是从si向Ni转移,Ni原子所带的电荷数不仅与Ni原子的配位数有关,还与Nisin团簇的基态结构密切相关;n=1～2时,团簇的自旋总磁矩为2 μB,当n≥3时,团簇的磁性消失,这可能与Ni原子内部较强的sp-d杂化以及si原子内部的s-p杂化有关.     

Diffusion and growth of nickel,iron and magnesium adatoms on the aluminum truncated octahedron: A molecular dynamics simulation

The structure of nickel (Ni), iron (Fe), and magnesium (Mg) adatoms on the aluminum (Al) truncated octahedron is studied using molecular dynamics and the analytic embedded atom method. First, the energy barriers of several typical diffusion processes of Ni, Fe, and Mg adatoms on the Al truncated octahedral cluster were calculated using the nudged elastic band method. The calculated energy barriers were found to be related to the surface energy and atomic radius of the adatom and substrate atom. The result shows that the incorporation of Ni and Fe atoms into Al core easily occurs, and the Mg atom should segregate at the surface of the Al cluster. Thus, the growth of Ni, Fe and Mg on the Al truncated octahedron with 1289 atoms was simulated at several temperatures. In the Ni–Al and Fe–Al cases, the core-shell structure was not obtained. For the Mg–Al system, a good Mg shell on the Al core was found at lower temperatures, and an almost perfect truncated octahedron with more Al shells emerged with an increase in temperature.     

AlnCl(n=2-14)团簇结构与电子性质的理论研究

运用卡里普索(CALYPSO)结构预测方法,在杂化密度泛函B3LYP/6-311G+(d)基组水平上,对AlnCl(n=2-14)团簇的几何结构与电子性质进行优化计算,并讨论了团簇的平均结合能、能隙、二阶能量差分、电离能、亲和能以及电子自然布局和极化率。研究结果表明：AlnCl(n=2-14)团簇的基态构型由简单平面几何结构向立体结构演化,形成Cl原子戴帽Aln-1Cl团簇结构;Cl原子的掺杂增大了Aln团簇的平均结合能;二阶能量差分、能隙、电离能、亲和能的变化表明Al7Cl是幻数团簇结构;团簇中的电荷总是由Al原子向Cl原子转移,原子之间的成键作用随着团簇尺寸的增大而增强。     

Al_nCl(n=2-14)团簇结构与电子性质的理论研究

运用卡里普索(CALYPSO)结构预测方法,在杂化密度泛函B3LYP/6-311G+(d)基组水平上,对Al_nCl(n=2-14)团簇的几何结构与电子性质进行优化计算,并讨论了团簇的平均结合能、能隙、二阶能量差分、电离能、亲和能以及电子自然布居和极化率.研究结果表明:Al_nCl(n=2-14)团簇的基态构型由简单平面几何结构向立体结构演化,形成Cl原子戴帽Al_n-1Cl团簇结构;Cl原子的掺杂增大了Al_n团簇的平均结合能;二阶能量差分、能隙、电离能、亲和能的变化表明Al_7Cl是幻数团簇结构;团簇中的电荷总是由Al_原子向Cl原子转移,原子之间的成键作用随着团簇尺寸的增大而增强.     

低能Cu13团簇沉积薄膜的分子动力学模拟研究

利用分子动力学模拟方法对Cu₁₃团簇在Fe（001）表面上沉积薄膜进行了研究,分析了不同沉积条件对薄膜生长模式的影响,对比分析了不同沉积条件下表面粗糙度、缺陷分布和外延度等薄膜性质的差异。Cu₁₃团簇的初始沉积能量范围为0.1~10.0 eV/atom,沉积率为1.0 clusters/ps,衬底温度分别为300,700和1 000 K。模拟结果表明:团簇初始沉积能量主要影响薄膜生长模式,当初始沉积能量为7.5 eV/atom的Cu₁₃团簇沉积到温度为300 K的Fe（001）表面时,可形成表面光滑、内部缺陷少和较好外延度的高质量Cu薄膜。     

Aln(n=2～7)团簇的结构和能级分布

采用密度泛函B3LYP的方法研究了小原子团簇Al2～7的几何结构和能级分布,分析了随团簇原子数的增加,团簇的几何结构和费米能级的变化情况.研究结果表明:Al2～7的团簇的几何结构在5个原子以前为平面结构,而从六个原子开始为空间立体的稳定结构.电子壳层结构表明,在铝团簇中没有出现非常明显的象碱金属那样的稳定幻数结构.在Al2～Al7团簇中,能级结构呈现明显的分立特征,费米能级随着原子个数的增加而减小,到Al7时又有所增加,且团簇的能量间隙最小.     

Deposition of Ni<Subscript> 13</Subscript> and Cu<Subscript> 13</Subscript> clusters on Ni(111) and Cu(111) surfaces

The soft deposition of Ni₁₃ and Cu₁₃ clusters on Ni(111) and Cu(111) surfaces is studied by means of constant-energy molecular-dynamics simulations. The atomic interactions are described by the Embedded Atom Method. It is shown that the shape of the nickel clusters deposited on Cu(111) surfaces remains rather intact, while the copper clusters impacting on Ni(111) surfaces collapse forming double and triple layered products. Furthermore, it is found that for an impact energy of 0.5 eV/atom the structures of all investigated clusters show the lowest similarity to the original structures, except for the case of nickel clusters deposited on a Cu(111) surface. Finally, it is demonstrated that when cluster and substrate are of different materials, it is possible to control whether the deposition results in largely intact clusters on the substrate or in a spreading of the clusters. This separation into hard and soft clusters can be related to the relative cohesive energy of the crystalline materials.     

NiMgn(n=1—12)团簇的第一性原理研究

采用基于密度泛函理论(DFT)中的广义梯度近似 (GGA),在考虑自旋多重度的情况下,对NiMg_n(n=1—12)团簇进行了构型优化,频率分析和电子性质计算.结果表明：n=1,2时,体系的基态为自旋三重态,n≥3时,为单重态;Ni原子掺杂使主团簇结构发生了明显变化. n≤8时,三角双锥,四角双锥结构主导着NiMg_n基态团簇的生长行为; n在9—12之间时,主团簇Mg_n+1(n=1—12)的基于三棱柱构型的基态演化行为发生了一定程度的改变;n≥6时,Ni原子陷入了主团簇内部;掺杂使体系的平均结合能增大,能隙减小;n=4,6,10是团簇的幻数;不同尺寸团簇的s, p, d轨道杂化中,Ni原子3d, 4p成分所起作用不同; NiMg₆基态结构具有很高的对称性(O_h),很好的稳定性和化学活性,能隙仅为0.25eV. 关键词： n团簇')" href="#">NiMg_n团簇 几何结构 稳定性 化学活性     

金团簇碰撞并合过程的分子动力学模拟研究

采用微正则系综理论和嵌入原子模型计算了原子数为16-152之间7个Aun团簇碰撞合并过程。结果表明：碰撞能量大于0.3eV时,团簇处于熔融态,碰撞能量是合并过程的主导因素,碰撞能量小于0.3eV时,团簇处于凝固态,碰撞能量和结合能共同影响合并过程;碰撞能量小于0.1eV时,在结合能作用下,团簇结构随碰撞能量增加发生变化,其中壳层结构势能较低,稳定性较好。     

小尺寸铝纳米团簇的相变行为

采用分子动力学方法模拟了半径从0.3–1.3 nm变化的小尺寸铝纳米团簇的熔化、凝固行为. 基于势能-温度曲线、热容-温度曲线分析, 获得了熔点、凝固点与尺寸的依变关系, 并利用表面能理论、小尺寸效应开展了现象分析.研究表明, 铝团簇原子数小于80时, 熔点和凝固点的尺寸依赖性出现无规律的异常变化; 而大于该原子数, 熔、凝固点则随着团簇尺寸的减小而单调下降; 当原子数为27时, 团簇熔点高于块材熔点近40 K. 同时, 铝纳米团簇呈现出凝固滞后现象, 即凝固点低于熔点. 关键词： 纳米团簇 熔点 凝固点 分子动力学     

Cluster formulae for metallic glasses derived from devitrification phases

Bulk metallic glasses are known to have a composition formula [cluster](glue atom)_1,3 within the framework of the cluster-plus-glue-atom model. The key issue in applying the cluster formula is the determination of the right clusters and glue atoms. As examples, alloy phases in the glass-forming systems Al–Ni–Zr and B–Co–Si are analysed from the viewpoint of nearest coordination polyhedral clusters. These alloy phases are described by [effective cluster](glue atom) _x , where the effective cluster refers to true cluster composition after taking account of cluster-sharing in the phase structure. For each alloy phase, a principal cluster can be identified that features the local short-range order of that phase. It is pointed out that the principal clusters can express compositions with high glass-forming abilities, as verified by our experiments in Al–Ni–Zr and B–Co–Si–Ta.     

不同热浴下Co团簇熔化行为的分子动力学模拟

运用分子动力学方法,采用Berendsen热浴和Nose-Hoover热浴分别研究了Con (n=13,55,147)团簇的熔化特性,模型采用Gupta相互作用势.模拟结果表明：两种热浴对钴团簇熔点及预熔化区间给出了基本一致的描述.所研究团簇体系在给定温度下长时间内各Co团簇中单个原子的速率(速度)分布与麦克斯韦速率(速度)分布曲线符合很好.     

氧原子在X@Al12 (X=C，Si，P+ )团簇吸附性质的研究

利用密度泛函理论了对X@Al12 (X=C, Si, P+)团簇吸附氧原子特性进行了研究,分别分析了氧原子在三个吸附位置(顶位置、桥位置和空位置)的氧原子吸附能,HOMO-LOMO能隙和电荷转移等性质。结果表明氧原子更倾向于吸附在空位位置。C@Al12和Si@Al12, P+@Al12 在对氧的吸附性质上比较接近。氧在X@Al12(X= C, Si, P+)的吸附性质与Al13−有显著不同。在Al13−O中O是倾向于以桥位形式吸附在Al13−上的。X@Al12(X=C, Si, P+)和Al13−的电子结构有着明显的差异,这导致了O在这些团簇表面不同的吸附形式,表明掺杂可以有效的改变团簇的性质。     

不同热浴下Co团簇熔化行为的分子动力学模拟

运用分子动力学方法,采用Berendsen热浴和Nose-Hoover热浴分别研究了Co_n(n=13,55,147)团簇的熔化特性,模型采用Gupta相互作用势.模拟结果表明:两种热浴对钴团簇熔点及预熔化区间给出了基本一致的描述.所研究团簇体系在给定温度下长时间内各Co团簇中单个原子的速率(速度)分布与麦克斯韦速率(速度)分布曲线符合很好.     

氧原子在X@Al12 (X=C，Si，P+ )团簇吸附性质的研究

利用密度泛函理论了对X@Al12 (X=C, Si, P+)团簇吸附氧原子特性进行了研究,分别分析了氧原子在三个吸附位置(顶位置、桥位置和空位置)的氧原子吸附能,HOMO-LOMO能隙和电荷转移等性质。结果表明氧原子更倾向于吸附在空位位置。C@Al12和Si@Al12, P+@Al12 在对氧的吸附性质上比较接近。氧在X@Al12(X= C, Si, P+)的吸附性质与Al13−有显著不同。在Al13−O中O是倾向于以桥位形式吸附在Al13−上的。X@Al12(X=C, Si, P+)和Al13−的电子结构有着明显的差异,这导致了O在这些团簇表面不同的吸附形式,表明掺杂可以有效的改变团簇的性质。