Molecular dynamics calculation of half-lives for thermal decay of Lennard-Jones clusters

Molecular dynamics has been used with a Lennard-Jones (6–12) potential in order to study the decay behavior of neutral Argon clusters containing between 12 and 14 atoms. The clusters were heated to temperatures well above their melting points and then tracked in time via molecular dynamics until evaporation of one or more atoms was observed. In each simulation, the mode of evaporation, energy released during evaporation, and cluster lifetime were recorded. Results from roughly 2000 simulation histories were combined in order to compute statistically significant values of cluster half-lives and decay energies. It was found that cluster half-life decreases with increasing energy and that for a given value of excess energy (defined asE=(E _tot ?E _gnd)/n), the 13 atom cluster is more stable against decay than clusters containing either 12 or 14 atoms. The dominant decay mechanism for all clusters was determined to be single atom emission.