Melting of sodium clusters

Abstract

Thermal stability properties and the melting-like transition of Na$_N$, $N$ = 13-147, clusters are studied through microcanonical molecular dynamics simulations. The metallic bonding in the sodium clusters is mimicked by a many-body Gupta potential based on the second moment approximation of a tight-binding Hamiltonian. The characteristics of the solid-to-liquid transition in the sodium clusters are analyzed by calculating physical quantities like caloric curves, heat capacities, and root-mean-square bond length fluctuations using simulation times of several nanoseconds. Distinct melting mechanisms are obtained for the sodium clusters in the size range investigated. The calculated melting temperatures show an irregular variation with the cluster size, in qualitative agreement with recent experimental results. However, the calculated melting point for the Na$_{55}$ cluster is about 40 $%$ lower than the experimental value.