液态铝基合金微观结构及其性质的分子动力学模拟研究

利用分子动力学研究了液态铝基合金Al80Fe20和Al80Cu20在 800~2200K温度范围内的微观结构和性质的变化。扩散系数的计算结果显示两合金有着不同的动力学行为。温度大约1400K以下液态Al80Cu20的扩散系数随温度的变化比Arrhenius关系所预期的要慢些,所得结果与Brillo研究小组的实验结果基本一致。当液态铝基合金的动力学发生异常变化时,其微观二十面体短程序的含量和作用也发生相应的变化。因此,熔体的奇特动力学行为可能由于其微观局域结构短程序随着温度的变化而异常地改变而引起的。     

液态铝基合金微观结构及其性质的分子动力学模拟研究(英文)

利用分子动力学研究了液态铝基合金Al80Fe20和Al80Cu20在800~2200K温度范围内的微观结构和性质的变化.扩散系数的计算结果显示两合金有着不同的动力学行为.温度大约1400K以下液态Al80Cu20的扩散系数随温度的变化比Arrhenius关系所预期的要慢些,所得结果与Brillo研究小组的实验结果基本一致.当液态铝基合金的动力学发生异常变化时,其微观二十面体短程序的含量和作用也发生相应的变化.因此,熔体的奇特动力学行为可能由于其微观局域结构短程序随着温度的变化而异常地改变而引起的.     

One-dimensional atomic transport by clusters of self-interstitial atoms in iron and copper

Atomic-scale computer simulation has been used to study the thermally activated atomic transport of self-interstitial atoms (SIAs) in the form of planar clusters in pure Cu and f-Fe. There is strong evidence that such clusters are commonly formed in metals during irradiation with high-energy particles and play an important role in accumulation and spatial distribution of surviving defects. An extensive study of the mobility of SIA clusters containing two to 331 interstitials has been carried out using the molecular dynamics simulation technique for the temperature range from 180 to 1200 K. The results obtained show that clusters larger than three to four SIAs are one-dimensionally mobile in both Cu and Fe. Large clusters of more than 100 SIAs in Cu and 300 SIAs in Fe have significantly reduced mobility. The problem of describing one-dimensional (1D) motion in three-dimensional space is discussed. An attempt is made to describe the mobility of SIA clusters within the approximation of 1D diffusion. For clusters in both metals the effective migration energy of 1D diffusion as estimated via the jump frequency of the cluster centre of mass is found to be independent of the number of SIAs in the clusters, although the cluster jump frequency decreases with increasing cluster size. Mechanisms of 1D mobility of interstitial clusters are discussed.     

Cluster size and substrate temperature affecting thin film formation during copper cluster deposition on a Si （001） surface

The soft deposition of Cu clusters on a Si(001) surface was studied by molecular dynamics simulations.The embedded atom method,the Stillinger-Weber and the Lennar-Jones potentials were used to describe the interactions between the cluster atoms,between the substrate atoms,and between the cluster and the substrate atoms,respectively.The Cu13,Cu55,and Cu147 clusters were investigated at different substrate temperatures.We found that the substrate temperature had a significant effect on the Cu147 cluster.For smaller Cu13 and Cu55 clusters,the substrate temperature in the range of study appeared to have little effect on the mean center-of-mass height.The clusters showed better degrees of epitaxy at 800 K.With the same substrate temperature,the Cu55 cluster demonstrated the highest degree of epitaxy,followed by Cu147 and then Cu13 clusters.In addition,the Cu55 cluster showed the lowest mean center-of-mass height.These results suggested that the Cu55 cluster is a better choice for the thin-film formation among the clusters considered.Our studies may provide insight into the formation of desired Cu thin films on a Si substrate.     

Effect of Ni and vacancy concentration on initial formation of Cu precipitate in Fe–Cu–Ni ternary alloys by molecular dynamics simulation

In the present work, the effects of Ni atoms and vacancy concentrations(0.1%, 0.5%, 1.0%) on the formation process of Cu solute clusters are investigated for Fe–1.24%Cu–0.62%Ni alloys by molecular dynamics(MD) simulations. The presence of Ni is beneficial to the nucleation of Cu precipitates and has little effect on coarsening rate in the later stage of aging. This result is caused by reducing the diffusion coefficient of Cu clusters and the dynamic migration of Ni atoms. Additionally, there are little effects of Ni on Cu precipitates as the vacancy concentration reaches up to 1.0%,thereby explaining the embrittlement for reactor pressure vessel(RPV) steel. As a result, the findings can hopefully provide the important information about the essential mechanism of Cu cluster formation and a better understanding of ageing phenomenon of RPV steel. Furthermore, these original results are analyzed with a simple model of Cu diffusion, which suggests that the same behavior could be observed in Cu-containing alloys.     

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

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

Mass and thermal transport in liquid Cu-Ag alloys

In this paper, the diffusion, thermodynamic and thermotransport properties in Cu–Ag liquid alloys are extensively investigated with molecular dynamics over a wide composition and temperature range. The simulations are performed with the most reliable EAM potential. The Green-Kubo formalism is employed for calculating transport properties. It is found that the reduced heat of transport in Cu–Ag is very small (about 0.10?eV in absolute value) and almost temperature independent. Further it is found that the interdiffusion coefficient together with both self-diffusion coefficients are almost composition independent. In Cu–Ag, the thermodynamic factor is found to be less than unity whereas the Manning factor is greater than unity (with significant composition and temperature dependence) and their product is very close to 1.     

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

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

The molecular dynamic study of anharmonic effects at Cu(111) and Ag(111) surfaces in the presence of Cu- and Ag-trimer island

The molecular dynamics (MD) technique based on semi-empirical potentials, is used to carry out the diffusion of Cu- and Ag-trimer on Cu- and Ag(111) surface at 300, 500 and 700 K temperatures. The constant energy MD simulation elaborates the anharmonic effects at the surface such as fissures, dislocations and vacancy creation, in the presence of island. The fissures and dislocations formed are in the range of 1.5–4 Å and 1–7 Å, respectively, from the island's position. The Cu and Ag islands both diffuse easily on Cu(111) surface, manipulate that the trend of diffusion is faster on Cu surface as compared to Ag surface. The process of breaking and opening of the island has also been observed. Moreover, a surface atom popped-up at 700 K by creating a vacancy near the Cu island on Ag surface. The rate of diffusion increases with the increase in temperature, both for homo- and hetero-cases.     

Surface diffusion of Pb clusters on Cu(1 1 1) influenced by size dependent cluster-substrate misfit: Monte Carlo tight binding simulation model

The migration of two-dimensional Pb clusters on Cu(1 1 1) surface is studied in the framework of three-dimensional (3D) continuum space tight binding (TB) computational model. Monte Carlo simulations based on many-body TB interactions reveal significant influence of cluster-substrate misfit on diffusion behavior of two-dimensional (2D) clusters. The analysis of pair distribution function (PDF) demonstrates that cluster lattice constant depends on the number of atoms N for 1 < N < 10. The observed effect is more pronounced for heteroepitaxial than homoepitaxial systems. It can be explained in the framework of model accounting for enhanced charge transfer in heteroepitaxy and strong influence of surface potential on cluster atomic arrangement. The variation of lattice constant leads to variable misfit which affects the island migration behavior. The results are discussed in a physical model that implies cluster diffusion with size dependent cluster-substrate misfit.     

Ni diffusion in near-equiatomic Ni-Ti and Ni-Ti(-Cu) alloys

Tracer volume diffusion of 63 Ni in Ni-50.07 at.% Ti binary and Ni-48.83 at.% Ti-9.29 at.% Cu ternary alloys have been measured between 783 and 1288 K. The temperature dependence of the diffusion coefficients can be well described by a straight Arrhenius function This illustrates that a diffusional anomaly (typical in some bcc metallic alloys and related to the well-known phonon softening), if there is any in this system, is negligible. On the other hand, the small activation energy (about half the value expected from simulations for the commonly accepted mechanism with thermally activated vacancies) shows that the mechanism of diffusion is probably mediated by structural vacancies.     

Dynamical heterogeneity and diffusion in high-density Al2O3·2SiO2 melts

Dynamical heterogeneity (DH) in high-density Al₂O₃·2SiO₂ melts has been studied in a model containing 3025 atoms via molecular dynamics (MD) simulation and at the fixed density of 4.0 g/cm³. Non-Gaussian parameter of atomic species in the system has been found and discussed. We found a clear evidence of the existence of DH in high-density Al₂O₃·2SiO₂, which has specific features differed from those observed in the lower-density one. The most mobile and immobile atoms in the system have a tendency to form clusters and temperature dependence of their mean cluster size was found. On the other hand, diffusion constant of atomic species in the system has been calculated at temperatures ranged from 3150 to 7000 K. Calculations show that at relatively not high temperatures, temperature dependence of diffusion constant shows an Arrhenius law and at higher temperatures it shows a power law: D∝(T−T_C)^γ. Diffusion data of high-density melts have been compared with those for the low-density ones. Diffusion mechanism in the system has been discussed via the temperature dependence of diffusion constant ratio and activation energy. And we found the existence of cooperative diffusion mechanism in the system.     

低能单个Cu原子与Cu13团簇沉积到Fe（001）表面的比较

本文利用分子动力学模拟方法对相同初始沉积条件下的单个Cu原子和Cu13团簇与Fe(001)表面的相互作用分别进行了模拟研究, 并将两者的模拟结果进行了比较分析. 单个Cu原子和Cu13团簇的初始入射能量范围均为1eV/atom、3eV/atom、5eV/atom和10eV/atom, 初始入射角度均为0o、10o、30o和45o, 衬底温度分别为100K、300K和800K. 对单个Cu原子和Cu13团簇的原子动能、质心高度、迁移距离和最终沉积形貌进行了分析, 对比研究了相同初始沉积条件下单个Cu原子和Cu13团簇在沉积过程中和沉积效果上的具体差异. 模拟结果表明: 单个Cu原子和Cu13团簇与Fe(001)表面的相互作用机制存在差异, Cu13团簇表现出显著的集体效应. 在特定沉积条件下, 由于Cu13团簇的集体效应, 导致Cu13团簇与Fe(001)表面的结合能力和在Fe(001)表面上的扩散能力均强于单个Cu原子.     

An investigation of the internal temperature dependence of Pd-Pt cluster beam deposition: A molecular dynamics study

We investigated the internal temperature dependence of the Pd_1−aPt_a cluster beam deposition in the present study via the molecular dynamics simulations of soft-landing. By analysis of the velocity distribution and diffusion coefficient of the bimetallic cluster, Pd atoms with better mobility improved the diffusibility of Pt atoms. The radial composition distribution showed that a Pt-core/Pd-shell structure of the cluster formed at high internal temperatures through migrations of the Pd atoms from inner to surface shells. In the soft-landing process, the diffusing and epitaxial behaviors of the deposited clusters mainly depended on the internal temperature because the incident energy of the cluster was very small. By depositing clusters at high internal temperatures, we obtained a thin film of good epitaxial growth as the energetic cluster impact. Furthermore, nonepitaxial configurations such as scattered nonepitaxial atoms, misoriented particles, and grain boundaries of (1 1 1) planes were produced in the growth of the cluster-assembled film. As the size of the incident cluster increased, the internal temperature of the cluster needed for better interfacial diffusion and contact epitaxy on the substrate also rose.     

Magnetic dynamics of dilute iron nano-clusters in silver films from Mössbauer spectroscopy and muon spin rotation

A silver film containing nanometer size clusters of iron (nominal conc. 1 at%) has been studied by Mössbauer spectroscopy and Low-Energy Muon Spin Rotation. Below about 20 K spin glass freezing due to interparticle interactions is found from both methods. Whereas Mössbauer spectra are insensitive to the fast fluctuations of cluster moments above spin glass freezing temperature, muon spin rotation in magnetic fields applied perpendicular to the polarized muon spins allows tracing the fluctuations of superparamagnetic moments. The temperature dependence of the damping of the muon spin rotation signal shows Arrhenius behavior between 10 to 100 K. Depending on the assumed shape of damping the activation energy of superparamagnetic fluctuations of cluster moments ranges between about 20 K ·k _B and 40 K ·k _B . Above about 120 K muon spin depolarization indicates diffusion and trapping of muons.     

Misfit-dislocation-mediated heteroepitaxial island diffusion

Scanning tunneling microscopy combined with molecular dynamics simulations reveals a dislocation-mediated island diffusion mechanism for Cu on Ag(111), a highly mismatched system. Cluster motion is tracked with atomic precision at multiple temperatures and diffusion barriers and prefactors are determined from direct measurements of hop rates. The barrier to nucleate a dislocation is sensitive to island size and shape, resulting in a non-monotonic size dependence of the diffusion barrier.     

Simulation of nanostructure evolution under helium implantation in Fe–(2.5–12.5)at% Cr alloys at a temperature of 343 K

Characterization of helium-implanted Fe–(2.5–12.5)at% Cr alloys with the flux of 7?×?10^–6?dpa/s at a temperature of 343?K has been performed by means of cluster dynamics simulations. We have suggested a model for simulating an Fe–C–Cr system under helium implantation based on a selection of the latest data from atomistic studies and available experiments. Kinetics of carbon-vacancy and helium-vacancy complexes has been studied. Only one parameter is used to guarantee the best reproduction possible of experimental positron annihilation lifetime spectroscopy data for Fe–Cr alloys on dependences of vacancy cluster size on chromium content and irradiation dose via fitting. This is an effective binding energy of self-interstitial atoms to dislocation loops decorated by chromium atoms. It has a “snaky” dependence of chromium content with a minimum of about 9%Cr.     

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

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

Molecular dynamics simulation of a graphite-supported copper nanocluster: thermodynamic properties and gas adsorption

Molecular dynamics simulations were conducted for a cubic Cu cluster supported on a graphite bilayer. The Sutten–Chen and Lennard–Jones potentials were used for metal–metal and metal–graphite interactions, respectively. Heating and cooling processes were performed by NVT simulations at different temperatures in the range 200 to 1800?K. The melting point was identified on the basis of caloric and heat capacity curves. The calculated melting point was 770?K, far below the bulk melting point of crystalline copper. Several phenomena such as the appearance of a hysteresis (irreversibility) in caloric curves, surface melting, and cluster-induced surface wetting were justified from the results. The simulation of cluster in the presence of gas atmosphere showed that the CO gas is adsorbed more than H₂ and it has a greater impact on the cluster's structure.