In this paper, the melting behaviors of Rh–Ag–Au nanoalloys are investigated with MD simulation. For Rh–Ag–Au nanoalloys, icosahedron structure was considered. The local optimizations of Rh–Ag–Au nanoalloys were carried out with the BH algorithm. The interatomic interactions were modeled with the Gupta potential. The local optimization results of Rh–Ag–Au nanoalloys show that Au and Ag atoms prefer to locate on the surface, and Rh atoms prefer to locate in the inner shells. The bond order parameter result is compatible with the excess energy analysis. It is noted that structures with more Ag–Au bonds are more energetically stable. Caloric curve, heat capacity, Lindemann index, and RMSD methods were used for estimating the melting temperatures of Rh–Ag–Au nanoalloys. According to the simulation results, melting temperatures depend on the composition. Also, it is discovered that nanoalloys are generally melting in two stages. Surface melting of the third shell is occupied by Ag and Au atoms, and then homogeneous melting of the inner shells is occupied by Rh atoms. It is found that the difference between surface melting temperatures and homogeneous melting temperatures in Ag-poor compositions is more significant than that of Ag-rich nanoalloys. In addition, the melting temperatures of the nanoalloys are found to be increased as the size of nanoalloys increases.
Theoretical study of the structures of bimetallic Ag–Au and Cu–Au clusters up to 108 atoms