  Molecular dynamics simulations of nanometer-sized lead clusters have been
performed using the Lim, Ong and Ercolessi glue potential (Surf. Sci. {\bf
269/270}, 1109 (1992)). The binding energies of clusters forming crystalline
(fcc), decahedron and icosahedron structures are compared, showing that fcc
cuboctahedra are the most energetically favoured of these polyhedral model
structures. However, simulations of the freezing of liquid droplets produced a
characteristic form of ``shaved'' icosahedron, in which atoms are absent at the
edges and apexes of the polyhedron. This arrangement is energetically favoured
for 600-4000 atom clusters. Larger clusters favour crystalline structures.
Indeed, simulated freezing of a 6525-atom liquid droplet produced an imperfect
fcc Wulff particle, containing a number of parallel stacking faults. The
effects of temperature on the preferred structure of crystalline clusters below
the melting point have been considered. The implications of these results for
the interpretation of experimental data is discussed.
