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.

[1]  H. Reiss,et al.  Simulative determination of kinetic coefficients for nucleation rates , 2001 .

[2]  J. A. Barker,et al.  Theory and Monte Carlo simulation of physical clusters in the imperfect vapor , 1973 .

[3]  M. Kulmala,et al.  Dimers in nucleating vapors , 1998 .

[4]  E. Ruckenstein,et al.  Homogeneous nucleation in gases: A three‐dimensional Fokker–Planck equation for evaporation from clusters , 1991 .

[5]  B. C. Garrett,et al.  Dynamical nucleation theory: Calculation of condensation rate constants for small water clusters , 1999 .

[6]  J. Barrett The significance of cluster lifetime in nucleation theory , 2002 .

[7]  I. Ford Nucleation theorems, the statistical mechanics of molecular clusters, and a revision of classical nucleation theory , 1997 .

[8]  Dimo Kashchiev,et al.  Nucleation : basic theory with applications , 2000 .

[9]  D. Oxtoby,et al.  Nonclassical nucleation theory for the gas-liquid transition , 1988 .

[10]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[11]  H. Kramers Brownian motion in a field of force and the diffusion model of chemical reactions , 1940 .

[12]  H. Risken The Fokker-Planck equation : methods of solution and applications , 1985 .

[13]  M. Klein,et al.  Aggregation-volume-bias Monte Carlo simulations of vapor-liquid nucleation barriers for Lennard-Jonesium , 2001 .

[14]  I. Ford,et al.  A dynamical definition of quasibound molecular clusters , 2003 .

[15]  H. Risken Fokker-Planck Equation , 1984 .

[16]  P. Hänggi,et al.  Reaction-rate theory: fifty years after Kramers , 1990 .

[17]  E. Ruckenstein,et al.  A kinetic approach to the theory of nucleation in gases , 1991 .

[18]  Howard Reiss,et al.  A molecular theory of the homogeneous nucleation rate. I. Formulation and fundamental issues , 1999 .

[19]  J. Seinfeld,et al.  Direct evaluation of the equilibrium distribution of physical clusters by a grand canonical Monte Carlo simulation , 1998 .

[20]  Daan Frenkel,et al.  Computer simulation study of gas–liquid nucleation in a Lennard-Jones system , 1998 .

[21]  Fokker–Planck equations for nucleation processes revisited , 1998, cond-mat/9802299.

[22]  B. N. Hale Monte Carlo calculations of effective surface tension for small clusters , 1996 .

[23]  X. Zeng,et al.  Formation free energy of clusters in vapor-liquid nucleation: A Monte Carlo simulation study , 1999 .

[24]  R. Becker,et al.  Kinetische Behandlung der Keimbildung in übersättigten Dämpfen , 1935 .

[25]  R. Soto,et al.  Cluster birth–death processes in a vapor at equilibrium , 1999 .

[26]  Frank H. Stillinger,et al.  Rigorous Basis of the Frenkel-Band Theory of Association Equilibrium , 1963 .

[27]  B. C. Garrett,et al.  Variational transition state theory of vapor phase nucleation , 1999 .

[28]  B. C. Garrett,et al.  Dynamical benchmarks of the nucleation kinetics of water , 2002 .

[29]  New criteria for cluster identification in continuum systems , 2001, cond-mat/0110518.