快凝过程中液态Cu64Zr36合金二十面体团簇遗传与演化跟踪

采用分子动力学方法模拟研究了液态Cu₆₄Zr₃₆合金在冷速50 K/ns下 的快速凝固过程, 并通过双体分布函数、Honeycutt-Andersen (H-A) 键型指数和团簇类型指数对其微结构演变特性进行了分析. 液态与快凝玻璃合金的主要原子组态都是二十面体(12 0 12 0)及其变形结构 (12 8/1551 2/1541 2/1431), 其中比例最高的是Cu芯Cu₈Zr₅基本原子团, 其次是Cu₇Zr₆和Cu₉Zr₄团簇; 并且由这些二十面体基本原子团铰链形成的中程序, 其尺寸分布在液相和固相中分别呈现出13, 19, 25,···和13, 19, 23, 25, 29, 37,···的幻数特征. 团簇的演化与跟踪分析发现: 没有任何团簇能从液态直接遗传到固态合金, 遗传的起始温度出现在T_m–T_g过冷液相区. 二十面体团簇的遗传主要以完全和直接遗传为主, 并且一个明显的增加发生在T_g附近. 在玻璃化转变温度T_g以下, (12 0 12 0) 二十面体比 (12 8/1551 2/1541 2/1431) 变形二十面体具有更高的结构遗传能力, 但仅有少部分在遗传过程中能保持化学成分的恒定. 通过部分遗传, 某些二十面体中程序甚至也能从过冷液体中被遗传到玻璃合金. 关键词： 快速凝固 分子动力学 二十面体团簇 遗传     

液态Ni原子团簇演变的计算机模拟

采用常温常压分子动力学模拟技术,模拟了液态Ni中原子团簇在快速凝固条件下的演变过程,模型采用TB作用势.采用偶分布函数、键对和多面体等结构参数来描述快速凝固条件下团簇种类和数量的变化,并将团簇结构可视化.在2 000 K下,液态Ni中团簇数量较少,由一定数量的1551、1441及1661键对所形成的缺陷二十面体构成;在快速冷却条件下,团簇的数量随温度的降低不断增加,且出现由12个1551键型所形成的完整二十面体团簇,体系最终形成了由二十面体和缺陷二十面体团簇网络所组成的非晶结构.     

Surface structure and solidification morphology of aluminum nanoclusters

Classical molecular dynamics simulation with embedded atom method potential had been performed to investigate the surface structure and solidification morphology of aluminum nanoclusters Al_n (n=256, 604, 1220 and 2048). It is found that Al cluster surfaces are comprised of (1 1 1) and (0 0 1) crystal planes. (1 1 0) crystal plane is not found on Al cluster surfaces in our simulation. On the surfaces of smaller Al clusters (n=256 and 604), (1 1 1) crystal planes are dominant. On larger Al clusters (n=1220 and 2048), (1 1 1) planes are still dominant but (0 0 1) planes cannot be neglected. Atomic density on cluster (1 1 1)/(0 0 1) surface is smaller/larger than the corresponding value on bulk surface. Computational analysis on total surface area and surface energies indicates that the total surface energy of an ideal Al nanocluster has the minimum value when (0 0 1) planes occupy 25% of the total surface area. We predict that a melted Al cluster will be a truncated octahedron after equilibrium solidification.     

Geometric magic numbers of sodium clusters: Interpretation of the melting behaviour

Putative global minima of sodium clusters with up to 380 atoms have been located for two model interatomic potentials in order to identify the structures responsible for the size-dependence of the thermodynamic properties in experiments. Structures based upon the Mackay icosahedra predominate for both potentials, and the magic numbers for the Murrell-Mottram model show excellent agreement with the sizes at which maxima in the latent heat and entropy change at melting have been found in experiment. In particular, the magic numbers at sizes intermediate between the complete Mackay icosahedra are due to unusual twisted icosahedral structures.     

Icosahedral medium-range order formed in Mg70Zn30 metallic glass: a larger-scale molecular dynamics simulation

A larger-scale Mg70Zn30 alloy system including 100000 atoms has been simulated by using the molecular dynamics method to investigate the icosahedral medium-range order (IMRO) formed in the Mg70Zn30 metallic glass. It is found that the simulated pair distribution function of Mg70Zn30 metallic glass is in good agreement with the experimental results. The glass transition temperature Tg is near 450 K under the cooling rate of 1×1012 K/s. The icosahedral local structures play a critical role in the formation of metallic glass, and they are the dominant local configurations in the Mg70Zn30 metallic glass. The IMRO in the Mg70Zn30 metallic glass is characterized by certain types of extended icosahedral clusters combined by intercross-sharing atoms in the form of chains or dendrites. The size distributions of these IMRO clusters present a magic number sequence of 19, 23, 25, 27, 29, 31, 33, 35, 37, 39,. . . , and the magic clusters can be classified into three types according to their compactness. The IMRO clusters grow rapidly in a low-dimensional way with cooling, but this growth is limited near Tg.     

Formation and stabilization of Fe-induced magic clusters on Si(1 1 1)-(7 × 7) template

We demonstrate the growth of Fe-induced magic clusters on Si(1 1 1)-(7 × 7) template by in situ scanning tunneling microscopy (STM). These clusters form near a dimer row at one side of the half-unit cell (HUC); and with three different equivalent orientations. A cluster model comprising three top layer Si atoms bonded to six Fe atoms at the next layer in the 7 × 7 faulted-half template is proposed. The optimized cluster structure determined by first-principles total-energy calculation shows an inward-shifting of the three center Fe atoms. The clusters and the nearby center-adatoms of the next HUCs appear with a significantly reduced height below bias voltages 0.4 V in high resolution empty-state STM images, suggesting an energy gap opening near the Fermi level at these localized cluster and adatom sites. We explain the stabilization of the clusters on the 7 × 7 template using the gain in electronic energy as the driving force for cluster formation.     

稳态Cu-Zr二十面体团簇电子结构的密度泛函研究

采用第一原理对以Cu为心的低能稳态Cu_nZ_(r13-n)(n=6,7,8,9)二十面体团簇的电子结构进行计算,结果表明:同一化学组分下,以Cu为心的Cu-Zr二十面体团簇中出现的同类原子聚集现象可以增强团簇的稳定性,降低费米能级(EF)上的电子数N(EF),这为低能稳态团簇拥有较小的N(EF)提供了深层次的理论解释.进一步的差分电子密度与Mulliken布居分析得知,Cu-Zr二十面体中共价键与离子键共存,成键态与反键态共存,且团簇在形成时壳层Zr与中心Cu原子是电子的提供者,壳层Cu是电子的获得者.该电荷转移方向是金属玻璃中以Cu为心的Cu-Zr二十面体团簇普遍遵循的规律,不随团簇的化学序参数及化学组分的变化而变化.计算的红外振动谱为实验上准确表征不同二十面体原子团提供了一种新的思路.     

液态金属In凝固过程中微观结构转变的模拟研究

采用分子动力学方法对液态金属In的快速凝固过程进行计算机模拟跟踪研究.运用HA键型指 数法和原子成团类型指数法分析了金属原子In的成键类型和形成的原子团簇结构.发现：与通常的液态金属(如Al)相反，随着温度的降低，二十面体及与二十面体相关的1 551 键越来越少；与四方体，六角密堆积相关的1421，1422和1431键数目总和变化很小；而与菱 面体相关的132l，1311，1301和1201的数目却随着温度的降低而显著增加，逐渐占据优势 .最后形成一种新型的以菱面体结构为主、夹杂着立方体(fcc，bcc) 关键词： 液态金属In 微结构转变 团簇结构 分子动力学 计算机模拟     

利用Gupta势结合遗传算法研究ConCu55－n(n=0—55)混合团簇的结构演化及基态能量

采用半经验的Gupta多体势结合遗传算法对Co_nCu_55－n(n=0—55)混合团簇的基态结构和能量进行了研究,发现这些混合团簇的基态结构是在Co₅₅,Cu₅₅单质团簇(Mackay二十面体)的基础之上发生的畸变;从n=0(Cu₅₅)开始,Co原子从中心到表面,从棱到顶点依次、连续替换Cu原子;基态结构与键能较大键的数目及其平均键长有关;Co₁₃Cu₄₂具有最稳定的结构,13个Co原子全部位于团簇内部形成Mackay二十面体对整个团簇的稳定性有显著影响. 关键词： 团簇 结构和能量 Gupta势 遗传算法     

Optical switching in a Ξ system: A comparative study on DROP and EIT

Liquid Lennard-Jones clusters with magic number of atoms N = 55, 147, 309, 561 and 923 were cooled down in Monte Carlo simulations until freezing. Structural properties of the clusters, including the radial dependence of atomic concentration/density and the local regular structure in arrangement of atoms, just before freezing were analysed. Existence of spherical layers in atomic density around the centre of mass of liquid LJ clusters was confirmed. Formation of layers is explained by central net forces acting on every cluster atom and leading to positioning an atom close to the cluster centre of mass. The strong layering in small clusters of N = 55 and 147 affects atomic diffusion in radial and tangential directions inside the cluster, leading to easier movement of atoms on the layer surface. Analysis of radial profiles of four types of structural units detected in liquid clusters reveals that icosahedral units are the most numerous and are located mainly near cluster surface of all clusters and also in the centre of small clusters.     

Deformation of slow liquid and solid clusters upon deposition: A molecular-dynamics study of Al cluster impact on an Al surface

Using molecular-dynamics simulation, we investigate the self-deposition of Al_n clusters (n < 4000) on an Al substrate at velocities below the velocity of sound. Both cold crystalline and hot liquid clusters are studied. We examine the cluster deformation after impact on the surface, which we quantify by its height and base radius. At a given cluster velocity, the shape of deposited crystalline clusters is rather independent of the cluster size; only at small cluster sizes, n ? 40, the clusters are less strongly deformed. With increasing cluster size, liquid clusters are more strongly deformed than crystalline clusters. Faster projectiles become more strongly flattened by the deposition process. Even clusters depositing with vanishing velocity show a finite deformation, which is considerable for smaller clusters. At large cluster speed, clusters penetrate deeper into the (1 0 0) surface than into the (1 1 1) surface and also deform more strongly.     

Local icosahedral order and thermodynamics of simulated amorphous Fe

Local icosahedral order and thermodynamics of amorphous Fe have been analyzed in detail for models containing 3000 atoms, which were obtained by the molecular dynamics (MD) method. Models were obtained by cooling from the melt. Local order in models has been analyzed by using the technique proposed by Honeycutt and Andersen; we found an existence of icosahedral order in the system. Moreover, structural properties of models were also studied via radial distribution function (RDF), static structure factor, mean atomic distances, coordination number and bond-angle distributions. Glass transition temperature, heat capacity and potential energy of the system were found in addition to the evolution of structure and mean-squared displacement (MSD) of atoms upon cooling from the melt toward the glassy state. We found the glass transition temperature of simulated liquid Fe via temperature dependence of potential energy and it is close to that observed previously in the literature, i.e. T_g≈1070 K. Calculations showed that structural properties of amorphous Fe models with the Pak-Doyama interatomic potential agreed well with the experimental data.     

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.     

Experimental and theoretical studies of interactions between Si<Subscript>7</Subscript> clusters

The possibility of using magic Si₇ clusters to form a cluster material was studied experimentally and theoretically. In experiments Si₇ clusters were deposited on carbon surfaces, and the electronic structure and chemical properties of the deposited clusters were measured using X-ray photoelectron spectroscopy (XPS). A non bulk-like electronic structure of Si₇ was found in the Si 2p core level spectra. Si₇ is suggested to form a more stable structure than the non-magic Si₈ cluster and Si atoms upon deposition on carbon surfaces. Theoretically it was possible to study the interaction between the clusters without the effect of a surface. Density functional theory (DFT) calculations of potential curves of two free Si₇ clusters approaching each other in various orientations hint at the formation of cluster materials rather than the fusion of clusters forming bulk-like structures.     

Scanning tunneling microscopy study of growth of Pt nanoclusters on thin film Al2O3/NiAl(1 0 0)

The growth of Pt nanoclusters on thin film Al₂O₃ grown on NiAl(1 0 0) was studied by using scanning tunneling microscopy (STM). The samples were prepared by vapor depositing various amounts of Pt onto the Al₂O₃/NiAl(1 0 0) at different substrate temperatures in ultra high vacuum (UHV). The STM images show that sizeable Pt nanoclusters grow solely on crystalline Al₂O₃ surface. These Pt clusters appear to be randomly distributed and only a few form evident alignment patterns, contrasting with Co clusters that are highly aligned on the crystalline Al₂O₃. The size distributions of these Pt clusters are rather broader than those of the Co clusters on the same surface and the sizes are evidently smaller. With increasing coverage or deposition temperature, the number of larger clusters is enhanced, while the size of the majority number of the clusters remains around the same (0.4 nm as height and 2.25 nm as diameter), which differs drastically from the Pt clusters on γ-Al₂O₃/NiAl(1 1 0) observed earlier. These Pt cluster growth features are mostly attributed to smaller diffusion length and ease to form stable nucleus, arising from strong Pt-Pt and Pt-oxide interactions and the peculiar protrusion structures on the ordered Al₂O₃/NiAl(1 0 0). The thermal stability of Pt nanoclusters was also examined. The cluster density decreased monotonically with annealing temperature up to 1000 K at the expense of smaller clusters but coalescence is not observed.     

Pressure effect on the structural transition of liquid Au

We present a molecular dynamics simulation within the framework of empirical tight-binding potential on the liquid structure of Au under different pressures during the rapid cooling process. The pair correlation function (PCF) and the pair analysis (PA) technique are used to reveal the structural characteristics of liquid Au under normal and high pressures. The split of the second peak of g(r) is associated with the glass transition at the cooling rate of 2.06×10¹³ K/s under different pressures. The nearest distance is shortened under high pressures. High-pressure is in favor of FCC-type and BCC-type atomic clusters. The number of icosahedral structures increases in the high temperature region as pressure increase, while high-pressure is not preferable for icosahedra at low temperature. With the increase of the pressure, the possibility that the system forms crystalline structure is enhanced. The influence of a higher pressure on the microstructure of liquid metal is much clearer than that of a lower pressure.     

Molecular dynamics simulation of the solidification of liquid gold nanowires

The solidification behavior of liquid gold nanowires with about 1.84 nm in diameter has been studied by using molecular dynamics simulation with an embedded atom potential. It is found the cooling rate has great effect on the final structure of the gold nanowires during solidification from liquid. With the decrease of cooling rates, the final structure of the gold nanowires varies from amorphous to crystalline via helical multi-shelled structure. The face-centered cubic structure of the gold nanowires is proven energetically the most stable form.     

Scanning tunneling microscopy and abinitio studies of precursor states of Ga-induced cluster on Si(0 0 1) surface

The growth processes and structures of Ga layers formed on a Si(0 0 1) surface have been studied by scanning tunneling microscopy (STM) and abinitio calculation. The precursor states of the Ga clusters that compose the Si(0 0 1) 8 × n-Ga (n = 4, 5, 6) structures are observed in addition to the 2 × 2- and 4 × 2-Ga structures at a Ga coverage of 0.55 ML at 300 °C. There are two types of precursor clusters whose protrusions are observed as different shapes in filled-state STM images. We compare the observed STM images with the simulated ones to identify the possible structural models. From the results, we determine the structure of each precursor cluster. On the basis of the results, the formation processes of the cluster are discussed.     

二十面体团簇的遗传:一个与快凝Cu56Zr44合金玻璃形成能力有关的动力学参数

采用分子动力学方法模拟研究了液态Cu₅₆Zr₄₄合金在不同冷速γ与压力P下的快速凝固过程, 并通过基于Honeycutt-Andersen键型指数的扩展团簇类型指数法对其微结构演变特性进行了分析. 结果表明: 快凝玻璃合金的局域原子组态主要是(12 12/1551)规则二十面体、以及 (12 8/1551 2/1541 2/1431)与(12 2/1441 8/1551 2/1661) 缺陷二十面体. 通过原子轨迹的逆向跟踪分析发现: 从过冷液体中遗传下来的二十面体对快凝合金的玻璃形成能力(GFA)具有重要影响, 不仅其可遗传分数F_i =N_{300 K←Tg}ⁱ/N_Tg 与GFA密切相关, 而且其遗传起始温度(T_onset)与合金约化玻璃转变温度T_rg = T_g/T_m也存在很好的对应关系.     

Density functional study of TaSin (n = 1-3, 12) clusters adsorbed to graphene surface

A plane-wave density functional theory (DFT) calculations have been performed to investigate structural and electronic properties of TaSi_n (n = 1-3, 12) clusters supported by graphene surface. The resulting adsorption structures are described and discussed in terms of stability, bonding, and electron transfer between the cluster and the graphene. The TaSi_n clusters on graphene surface favor their free-standing ground-state structures. Especially in the cases of the linear TaSi₂ and the planar TaSi₃, the graphene surface may catalyze the transition of the TaSi_n clusters from an isomer of lower dimensionality into the ground-state structure. The adsorption site and configuration of TaSi_n on graphene surface are dominated by the interaction between Ta atom and graphene. Ta atom prefers to adsorb on the hollow site of graphene, and Si atoms tend to locate on the bridge site. Further, the electron transfer is found to proceed from the cluster to the surface for n = 1 and 2, while its direction reverses as n > 2. For the case of TaSi, chemisorption is shown to prevail over physisorption as the dominant mode of surface-adsorbate interaction by charge density analysis.