We find that the potential energy distribution of atoms in clusters can consistently explain many important phenomena related to the phase changes of clusters, such as the nonmonotonic variation of melting temperature with size, the dependence of melting, boiling, and sublimation temperatures on the interatomic potentials, the existence of a surface-melted phase, and the absence of a premelting peak in heat capacity curves. We also find a new type of premelting mechanism in double icosahedral Pd 19 clusters, where one of the two internal atoms escapes to the surface at the premelting temperature. DOI: 10.1103/PhysRevLett.86.999 PACS numbers: 36.40.Ei, 61.46.+w, 64.70.Dv With the possibility that metal clusters could be used as building blocks for new nanomaterials [1], the understanding of the structures and the thermal stabilities of metal clusters has become important. Recent experimental measurements of the heat capacities for unsupported metal clusters by Schmidt and co-workers [2,3] have given a strong motivation for the atomistic understanding of the melting of metal clusters. It has been observed both in experiments [4] and in theoretical simulations [5] that clusters melt at temperatures lower than those of bulk melting due to the high proportion of surface atoms with lower binding energies. As the cluster size decreases, the melting temperature monotonically decreases [4,5]. However, when the cluster size is small enough (# 200 atoms), the melting temperature does not vary monotonically with the size of cluster [3,6,7].Multimodal heat capacity curves observed in small clusters indicate that thermal energy can be dissipated by a "premelting" mechanism before overall melting [7]. Premelting has been attributed to surface melting [8], partial melting [9], orientational disordering in molecular cluster [10], and isomerization [7,11]. It has been observed that clusters can dissipate thermal energy even by boiling [12] and sublimation [13]. In some clusters, the premelting peak in the heat capacity curve is not observed [9,[14][15][16].Previously published simulations of cluster melting have used either Monte Carlo (MC) methods at constant temperature (canonical ensemble) or molecular dynamics (MD) at constant energy (microcanonical ensemble). There has been some controversy concerning the dependence of the simulated thermodynamic properties on the ensemble until the work by Calvo and Labastie [17]. These authors [17] showed that if the sampling of the phase space is accurate enough (a large enough number of MC configurations and MD trajectories), the thermodynamic results are identical in both ensembles because the thermodynamics is determined by the configurational density of states [7,8,17,18].In this Letter, we show that the potential energy distribution of atoms in the clusters can consistently explain many of the important phenomena which occur during phase changes of small clusters, such as the nonmonotonic variation of melting temperature with the size of clusters [3,6,7], the dependence of meltin...