The thermodynamic data, such as the formation enthalpies of disordered solid solutions and intermetallic compounds, and the lattice parameters in the RE–Au (RE = Sc, Y, Ce, Pr, Nd, Gd, Tb, Dy, Ho, Er) systems are calculated with the modified analytical embedded atom method (EAM). This work has been undertaken to investigate systemically in RE–Au alloying energies, and to augment available calorimetric data for enthalpies of formation in support of the development of accurate multicomponent thermodynamics databases for these technologically interesting systems. The accuracy of our calculations is assessed through comparisons with the experimental measurements and the theoretical results from the first-principles VASP as well as Miedema’s theory. The composition dependence of the heats of formation for these 10 binary systems is similar, and there is a minimum at 50 at.% Au for REAu (CsCl-type) ordered phase. In all 10 binary systems, the calculated zero-temperature intermetallic formation energies generally agree well with the calorimetric data obtained by the direct reaction synthesis and/or VASP results except for REAu2 (MoSi2type) phase. For the intermetallic phase, the calculated EAM and VASP zero-temperature lattice parameters agree well with the experimental data at ambient temperature. The largest discrepancy between the EAM calculated lattice parameters and experimental data [K. Fitzner,W.G. Jung, O.J. Kleppa, Metall. Trans. A 22 (1991) 1103 [16]; K. Fitzner, O.J. Kleppa, Metall. Trans. A 24 (1993) 1827 [17]; K. Fitzner, O.J. Kleppa, Metall. Trans. A 25 (1994) 1495 [18]; R. Ferro, G. Borzone, N. Parodi, J. Alloys. Compd. 321 (2001) 248 [19]] is approximately 10-18% for some intermediate phases, such as REAu (CrB-type) and RE3Au4 (Pu3Pd4-type).

Analytic modified embedded atom method type many-body potentials have been constructed for alkali metals. The bcc lattice is shown to be energetically most stable when compared with fcc and ideal hcp structures. The phonon dispersion curves, density of states, Debye temperature, heat capacity, surface energy and thermal expansion properties of these metals are calculated, which are comparable to experimental and other theoretical results. The properties of point defects, such as vacancy, divacancy, and self-interstitial atoms, have also been evaluated with these potentials, and the diffusion mechanism in alkali metals are discussed.

Using molecular dynamics simulations and a modified analytic embedded-atom method (MAEAM), the relaxations and vibrations of Ni(977) surface are studied in the temperature range of 100–1500K, The calculated results for the temperature dependence of phonon frequencies, line-width and the mean square amplitude show that the anharmonic effects are small for temperatures up to 900K. The calculated interlayer separation for the first and second surface decreases, and the distance along the x direction increases between the topmost and second surface, with increasing temperature. In addition, the calculated layer structure factor indicates that the Ni(977) surface is well ordered and does not premelt up to a temperature of 1700K.

In the present paper, the surface and size effects on the alloying ability and phase stability of immiscible alloy nanoparticles have been studied with calculating the heats of formation of Au-Pt alloy nanoparticles from the single element nanoparticles of their constituents (Au and Pt) with a simple thermodynamic model and an analytic embedded atom method. The results indicated that, besides the similar compositional dependence of heat of formation as in bulk alloys, the heat of formation of alloy nanoparticles exhibits notable size-dependence, and there exists a competition between size effect and compositional effect on the heat of formation of immiscible system. Contrary to the positive heat of formation for bulkimmiscible alloys, a negative heat of formation may be obtained for the alloy nanoparticles with a small size or dilute solute component, which implies a promotion of the alloying ability and phase stability of immiscible system on a nanoscale. The surface segregation results in an extension of the size range of particles with a negative heat of formation. The molecular dynamics simulations have indicated that the structurally and compositionally homogeneous Au Pt nanoparticles tend to form a core-shell structure with temperature increasing.

Using the modified analytic embedded atom method and molecular dynamics, the binding energies and their second order finite differences stability functionsd of icosahedral Ni clusters with shell and subshell periodicity are studied in detail via atomic evolution. The results exhibit shell and sub shell structures of the clusters with atoms from 147 to 250 000, and the atomic numbers corresponding to shell or subshell structures are in good agreement with the experimental magic numbers obtained in time-of-flight mass spectra of threshold photoionization, and Martin’s theoretical proposition of progressive formation of atomic umbrellas. Clusters with size from 147 to 561 atoms are energetically investigated via one-by-one atomic evolution and their magic numbers are theoretically proved. For medium-size Ni clusters with 561 to 2057 atoms, the prediction of magic numbers with atomic numbers is performed on the basis of umbrellalike subshell growth in near face-edge-vertex order. The similarity of the energy curves makes it possible to extend the prediction to even larger Ni nanoclusters in hierarchical Mackay icosahedral configurations.