Improvement of modified analytic embedded atom method potentials for noble metals and Cu

The modified analytic embedded atom model (EAM) potentials considering farther neighbor atoms are improved for the noble metals (Ag, Au, Pt, Pd, Rh) and Cu. We not only adopt an end processing function and an enhanced smooth continuous condition for the pair potential, but also adjust the model parameters of multi-body potential by fitting a cohesive energy, a mono-vacancy formation energy, the Rose equation curve for the cohesive energy as a function of lattice parameter, a structure energy difference, elastic parameters and an equilibrium condition of crystal. The calculation results of structure energy differences misfit the experiment data for the noble metals and Cu in the unimproved EAM, because anyone of these differences have not been considered in the calculation of its model parameters. After the modification, the model showed better simulation results for the noble metals and Cu.     

Vibrational properties of vacancy in bcc transition metals using embedded atom method potentials

The embedded atom method (EAM) potentials, with the universal form of the embedding function along with the Morse form of pair potential, have been employed to determine the potential parameters for three bcc transition metals: Fe, Mo, and W, by fitting to Cauchy pressure (C ₁₂???C ₄₄)/2, shear constants $G_\textrm{v} =({C_{11} -C_{12} +3C_{44}})/5$ and C ₄₄, cohesive energy and the vacancy formation energy. The obtained potential parameters are used to calculate the phonon dispersion spectra of these metals. Large discrepancies are found between the calculated results of phonon dispersion using the EAM and the experimental phonon dispersion results. Therefore, to overcome this inadequacy of the EAM model, we employ the modified embedded atom method (MEAM) in which a modified term along with the pair potential and embedding function is added in the total energy. The phonon dispersions calculated using potential parameters obtained from the MEAM show good agreement with experimental results compared to those obtained from the EAM. Using the calculated phonons, we evaluate the local density of states of the neighbours of vacancy using the Green’s function method. The local frequency spectrum of first neighbours of vacancy in Mo shows an increase at higher frequencies and a shift towards the lower frequencies whereas in Fe and W, the frequency spectrum shows a small decrease towards higher frequency and small shift towards lower frequency. For the second neighbours of vacancy in all the three metals, the local frequency spectrum is not much different from that of the host atom. The local density of states of the neighbours of the vacancy has been used to calculate the mean square displacements and the formation entropy of vacancy. The calculated mean square displacements of the first neighbours of vacancy are found to be higher than that of the host atom, whereas it is lower for the second neighbours. The calculated results of the formation entropy of the vacancy compared well with other available results.     

Analysis of semi-empirical interatomic potentials appropriate for simulation of crystalline and liquid Al and Cu

We investigate the application of embedded atom method (EAM) interatomic potentials in the study of crystallization kinetics from deeply undercooled melts, focusing on the fcc metals Al and Cu. For this application, it is important that the EAM potential accurately reproduces melting properties and liquid structure, in addition to the crystalline properties most commonly fit in its development. To test the accuracy of previously published EAM potentials and to guide the development of new potential in this work, first-principles calculations have been performed and new experimental measurements of the Al and Cu liquid structure factors have been undertaken by X-ray diffraction. We demonstrate that the previously published EAM potentials predict a liquid structure that is too strongly ordered relative to measured diffraction data. We develop new EAM potentials for Al and Cu to improve the agreement with the first-principles and measured liquid diffraction data. Furthermore, we calculate liquid-phase diffusivities and find that this quantity correlates well with the liquid structure. Finally, we perform molecular dynamics simulations of crystal nucleation from the melt during quenching at constant cooling rate. We find that EAM potentials, which predict the same zero-temperature crystal properties but different liquid structures, can lead to quite different crystallization kinetics. More interestingly, we find that two potentials predicting very similar equilibrium solid and liquid properties can still produce very different crystallization kinetics under far-from-equilibrium conditions characteristic of the rapid quenching simulations employed here.     

Liquid/vapor surface tension of metals: embedded atom method with charge gradient corrections

Molecular dynamics simulations for three embedded atom method (EAM) function sets are used to determine the liquid/vapor surface tension gamma for Al, Ni, Cu, Ag, and Au. The three EAM models differ in both the functional forms employed and the fitting procedure used. All the EAM potentials underestimate gamma but one of the models performs consistently better than the others. We show that including a correction to the local charge density associated with gradients in the density together with exploiting the invariance of the EAM potentials to appropriate transformations in the charge density can lead to improved values for gamma, as well as for solid free surface energies, within existing EAM function sets.     

金和银的晶格反演势的构建及应用

如何确定精确的原子间作用势一直是模拟计算的重要基础问题. 以面心立方金属金和银为对象, 采用第一性原理方法, 分别得到金和银的晶格内聚能-原子距离曲线及基态原子能曲线. 根据陈-莫比乌斯(Chen-Möbius)晶格反演理论和自编程序, 得到了精确的反演对势曲线. 对该曲线进行拟合, 通过比较不同势函数的拟合结果, 提出了双指数型的势函数解析式, 并得到具有全局性且高精度的拟合效果. 为了验证反演势的有效性, 利用反演势结果计算了金和银的声子谱并与Sutton-Chen提出的嵌入原子势和第一性原理计算得到的声子谱做比较. 分析表明, 反演势能够合理反映原子间的相互作用. 最后, 利用得到的结果, 计算了金和银的热膨胀系数和弹性模量等物理量, 计算结果与实验数据基本符合, 表明构建的金和银的反演势是准确有效的.     

Dynamics of Liquid Lithium Atoms. Pseudopotential and EAM-Type Potentials

It is generally accepted that the complicated character of the interparticle interaction in liquid metals is reproduced most correctly by many-particle potentials of the EAM-type (embedded atom model) interparticle interaction. It is shown that in the case of liquid lithium near the melting temperature (T_m = 453.65 K), the spherical pseudopotential provides a better agreement with experimental data on elastic and inelastic X-ray scattering as compared to the known EAM potentials. The calculations of the dynamic structural factor and spectral densities of the longitudinal and transverse atomic currents lead to the conclusion that although the pseudopotential and EAM potentials generate a certain qualitative correspondence in the features of collective dynamics, the interparticle interaction of the spherical type reproduces correctly the general form of the dynamic structure factor in a certain wavenumber range, as well as the dispersion relation for collective excitations.     

Thermodynamic and physical properties of FeAl and Fe3Al: an atomistic study by EAM simulation

With this work we present a newly developed potential for the Fe–Al system, which is based on the analytical embedded atom method (EAM) with long range atomic interactions. The potential yields for the two most relevant phases B2-FeAl and D0₃-Fe₃Al lattice constants, elastic constants, as well as bulk and point defect formation enthalpies, which are in good agreement with experimental and other theoretical data. In addition, the phonon dispersions for B2-FeAl and D0₃-Fe₃Al show a good agreement with available experiments. The calculated lattice constants and formation enthalpy for disordered Fe–Al alloys are in good agreement with experimental data or other theoretical calculations. This indicates that the present EAM potentials of Fe–Al system is suitable for atomistic simulations of structural and kinetic properties for the Fe–Al system.     

A reliable single parameter interatomic potential for argon

A simple semiempirical approximation previously proposed for the isotropic intermolecular forces between two closed shell systems is tested in detail for the argon-argon interaction. The potential is based on the knowledge of the first-order coulomb interaction energy, a suitably damped three term long range asymptotic expansion of the second order coulomb energy, and a semiempirical representation of the exchange interaction energy which contains one adjustable parameter. The single adjustable parameter can be reliably determined by fitting the second virial coefficient for argon in the 130–773 K temperature range with the long range interaction coefficients being constrained within the theoretical bounds specified by Tang, Norbeck and Certain. The reliability of the potential is compared with that of several literature potentials by comparing the theoretical predictions obtained from the potentials with experimental results for the second virial coefficient, viscosity, thermal conductivity and thermal diffusion ratios for dilute argon gas, and with spectroscopic data for the dimer, and with SCF calculations of the Ar-Ar potential at small interatomic separations. Our best potential predicts these properties with a precision as good as or better than other recent potentials which generally contain more adjustable parameters and/or involve more input data. The results confirm earlier work that suggested that the scheme tested is capable of yielding reliable isotropic potentials for the interaction of closed shell systems for 0·3 ? R/R_m ? ∞ where R_m is the intermolecular distance at the van der Waals minimum. The scheme appears to offer a method for obtaining reliable potentials while avoiding problems associated with optimizing many parameters with respect to fitting experimental constraints.     

Construction of embedded-atom-method interatomic potentials for alkaline metals （Li,Na,and K） by lattice inversion

The lattice-inversion embedded-atom-method interatomic potential developed previously by us is extended to alkaline metals including Li,Na,and K.It is found that considering interatomic interactions between neighboring atoms of an appropriate distance is a matter of great significance in constructing accurate embedded-atom-method interatomic potentials,especially for the prediction of surface energy.The lattice-inversion embedded-atom-method interatomic potentials for Li,Na,and K are successfully constructed by taking the fourth-neighbor atoms into consideration.These angular-independent potentials markedly promote the accuracy of predicted surface energies,which agree well with experimental results.In addition,the predicted structural stability,elastic constants,formation and migration energies of vacancy,and activation energy of vacancy diffusion are in good agreement with available experimental data and first-principles calculations,and the equilibrium condition is satisfied.     

Phonon dispersion for the bcc metals Mo and Cr

In this paper the phonon dispersion for the bcc metals Mo and Cr is calculated based on the pair potentials obtained from cohesive energies and the Slater-Kirkwood- type three body interaction. In the calculation of the pair potentials the Möbius transform in the number theory is used and the cohesive energy is evaluated by the LMTO method. The results show a good agreement with inelastic neutron scattering data and indicate that the three-body interaction is necessary to account for the phonon dispersion.     

Structural transition of sheared-liquid metal in quenching state

In this Letter, the effects of shear rate on structural properties of liquid Al in quenching process were investigated via molecular dynamics (MD) simulations based on the EAM potential. Analyses in internal energy and pair correlation functions (PCF) reveal an increasing structural transition temperature as the shear rate is enhanced in the liquid. Results of pair analysis indicate that for liquid Al under normal condition, face center cubic (FCC) structure is clearly detected upon cooling; while in sheared liquid, structural transition from FCC to body center cubic (BCC) at temperature of 800 K is manifested, leading to the dominance of BCC structural order at low temperatures.     

自由能方法与零压下Al的熔化温度

利用自由能方法的分子动力学模拟,计算了零压下Al的熔化温度.在计算液相自由能的过程中,采用勒纳-琼斯(LJ)液体作为参考系统,同时将计算结果与Mei和Davenport等人的计算结果进行了比较,计算结果表明：1)选用LJ参考系统使液相自由能的计算时间节省一半,并且不影响熔化温度的计算结果;2)采用不同的埋入原子势(EAM)的分子动力学模拟计算得到的熔化温度与实验值都存在偏差,而就金属Al而言,采用Cai等人的EAM势的熔化温度的计算结果比Mei和Davenport及Morris等人采用的势模型的结果略有改 关键词： 熔化温度 自由能方法 分子动力学模拟     

Static and thermophysical properties of chalcogenide crystals with NaCl structure

In the present communication an analysis of interionic potentials in fourteen chalcogenide crystals has been performed. This interionic potential has been used to predict the values of cohesive energy, isothermal bulk modulus and the pressure derivatives of bulk modulus in the solids under study. The many body interaction (MBI) effects have been taken into account within the framework of Hafemeister Zarht potential. Instead of using BM potential the Hafemeister-Zarht (HZ) type short range overlap potential has been considered between nearest as well as between next nearest neighbour ions. The short-range interactions, effective up to second neighbours are treated by considering the hardness parameter as an ionic property. The hardness parameter ρ_ij is evaluated using the data on overlap integrals. The results achieved in the present study are generally in good agreement with available experimental data. Values of cohesive energy, bulk modulus and its pressure derivatives calculated by previous investigators have also been shown for the sake of comparison.     

Molecular dynamics simulation of self-diffusion coefficients for liquid metals

The temperature-dependent coefficients of self-diffusion for liquid metals are simulated by molecular dynamics meth ods based on the embedded-atom-method （EAM） potential function. The simulated results show that a good inverse linear relation exists between the natural logarithm of self-diffusion coefficients and temperature, though the results in the litera ture vary somewhat, due to the employment of different potential functions. The estimated activation energy of liquid metals obtained by fitting the Arrhenius formula is close to the experimental data. The temperature-dependent shear-viscosities obtained from the Stokes-Einstein relation in conjunction with the results of molecular dynamics simulation are generally consistent with other values in the literature.     

Calculation of Vacancy Induced Properties in FCC Metals Using an Alternative Lattice Dynamical Model

Combining the conventional Kanzaki's lattice statics method with an alternative lattice dynamical model called decoupling transformation, a new approach is proposed to calculate the vacancy induced properties in FCC metals. This approach has the advantage over the traditional least-square fit approach in the way that all the model parameters are linearly independent and so it can avoid the problem of non-uniqueness in dynamical model parameters when interactions with more distant atoms are taken into account by fitting to the experimental phonon frequencies only. Numerical results on the lattice dynamical properties such as phonon dispersion curves and interatomic interaction force constants and the vacancy induced properties such as the atomic displacements, lattice relaxation energy, divacancy interaction energy and relaxation volume are obtained specifically for two similar FCC metals, to wit, Cu and Ni and are compared with previous calculations. It is concluded that the interatomic interactions up to the fourth nearest neighbours are already good enough for describing both the lattice vibrational and vacancy induced properties for both metals.     

Reconstructed (1 1 0) surfaces of FCC transition metals

The energies of the ideal, missing row (MR) and missing column (MC) (1 1 0) surfaces have been calculated by using modified embedded atom method (MEAM) for seven face centered cubic (FCC) transition metals Au, Pt, Ag, Pd, Rh, Cu and Ni. The results, that the MC reconstruction can not be formed for all metals, while the MR reconstruction can be formed naturally for Au and Pt, inductively for Ag, Pd, Rh and Cu and difficultly for Ni, are better than EAM calculated results in comparing with experimental results. In addition to the surface energy explanation, the results are also related to the surface topography and valence electron structure.