Andreas Dullweber, Matt P. Hodges and David J. Wales*
University Chemical Laboratories, Lensfield Road, Cambridge CB2 1EW, UK
We use a combination of Molecular Dynamics, Monte Carlo and geometry optimisation techniques to study benzene-Arn clusters for 1<=n<=19, with particular emphasis on BzAr19. In particular, we discuss the difficulties which arise in the accurate simulation of heterogeneous clusters due to problems of ergodicity and non-additive contributions to the energy. The sensitivity of minima, transition states and reaction pathways to parameters of the potential and the induction energy is also considered. An efficient integration scheme with adaptive step size due to Bulirsch and Stoer is employed in the Molecular Dynamics simulations. Both geometry optimisation and Molecular Dynamics are considered to evaluate the usefulness of the Jump-Walking Monte Carlo method proposed by Frantz, Freeman and Doll. This approach improves the ergodicity of canonical simulations using data from different temperatures which we achieved using multiple parallel runs. We then apply a multiple histogram method to calculate the relative number of states in phase space and various thermodynamic properties covering the full temperature range in the canonical and the microcanonical ensembles. The Monte Carlo and Molecular Dynamics simulations result in a consistent picture of transitions between minima and escape times. Finally, we consider the rate of side-crossing by Ar atoms in BzAr1 and BzAr19 and compare statistical theories with rates obtained from simulations.
© 1996 American Institute of Physics, The Journal of Chemical Physics, 47, 1996 (in press)
This versions was revised on the 20th October 1996.