Calculation of the surface energy of fcc metals with modified embedded-atom method

The surface energies for 38 surfaces of fcc metals Cu, Ag, Au, Ni, Pd, Pt, A1, Pb, Rh and Ir have been calculated by using the modified embedded-atom method. The results show that, for Cu, Ag, Ni, A1, Pb and Ir, the average values of the surface energies are very close to the polycrystalline experimental data. For all fcc metals, as predicted, the close-packed (111) surface has the lowest surface energy. The surface energies for the other surfaces increase linearly with increasing angle between the surfaces (hkl) and (111). This can be used to estimate the relative values of the surface energy.     

fcc金属表面能的各向异性分析及表面偏析的预测

本文将元素变量(φ^*和n_WS)和MAEAM相结合,从原子尺度上对10种fcc金属Cu,Ag,Au,Ni,Pd,Pt,Rh,Al,Ir和Pb的38个不同晶面的表面能进行模拟计算及各向异性分析. 结果表明,fcc金属的密排面(111)的表面能最小,则该晶粒取向优先生长,与实验结果和第一原理的LMTO-ASA计算结果一致;各个晶面的表面能均随着其他晶面与(111)晶面的夹角cosθ_(hkl)的增长而呈线性 关键词： FCC金属 MAEAM 表面能 表面偏析     

Molecular dynamics simulation of surface melting behaviours of the V（110） plane

The modified analytic embedded-atom method and molecular dynamics simulations are applied to the investigation of the surface premelting and melting behaviours of the V（110） plane by calculating the interlayer relaxation, the layer structure factor and atomic snapshots in this paper. The results obtained indicate that the premelting phenomenon occurs on the V（110） surface at about 1800K and then a liquid-like layer, which approximately keeps the same thickness up to 2020K, emerges on it. We discover that the temperature 2020K the V（110） surface starts to melt and is in a completely disordered state at the temperature of 2140K under the melting point for the bulk vanadium.     

用改进嵌入原子法计算Cu晶体的表面能

用改进嵌入原子法(MEAM)计算了Cu晶体12个晶面的表面能.结果表明，密排面(111)的表面能最小.其他晶面的表面能随其晶面与(111)晶面夹角的增加而增加，据此可以粗略地估计各晶面表面能的相对大小.给出的几何结构因子的确定方法及结果可以直接用于计算其他面心立方晶体的表面能及其他特性.在Cu,Ag等面心立方薄膜中出现(111)择优取向或织构的机理是表面能的最小化. 关键词： 改进嵌入原子法 铜 表面能 计算     

Atomic simulation of the vacancies in BCC metals with MAEAM

The formation energy of the mono-vacancy and both the formation energy and binding energy of the di-and tri-vacancy in BCC alkali metals and transition metals have been calculated by using the modified analytical embedded-atom method (MAEAM). The formation energy of each type of configuration of the vacancies in the alkali metals is much lower than that in the transition metals. From minimum of the formation energy or maximum of the binding energy, the favorable configuration of the di-vacancy and tri-vacancy respectively is the first-nearest-neighbor (FN) or second-nearest-neighbor (SN) di-vacancy and the [112] tri-vacancy constructed by two first-and one second-nearest-neighbor vacancies. It is indicated that there is a concentration tendency for vacancies in BCC metals.     

Calculation of the formation energies of isolated vacancy and adatom-vacancy pair at low-index surfaces of fcc metals with MAEAM

The formation energies of isolated vacancy and adatom-vacancy pair (where the two are separated by a large distance) at low-index surfaces of fcc metals calculated by using the modifies analytical embedded atom method (MAEAM). The results predict the prevailing formation of vacancies on the surfaces (1 1 1), (1 0 0) (but Pd), Cu and Ni on the (1 1 0) surfaces at low temperature, and the defect formation energies consistently create in the sequence (1 1 0) → (1 0 0) → (1 1 1). With good accuracy, the calculated energy values coincide with those obtained by the embedded atom methods (EAM) and from experiments. The correctness of the method by which calculated the formation energies of point defects on the surface was proved.     

Mo纳米薄膜热力学性质的分子动力学模拟

采用改进嵌入原子法(MAEAM),通过经典的分子动力学(MD)模拟计算了高熔点过渡金属体心立方(bcc) Mo块体Gibbs自由能和表面能. 对于纳米薄膜的热力学数据,比如Gibbs自由能等,可以看成是薄膜内部原子和表面原子两部分数据之和,然后根据薄膜的体表原子之比就可以直接计算出总的自由能,并由此可以得到热力学性质与薄膜尺寸及温度的定量关系式. 分别计算了bcc Mo块体及其纳米尺寸薄膜的自由能和热容,结果表明,Mo纳米薄膜的热力学性质具有尺寸效应,并且在薄膜尺寸小于15—20nm时,这种效应变得非常明 关键词： 改进嵌入原子法 Mo纳米薄膜 表面自由能 热容     

In在Al(001)表面偏析的模拟

通过一种空位模型详细的描述了In在Al(001)表面的扩散偏析过程,利用周期性密度泛函理论方法计算了这个偏析过程中每步构型的能量和In原子扩散的能量壁垒,并对可能的偏析机理进行分析.结果表明:In原子从Al(001)表面第二层扩散偏析至表面层时,系统的能量降低了0.64 eV,最大的扩散迁移壁垒为0.34 eV;而从表面更内层向表面第二层扩散时系统能量基本保持不变,扩散需要克服的能量壁垒为0.65 eV,说明In原子在Al(001)表面只能由体内向表面扩散偏析.In在Al(001)的清洁表面具有强烈的偏析趋势,在热力学上是容易进行的. 关键词： 密度泛函理论 表面偏析 扩散 Al合金     

金属fcc(100)表面吸附原子的自扩散现象

通过静态计算和分子动力学方法 ,系统地研究了吸附原子在金属fcc( 10 0 )表面的自扩散现象 ,所采用的模型势是一系列由原子嵌入方法 (EAM )得到的。结果表明金属fcc( 10 0 )表面可分为三种不同的类型 ,即硬、中性和软。而且发现在大多数表面 ,如硬的和许多中性表面的简单交换机制都是通过“应力产生和释放”的模型 ,而不是以前所认为“合作运动”模型 ,只有在软表面“合作运动”模型才变得重要。另外 ,在软表面上还观察到了另一种形式的复杂交换机制 ,它不同于“应力产生和释放”模型基础上的复杂交换机制。     

Surface segregation of Al substrate metal on Zr film surface

Segregation of substrate Al on thin Zr film, Zr/Al/Al system was investigated by heating the specimen in a UHV chamber. Dual-cathode magnetron-sputtering source was used for deposition of Zr film as well as thin Al film to avoid aluminum oxide formation at Zr/Al interface. Al segregates on Zr film surface at 730 K. It was found that oxide-free interface between film and substrate is important for segregation in Al system. The diffusion coefficient calculated for surface segregation and inter-metallic compound showed that the grain boundary diffusion and bulk diffusion are very close in Zr/Al/Al system. Hence, it is important to control specimen heating to cause surface segregation by grain boundary diffusion.     

Phonon dispersion on Ag(100) surface:A modified analytic embedded atom method study

Within the harmonic approximation, the analytic expression of the dynamical matrix is derived based on the modified analytic embedded atom method(MAEAM) and the dynamics theory of surface lattice. The surface phonon dispersions along three major symmetry directionsΓˉXˉ, Γˉ Mˉ, and Xˉ Mˉ are calculated for the clean Ag(100) surface by using our derived formulas. We then discuss the polarization and localization of surface modes at points Xˉ and Mˉ by plotting the squared polarization vectors as a function of the layer index. The phonon frequencies of the surface modes calculated by MAEAM are compared with the available experimental and other theoretical data. It is found that the present results are generally in agreement with the referenced experimental or theoretical results, with a maximum deviation of 10.4%. The agreement shows that the modified analytic embedded atom method is a reasonable many-body potential model to quickly describe the surface lattice vibration. It also lays a significant foundation for studying the surface lattice vibration in other metals.     

用分析型嵌入原子方法计算Pd-Au合金的热力学性质

用修正的分析型嵌入原子方法计算了PdAu二元合金系的稀溶解热、形成焓、表面能和表面聚集能,并将计算的热力学性质与实验结果或第一性原理计算结果进行了比较.     

Surface segregation energies in low-index open surfaces of bimetallic transition metal alloys

We present a database of 24 × 24 segregation energies of single transition metal impurities in low-index surfaces of transition metal hosts, calculated using the localized self-consistent Green’s function (LSGF) method, in combination with the atomic sphere approximation including a multipole correction to the electrostatic potential and energy. The surface energies of facets for fcc and bcc transition metals, and the more stable of the two facets of hcp transition metals are also calculated and compared with available theoretical results. Insights derived should be useful for determining the nature of active sites in a variety of catalytic reactions employing bimetallic catalysts.     

Influence of surface stress in the missing row reconstruction of fcc transition metals

Surface stress relief is often invoked as the driving force for surface reconstruction. The present ab initio total energy calculations show that it is not the case, at least for the missing row reconstruction undergone by the (1 1 0) surface of the face centered cubic transition metals. Instead, they show that the tendency to reconstruct or not is closely related to the sign and magnitude of the strain derivative of surface energy. Finally, they unambiguously confirm that relativistic effects are at the origin of this reconstruction.     

Effects of rapid thermal annealing on crystallinity and Sn surface segregation of Ge_(1-x)Sn_x films on Si(100) and Si(111)

Germanium-tin films with rather high Sn content(28.04% and 29.61%) are deposited directly on Si(100) and Si(111)substrates by magnetron sputtering. The mechanism of the effect of rapid thermal annealing on the Sn surface segregation of Ge_(1-x)Sn_x films is investigated by x-ray photoelectron spectroscopy(XPS) and atomic force microscopy(AFM). The x-ray diffraction(XRD) is also performed to determine the crystallinities of the Ge_(1-x)Sn_x films. The experimental results indicate that root mean square(RMS) values of the annealed samples are comparatively small and have no noticeable changes for the as-grown sample when annealing temperature is below 400℃. The diameter of the Sn three-dimensional(3 D) island becomes larger than that of an as-grown sample when the annealing temperature is 700℃. In addition, the Sn surface composition decreases when annealing temperature ranges from 400℃ to 700℃. However, Sn bulk compositions in samples A and B are kept almost unchanged when the annealing temperature is below 600℃. The present investigation demonstrates that the crystallinity of Ge_(1-x)Sn_x/Si(111) has no obvious advantage over that of Ge_(1-x)Sn_x/Si(100) and the selection of Si(111) substrate is an effective method to improve the surface morphologies of Ge_(1-x)Sn_x films. We also find that more severe Sn surface segregation occurs in the Ge_(1-x)Sn_x/Si(111) sample during annealing than in the Ge_(1-x)Sn_x/Si(100) sample.     

A method for the study of surface segregation in multicomponent alloys

A simple algorithm for the determination of segregation profiles in multicomponent systems based on a mean field formalism and a quantum approximate method for the energetics is introduced. The method is described and applied to two ternary systems, concentrating on the changes in segregation patterns relative to the corresponding binary cases.     

Stress development and impurity segregation during oxidation of the Si(1 0 0) surface

We have studied the segregation of P and B impurities during oxidation of the Si(1 0 0) surface by means of combined static and dynamical first-principles simulations based on density functional theory. In the bare surface, dopants segregate to chemically stable surface sites or to locally compressed subsurface sites. Surface oxidation is accompanied by development of tensile surface stress up to 2.9 Nm⁻¹ at a coverage of 1.5 monolayers of oxygen and by formation of oxidised Si species with charges increasing approximately linearly with the number of neighbouring oxygen atoms. Substitutional P and B defects are energetically unstable within the native oxide layer, and are preferentially located at or beneath the Si/SiO_x interface. Consistently, first-principles molecular dynamics simulations of native oxide formation on doped surfaces reveal that dopants avoid the formation of P-O and B-O bonds, suggesting a surface oxidation mechanism whereby impurities remain trapped at the Si/SiO_x interface. This seems to preclude a direct influence of impurities on the surface electrostatics and, hence, on the interactions with an external environment.     

Surface properties of 3d transition metals

Using the projector augmented wave method within density functional theory, we present a systematic study of the layer relaxation, surface energy and surface stress of 3d transition metals. Comparing the calculated trends for the surface energy and stress with those obtained for 4d and 5d metals we find that magnetism has a significant effect on the surface properties. Enhanced surface magnetic moments decrease the size of the surface relaxation, lower the surface energy and surface stress, leading to compressive stress in Cr and Mn.     

Interfacial segregation in Ag-Au,Au-Pd,and Cu-Ni alloys: I. (100) surfaces

Atomistic simulations of segregation to (100) free surface in Ag–Au, Au–Pd, and Cu–Ni alloy systems have been performed for a wide range of temperatures and compositions within the solid solution region of these alloy phase diagrams. In addition to the surface segregation profiles, surface free energies, enthalpies, and entropies were determined. These simulations were performed within the framework of the free energy simulation method, in which an approximate free energy functional is minimized with respect to atomic coordinates and atomic site occupation. The effects of the relaxation with respect to either the atomic positions or the atomic concentrations are discussed. For all alloy bulk compositions (0.05 C 0.95) and temperatures (400 T(K) 1,100) examined, Ag, Au, and Cu segregates to the surface in the Ag–Au, Au–Pd, and Cu–Ni alloy systems, respectively. The present results are compared with several theories for segregation. The resultant segregation profiles in Au–Pd and Ag–Au alloys are shown to be in good agreement with an empirical segregation theory, while in Cu–Ni alloys the disagreement in Ni-rich alloys is substantial. The width of the segregation profile is limited to approximately three to four atomic planes. The surface thermodynamic properties depend sensitively on the magnitude of the surface segregation, and some of them are shown to vary linearly with the magnitude of the surface segregation.