Magnetic anisotropy in icosahedral cobalt clusters

Magnetic anisotropy has been measured in multiply twinned, icosahedral cobalt clusters. It is found that the low-temperature magnetization of deposited cluster layers is well defined with the Stoner–Wohlfarth model by averaging over clusters with a range of anisotropy energy. Anisotropy energy calculation based on Néel's pair model shows that the icosahedral structure and the layer-by-layer growth of the clusters induce oscillations of the magnetic anisotropy as a function of the filling of the outer surface of the particle. The magnetization measurement at room temperature indicates a weakly correlated cluster glass, as deduced from the approach to saturation that is well described with 2D random anisotropy model.     

稳态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二十面体团簇普遍遵循的规律,不随团簇的化学序参数及化学组分的变化而变化.计算的红外振动谱为实验上准确表征不同二十面体原子团提供了一种新的思路.     

Melting Behaviour of Core-Shell Structured Ag-Rh Bimetallic Clusters

The melting behaviour of four typical core-shell structured 309-atom Ag-Rh bimetallic clusters, with decahedral and icosahedral geometric configurations, is investigated by using molecular dynamics simulation, based on the Sutton-Chen potential. The initial atomic configurations are obtained from semi-grand canonical ensemble Monte Carlo simulations. It is found that the melting point temperature Tm increases with the mole fraction of Rh in the bimetallic clusters, and Tm of Ag-Rh icosahedral clusters is higher than those of Ag-Rh decahedral clusters with the same Rh mole fraction. It is also found that the Ag atoms lie on the surface of Ag-Rh bimetallic clusters even after melting.     

Unsupported lead clusters and electron diffraction

A beam of Pb clusters is produced with the inert gas aggregation method and probed by electron diffraction. Analysis of the diffraction patterns indicates that average cluster size can vary between 3 and 7 nm, according to nucleation conditions. The diffraction patterns from beams with larger average cluster size are very similar to patterns calculated from model decahedron clusters, while those for smaller cluster size do not appear to have simple geometrical face-centred cubic, decahedral, or icosahedral structure. Received 30 November 2000     

Surface tensions and stress tensors of liquid and solid clusters by molecular dynamics

Surface tension and pressure (stress) tensors of Lennard-Jones clusters, in the size range 200 ~ 2700 atoms/cluster, formed from evaporating liquid droplets were calculated in a Molecular Dynamics simulation. Icosahedral clusters have a much larger surface tension than decahedral, fcc, and hcp ones, meanwhile asymmetric icosahedral clusters have a lower surface tension. Fcc and hcp clusters have a very small surface tension. Decahedral clusters have a surface tension closer to that of fcc and hcp ones than to that of icosahedral ones, though both icosahedral and decahedral structures have five fold symmetry axis. Binary component clusters have a higher surface tension than single component ones.     

二十面体Sc13，Sc13¬+1，Sc13¬-1团簇的稳定性与电子结构研究

采用密度泛函理论框架下的广义梯度近似(DFT/GGA),对Sc13团簇进行了几何结构优化,得到13原子钪团簇的基态结构为正二十面体(Ih),在此基础上对二十面体(Ih)Sc13, Sc13¬+1和Sc13-1团簇的稳定性、电子结构磁矩进行研究.结果发现Sc13,Sc13¬+1和Sc13¬-1团簇都在Ih结构时最稳定,该尺寸团簇的稳定性主要由二十面体密堆积构型决定;带电能使团簇的结构稍稍收缩从而使团簇的稳定性有所增强;团簇的键长和对称性对团簇的磁矩有很明显的影响.     

Molecular dynamics simulations of the freezing of gold nanoparticles

A set of molten gold clusters, each with 1157 gold atoms, was studied by molecular dynamics simulations as the clusters underwent freezing at three different temperatures. Most of the clusters attained an icosahedral structure upon freezing, a structure found to be stable to mild annealing. Other structures observed were imperfect truncated decahedral, truncated octahedral and hexagonal close packed structures. The role of kinetics in the process of cluster solidification is discussed. Received 6 November 2000     

Investigation of the charge distribution in small cluster ions Ar<Stack><Subscript>13</Subscript><Superscript>+</Superscript></Stack> and Ar<Stack><Subscript>19</Subscript><Superscript>+</Superscript></Stack>

The results of ab initio studies of the atomic and charge structure of small clusters and cluster ions formed by 13 and 19 argon atoms are reported. It was found that the icosahedral atomic structure is energetically the most favorable for such clusters. The calculations demonstrate that when a single electron is removed from a cluster, the excess positive charge is distributed primarily over the surface of the formed cluster ion.     

Competition between face-centered cubic and icosahedral cluster structures

The binding energy of atoms in icosahedral and face-centered-cubic clusters is calculated numerically for pairwise Morse-potential interactions between atoms and for clusters containing from 561 to 923 atoms, which corresponds to gradual filling of the sixth layer of the icosahedral cluster. Perturbation theory is used to calculate the cluster binding energy, in which the small parameter is the ratio of the interaction energy between non-nearest neighbor atoms to the interaction energy between nearest neighbors. Values of the Morse interaction potential parameter are found for which the energies of clusters with different structures coincide. Under the conditions used in these computations, the strain energy of a cluster can be neglected. Although the contribution of the interaction energy between non-nearest neighbors to the total cluster energy is small, it turns out to be important in finding the level crossing parameter. Zh. éksp. Teor. Fiz. 112, 1082–1090 (September 1997)     

Sc12X团簇(X=B,C,N,Al,Si,P)的电子结构和稳定性

基于第一性原理,在密度泛函理论框架下,用广义梯度近似(GGA)研究二十面体Sc12X(X=B、C、N、Al、Si、P)中性和荷电团簇的电子结构和稳定性,系统计算了它们的基态束缚能(BE)、原子间平衡间距、最高占据轨道(HOMO)与最低占据轨道(LUMO)之间的能隙、局域电荷以及HOMO电子构型.研究表明,用C、Si原子或荷电后的B、N、Al、P离子分别替代团簇Sc13中心原子可以使其成为稳定的结构.Sc12X团簇束缚能改变的原因在于掺杂改变了中心原子或离子与表面原子的轨道杂化.     

Three-dimensional simulation on explosions of hydrogen atomic clusters irradiated by an intense femtosecond laser pulse

Using classic particle dynamics simulations, the interaction process between an intense femtosecond laser pulse and icosahedral hydrogen atomic clusters H_{13}, H_{55} and H_{147} has been studied. It is revealed that with increasing number of atoms in the cluster, the kinetic energy of ions generated in the Coulomb explosion of the ionized hydrogen clusters increases. The expansion process of the clusters after laser irradiation has also been examined, showing that the expansion scale decreases with increasing cluster size.     

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

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

Novel metal-encapsulated caged clusters of silicon and germanium

We report the recent findings of metal (M) encapsulated clusters of silicon from computer experiments based on ab initio total energy calculations and a cage shrinkage and atom removal approach. Our results show that using a guest atom, it is possible to wrap silicon in fullerenelike (f) structures, as sp² bonding is not favorable to produce empty cages unlike for carbon. Transition M atoms have a strong bonding with the silicon cage that are responsible for the compact structures. The size and structure of the cage change from 14 to 20 Si atoms depending upon the size and valence of the M atom. Fewer Si atoms lead to relatively open structures. We find cubic, f, Frank-Kasper (FK) polyheral type, decahedral, icosahedral and hexagonal structures for M@Si_n with n = 12-16 and several different M atoms. The magic behavior of 15 and 16 atom Si cages is in agreement with experiments. The FK polyhedral cluster, M@Si₁₆ has an exceptionally large density functional gap of about 2.35 eV calculated within the generalized gradient approximation. It is likely to give rise to visible luminescence in these clusters. The cluster-cluster interaction is weak that makes such clusters attractive for cluster assembled materials. Further studies to stabilize Si₂₀ cage with M = Zr, Ba, Sr, and Pb show that in all cases there is a distortion of the f cage. Similar studies on M encapsulated germanium clusters show FK polyhedral and decahedral isomers to be more favorable. Also perfect icosahedral M@Ge₁₂ and M@Sn₁₂ clusters have been obtained with large gaps by doping with divalent M atoms. Recent results of the H interaction with these clusters, hydrogenated silicon fullerenes as well as assemblies of clusters such as nanowires and nanotubes are briefly presented.     

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

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

Composition-dependent electrical resistivity in an Al-Re-Si 1 /1-cubic approximant phase: an indication of electron confinement in clusters

We report the synthesis of alpha-AlReSi and show that it is a 1/1-cubic approximant phase of the icosahedral quasicrystal with a = 12.9 A. The trend of the resistivity of the new approximant phase shows a nonmetallic character, similar to those seen in the stable icosahedral phases. The resistivity depends sensitively on the Re concentration and the nonmetallic transport is observed only at the Re concentration close to 17.4 at. %, where the transition metal sites in the icosahedral cluster are exclusively occupied by Re atoms. In view of a recent ab initio calculation, the present result suggests strongly the formation of the virtual bound states, or confinement of electrons, in the icosahedral clusters of transition metal atoms.     

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.     

Independent variables in foam clusters

We find, by counting the degrees of freedom of two-dimensional bubble clusters (finite or periodic) of given topology and bubble areas, that the Plateau laws determine a unique configuration of a finite free cluster, but allow an infinite number of configurations of a periodic cluster. Each of these configurations is associated with a particular strain (stress) state of the cluster; there is in general one unstrained configuration, which corresponds to the minimum of the (surface) energy. Configurations of given topology that satisfy Plateau's laws may only exist in certain ranges of bubble area ratios and/or strains. Received 31 May 2001 and Received in final form 12 September 2001     

Density functional investigations for geometric and electronic properties of In4M and In12M （M = C,Si, In） clusters

First-principle calculations are performed to study geometric and electronic properties of both neutral and anionic In4M and In12M （M = C, Si, In） clusters. In4C and In4Si are found to be tetrahedral molecules. The icosahedral structure is found to be unfavourable for In12M. The most stable structure for In12C is a distorted buckled biplanar structure while for In12Si it is of an In-cage with the Si located in the centre. Charge effect on the structure of In12M is discussed. In4C has a significantly large binding energy and an energy gap between the highest-occupied molecularorbital level and the lowest unoccupied molecular-orbital level, a low electron affinity, and a high ionization potential, which are the characters of a magic cluster, enriching the family of doped-group-IIIA metal clusters for cluster-assembled materials.     

Molecular dynamics simulation about porous thin-film growth in secondary deposition

The thin film growth has been confirmed to be assembled by an enormous number of clusters in experiments of CVD. Sequence of clusters’ depositions proceeds to form the thin film at short time as gas fluids through surface of substrate. In order to grow condensed thin film using series of cluster deposition, the effect of initial velocity, substrate temperature and density of clusters on property of deposited thin film, especially appearance of nanoscale pores inside thin film must be investigated. In this simulation, three different cluster sizes of 203, 653, 1563 atoms with different velocities (0, 10, 100, 1000 and 3000 m/s) were deposited on a Cu(0 0 1) substrate whose temperatures were set between 300 and 1000 K. Four clusters and one cluster were used in primary deposition and secondary deposition, respectively. We have clarified that adhesion between clusters and substrate is greatly influenced by initial velocity. As a result, the exfoliation pattern of deposited thin film is dependent on initial velocity and different between them. One borderline dividing whole region into porous region and nonporous region are obtained to show the effect of growth conditions on appearance of nanoscale pores inside thin film. Moreover, we have also shown that the likelihood of porous thin film is dependent on the point of impact of a cluster relative to previously deposited clusters.