Molecular Dynamics Investigation of Differences in Melting Behaviors of Cu57 and Cu58 Clusters

应用基于嵌入原子方法的分子动力学计算研究了Cu57和Cu58团簇在升温过程中变为熔体的结构演化过程．两个团簇在熔化时表现出不同的结构变化行为,进而影响到它们能量变化的差异．在升温时,团簇不同区域的原子局域结构变化由原子密度分布函数确定．模拟表明,即使对于这两个仅相差一个原子的小尺寸团簇,结构变化也敏感于团簇的尺寸．     

Molecular dynamics study of structural changes in freezing clusters

Structural changes of Cu₅₁–Cu₅₆ clusters during freezing have been studied by employing molecular dynamics(MD) simulations. With the growth of the crystallization phase, local patterns in different regions determined by the atom density profiles present different characteristic for these studied clusters. While the atom movements in the even-number clusters result in the formation of initial two-shell, then three-shell, and finally four-shell structures, the atoms in the odd-number clusters, with one exception of the Cu₅₁ cluster, first move into the central region and the outermost parts, then follow the movement from the inner parts into the outer parts. The crystallization temperature of these clusters has an oscillatory increase behavior for odd- and even- numbered copper atoms with increasing cluster size.     

Local structure changes of 54-, 55-, 56-atom copper clusters on heating

The local structure changes among three types of Cu clusters containing 54-56 atoms at elevated temperatures have been studied by employing molecular-dynamics simulations. The simulations show sensitivities of structural changes to different structures of the three clusters with icosahedron-based geometries by removing or adding one atom from a complete icosahedron cluster, and how the structural changes can strongly cause internal energy to change accordingly. Because of the different structures of these clusters, the structure change processes of these types are quite different on heating.     

Local structure changes of Cu55 cluster during heating

The structural relaxation of a cluster containing 55 atoms at elevated temperatures is simulated by molecular dynamics. The interatomic interactions are given by using the embedded atom method （EAM） potential. By decomposing the peaks of the radial distribution functions （RDFs） according to the pair analysis technique, the local structural patterns are identified for this cluster. During increasing temperature, structural changes of different shells determined by atom density profiles result in an abrupt increase in internal energy. The simulations show how local structural changes can strongly cause internal energy to change accordingly.     

Evolution of three-shell onion-like and core-shell structures in （AgCo）201 bimetallic clusters

<正>This paper studies the structural evolution of（AgCo）₂₀₁ clusters with different Co concentrations under various temperature conditions by using molecular dynamics with the embedded atom method.The most stable position for Co atoms in the cluster is the subsurface layer at low temperature(lower than 200 K for the Ag₂₀₀Co₁ cluster).The position changes to the core layer with the increase of temperature,but there is an energy barrier in the middle layer. This makes the Ag-Co cluster form an Ag-Co-Ag three-shell onion-like configuration.When the temperature is high enough[higher than 800 K for（AgCo）₂₀₁ clusters with 50%Co],Co atoms can obtain enough energy to overcome the energy barrier and the cluster forms an Ag-Co core-shell configuration.Amorphization for the onion-like and core-shell clusters is induced by the large lattice misfit at Ag-Co interfaces.The structural evolution in the Ag-Co cluster is related to the release of excess energy.     

Fe基体中包含Cu团簇的Fe-Cu二元体系在升温过程中结构变化的原子尺度计算

采用基于嵌入原子方法的分子动力学方法模拟了具有体心立方晶格结构的Fe基体中包含小尺寸Cu纳米粒子的Fe-Cu二元体系在升温过程中的原子堆积结构变化.进行了Cu原子均方位移、Cu原子对分布函数和原子的径向密度分布函数的计算,并对纯Cu原子区、Fe-Cu界面区和纯Fe基体区的分区域原子堆积结构进行了分析.结果表明,Fe基体内Cu团簇的尺寸及其在Fe基体内所能占据区域的大小,对不同温度下的Cu团簇内原子堆积结构及Fe基体的原子堆积结构具有影响.升温过程中不同尺寸受基体约束Cu团簇对Fe基体结构改变的影响表现出很大差异.对于Fe_(bulk)-Cu_(135)体系,基体的应变临近Fe-Cu界面区,同时在团簇中间的基体区域出现大量空位缺陷和应变集中区;对于Fe_(bulk)-Cu_(141)体系,随温度升高,基体中出现的应变区域表现为小尺寸、数量多向大尺寸、小数量的变化.     

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     

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.     

熔融Cu55团簇在铜块体中凝固过程的分子动力学模拟

应用基于嵌入原子势函数的分子动力学方法,模拟了嵌入在具有面心立方结构同质块体中的熔融Cu₅₅团簇在不同急冷温度下微观结构的演变情况.通过计算熔融Cu₅₅团簇的均方位移和原子平均能量随时间步的变化,并应用键对分析技术,分析了急冷温度对熔融Cu₅₅团簇结构变化的影响.研究结果表明,由于受到块体结构的影响,在所研究的急冷温度范围内,熔融Cu₅₅团簇在凝固过程中形成了以面心立方结构为主的微观结构.结晶过程是原子不断交换其位置的过程,团簇原子位置的重排敏感于温度的变化.随着急冷温度的升高,原子的扩散范围增大.在100,300和500 K三个较低的温度下有利于形成稳定的面心立方结构,但当急冷到100 K时,团簇中的原子在没有找到其最佳位置之前就已经完成晶化.在急冷到500 K时,团簇中的原子在块体中扩散充分,与块体中的原子形成理想的面心立方结构.在700,900和1100 K三个较高的温度上,局域结构表现为随时间步波动性变化. 关键词： 团簇 分子动力学 计算机模拟 凝固     

A molecular dynamics study of the structural change differences between Au225 and Au369 clusters on MgO surfaces at low temperature

The differences in structural change between Au₂₂₅ and Au₃₆₉ clusters with their (111) facets supported on MgO(100) surfaces at 5 K are studied by using molecular-dynamics simulations with the atomic interchange potentials of the Au/MgO interface. The parameters are obtained from the ab initio energies using the Chen-Möbius inversion method. Analyses of the pair distribution functions show that the two Au clusters use different deformation processes to adjust the distances between the interface atoms, owing to the misfit between the atom distances among the clusters and the substrates. The local structural changes are identified by atomic density profiles.     

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

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

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

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

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

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

Processes involved in the formation of silver clusters on silicon surface

We analyze scanning electron microscopy measurements for structures formed in the deposition of solid silver clusters onto a silicon(100) substrate and consider theoretical models of cluster evolution onto a surface as a result of diffusion and formation of aggregates of merged clusters. Scanning electron microscopy (SEM) data are presented in addition to energy dispersive X-ray spectrometry (EDX) measurements of the these films. Solid silver clusters are produced by a DC magnetron sputtering source with a quadrupole filter for selection of cluster sizes (4.1 and 5.6 nm or 1900 and 5000 atoms per cluster in this experiment); the energy of cluster deposition is 0.7 eV/atom. Rapid thermal annealing of the grown films allows analysis of their behavior at high temperatures. The results exhibit formation of cluster aggregates via the diffusion of deposited solid clusters along the surface; an aggregate consists of up to hundreds of individual clusters. This process is essentially described by the diffusion-limited aggregation (DLA) model, and thus a grown porous film consists of cluster aggregates joined by bridges. Subsequent annealing of this film leads to its melting at temperatures lower than to the melting point of bulk silver. Analysis of evaporation of this film at higher temperatures gives a binding energy in bulk silver of ɛ₀= (2.74 ± 0.03) eV/atom. The text was submitted by the authors in English.     

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

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

Ag309团簇升温过程中结构转变的分子动力学模拟研究

采用分子动力学模拟研究了具有面心立方（fcc）晶格结构的截断八面体Ag309团簇升温过程中结构演变。对团簇的能量曲线变化、快照图演变和键对分析表明：无缺陷截断八面体Ag309团簇在410 K时转变为二十面体,在840 K时熔化;不同缺陷诱导二十面体结构转变温度异常变化,沿晶面滑移缺陷使二十面体转变温度升高,沿晶面旋转缺陷使二十面体结构转变温度降低;不同缺陷对团簇键型和势能产生的影响是使二十面体结构转变温度异常变化的主要诱导因素。这种通过缺陷控制团簇结构转变的研究为新型纳米结构的可控制备提供理论基础。     

Molecular dynamics calculation on three lithium atoms in planar cluster model C150H30 for glassy carbon

As the preliminary step to elucidate the charge-discharge mechanism in the glassy carbon cathode of lithium (Li) secondary battery, the molecular dynamics (MD) calculation was applied at molecular mechanics 2 (MM2) level for three Li atoms stabilized in the form of regular triangle in the central area of the hydrogen terminated planar cluster, C₁₅₀H₃₀. The stabilized Li sites are given on the basis of the structural optimization. Up to 10 K, the triangular Li aggregate kept at almost the same figure as that formed by three stabilization sites, rotating parallel to the cluster plane, goes around the central area of it, whereupon interatomic vibrational stretching is observed. Below 75 K, the aggregate of three Li atoms separates to a pair of two atoms and one Li atom arbitrarily, however, reformation of the pair occurs periodically among three atoms with the lapse of simulation time. Then, three Li atoms move correlatively irrespective of the long interatomic distance of 18 Å at maximum. However, at 100 K one Li atom goes out of the cluster model directly and the rest of two atoms continues the revolutional movement with the rotation as a pair in the central area of the cluster model. Thus, three Li atoms show the appreciably stable movements in the glassy carbon by forming the aggregate or the atomic pair, which will be responsible for the hysteresis in the charge-discharge cycle of lithium secondary battery.     

Effect of melting on ionization potential of sodium clusters

The effect of melting transition on the ionization potential has been studied for sodium clusters with 40, 55, 142, and 147 atoms, using ab initio and classical molecular dynamics. Classical and ab initio simulations were performed to determine the ionization potential of Na₁₄₂ and Na₁₄₇ for solid, partly melted, and liquid structures. The results reveal no correlation between the vertical ionization potential and the degree of surface disorder, melting, or the total energy of the cluster obtained with the ab initio method. However, in the case of 40 and 55 atom clusters, the ionization potential seems to decrease when the cluster melts. Received 1st November 2002 Published online 24 April 2003 RID="a" ID="a"e-mail: ar@phys.jyu.fi     

小尺寸铜团簇冷却与并合过程中结构变化的原子尺度研究

研究Cu_N（N=57,58,59）熔融铜团簇在冷却过程以及300 K时两个具有二十面体结构Cu₅₅团簇在并合过程中的结构变化.对这些小尺寸团簇的结构变化采用基于嵌入原子方法的正则系综分子动力学进行计算机模拟.通过对模拟结果的分析表明,小团簇的冷却和并合过程存在阶段变化的特点.降温过程中Cu_N（N=57,58,59）团簇的原子运动及其微观结构变化表现出较大差异,由此导致这三类团簇内原子排布的不同,其中Cu₅₉团簇结构的有序程度最低.在两个Cu₅₅团簇并合早期阶段,这两个团簇相接触后发生变形导致原子位置出现较大改变,在随后的并合过程中,原子扩散引起原子局部位置调整导致所并合体系的结构发生变化.远离两个团簇接触区的原子仍保持其并合前的结构. 关键词： 团簇 分子动力学 计算机模拟 表面     

Structure,energetic and phase transition of small nickel-palladium heterogeneous clusters

Molecular dynamics simulation (MD) with Sutton-Chen potential for palladium-palladium, nickel-nickel and palladium-nickel interactions has been used to generate the minimum energy structures and to study the thermodynamic and dynamic properties of mixed transition metal cluster motifs of Ni _n Pd_(13?n) for n ≤ 13. Thirteen particle icosahedral clusters of neat palladium and nickel atoms were first reproduced accordingly with the results in literature. Then in the palladium icosahedra, each palladium atom has been successively replaced by nickel atom. Calculation is repeated for both palladium-centered and nickel-centered clusters. It is found that the nickel-centered clusters are more stable than the palladium-centered clusters and cohesive energy increases along the palladium end to nickel end. Phase transition of each cluster from one end-species to the other end-species is studied by means of caloric curve, root mean square bond fluctuation and heat capacity. Trend in variation of melting temperature is opposite to the energy trend. Palladium-centered cluster shows a premelting at low temperature due to the solid-solid structural transition. Species-centric order parameters developed by Hewage and Amar is used to understand the dynamic behavior in the solid-solid transition of palladium-centered cluster to more stable nickel-centered cluster (premelting). This species-centric order parameter calculation further confirmed the stability of nickel-centered species over those of palladium-centered species and solid-solid structural transition at low temperature.