Density functional molecular dynamical simulations have been performed on Ga17 and Ga13 clusters to understand the recently observed higher-than-bulk melting temperatures in small gallium clusters [Phys. Rev. Lett. 91, 215508 (2003)]]. The specific-heat curve, calculated with the multiple-histogram technique, shows the melting temperature to be well above the bulk melting point of 303 K, viz., ...