The analytic embedded atom method (EAM) type many-body potentials of hcp rare earth metals (Dy, Er, Gd, Ho, Nd, Pr, and Tb) have been constructed. The hcp lattice is shown to be energetically most stable when compared with the fcc and bcc structure, and the hcp lattice with ideal c/a. The mechanical stability of the corresponding hcp lattice with respect to large change of density and c/a ratio is examined. The phonon spectra, stacking fault and surface energy are calculated. The activation energy for vacancy diffusion in these metals has been calculated and the most possible diffusion paths are predicted. Finally, the self-interstitial atom (SIA) formation energy and volume have been evaluated for eight possible sites. This calculation suggests that the crowdion and basal split are the most stable configurations. The SIA formation energy increases linearly with the increase of the melting temperature.

Analytic embedded-atom method (EAM) interatomic potentials for body-centered cubic (bcc) transition metals have been constructed. The total energy is regarded as consisting of a pair potential part, an embedding interaction part and a modification term. All these parts are analytic functions of the atomic separations only and are fitted to various bulk properties, such as cohesive energy, vacancy formation energy, lattice constant and elastic constants. The present potentials are shown to predict a more realistic pressure–volume relationship so that interactions at separations smaller than that of the first-nearest neighbors can be treated. Interstitial formation energies are calculated for various configurations, using quenched molecular dynamics. Vacancy, surface and phonon spectra have been also studied andsatisfactory agreement with available experimental data has been found.

Isothermal grain growth behaviors, including the effect of temperature and mean grain size, of nanocrystalline Ag are investigated using the molecular dynamics simulations with an analytical embedded atom potential. The atomistic simulations provide a level of atomic details on microscopic mechanism and process of grain growth that still remains inaccessible to experiments. The small grain size and high temperature accelerate the grain growth, it is the same as the conventional polycrystalline materials. The grain growth processes of nanocrystalline Ag are well characterized by a power-law growth curve, followed by a linear relaxation stage. Beside grain boundary migration and grain rotation mechanisms, the dislocations(or stacking faults)serve as the intermediate role in the grain growth process.

The surface premelting, melting behavior, melting competition, and structural transition of shell-closed icosahedral ICO and cuboctahedral CUB nickel clusters with atoms from 309 to 2057 were discussed extensively by using quantitative caloric curves based on the modified analytic embedded atom method and molecular dynamics, qualitative three-dimensional structural visualization of symmetric truncation, and the radial number distribution function. These studies reveal that smaller clusters melt at lower temperatures and a solid-to-solid structural transition occurs from CUB to ICO structure during melting process. The shell-closed ICO clusters could only be preferred until 923 atoms at temperatures no higher than 1380 K, which is in agreement with the experiments. The melting temperature of larger clusters would depend on their starting structures, which can be attributed to surface premelting.

By extrapolating the mean grain size of nanocrystal to an infinitesimal value, an amorphous phase has been obtained from the Voronoi construction. The molecular dynamics simulations indicated that for nanocrystal, the grain size variation of melting temperature exhibits two characteristic regions. As mean grain size above about 4 nm for Ag, the melting temperatures decrease with decreasing grain size. However, with grain size further shrinking, the melting temperatures almost keep a constant. This is because the dominant factor on the melting temperature of nanocrystal shifts from grain phase to grain boundary. As a result of fundamental difference in structure, the amorphous phase has a much lower solid-to-liquid transformation temperature than that of nanocrystal.