Molecular dynamics and DFT study of 38-atom coinage metal clusters.

[Display omitted] • Melting temperature of metal clusters is estimated from molecular dynamics simulations. • Structural changes are monitored as a function of temperature. • Trimetallic clusters exhibit the highest melting temperature. • DFT calculations predict metallic behavior of Cu 36 Au 1 Ag 1 cluster isomers. The thermal behavior of 38-atom mono-, bi-, and trimetallic clusters consisting of Cu, Ag, and Au atoms, is analyzed employing molecular dynamics simulations and DFT calculations for selected cluster compositions. Low-energy structures were singled out to perform NVT molecular dynamics simulations at several temperatures, using the Andersen thermostat for temperature control. The caloric curve is used to estimate the melting temperature and the specific heat. The pair distribution function g(r) of the solid and liquid-phase clusters is examined at different temperatures. When comparing the estimated melting points (T m) among the monatomic clusters, the order becomes T m Cu 38 > T m Ag 38 > T m Au 38 . For bimetallic clusters, an increase of T m is observed for Cu-Au compared to their monatomic counterparts, while the opposite occurs for Cu-Ag clusters. For trimetallic clusters, two low-energy isomers of the Cu 36 Ag 1 Au 1 cluster are investigated. In this case, T m is estimated to be 475 K, for the two isomers with the lowest-energy and second-to-lowest energy, respectively. For all the clusters studied, the pair distribution function g(r) shows that the first peak position is not shifted as an effect of temperature and its maximum value varies with composition, while the second peak essentially vanishes upon melting. The common-neighbor analysis (CNA) technique is used to analyze the local structural changes for the trimetallic clusters, again demonstrating a clear structural change upon melting. The HOMO-LUMO energy gap indicates that the trimetallic isomers' behavior is metallic, while the average binding energy show these clusters' energetic stability to be similar. [ABSTRACT FROM AUTHOR]