Molecular cluster decay viewed as escape from a potential of mean force

We show that evaporation from a quasistable molecular cluster may be treated as a kinetic problem involving the stochastically driven escape of a molecule from a potential of mean force. We derive expressions for the decay rate, and a relationship between the depth of the potential and the change in system free energy upon loss of a molecule from the cluster. This establishes a connection between kinetic and thermodynamic treatments of evaporation, but also reveals differences in the prefactor in the rate expression. We perform constant energy molecular dynamics simulations of cluster dynamics to calculate potentials of mean force, friction coefficients and effective temperatures for use in the kinetic analysis, and to compare the results with the directly observed escape rates. We also use the simulations to estimate the escape rates by a probabilistic analysis. It is much more efficient to calculate the decay rate by the methods we have developed than it is to monitor escape directly, making these approaches potentially useful for the assessment of molecular cluster stability.