next up previous
Next: BzAr Up: Molecular Dynamics Results Previous: Molecular Dynamics Results

BzAr tex2html_wrap_inline1451

The PES of BzAr tex2html_wrap_inline1451 is very simple, there being only one minimum and one transition state. We use two approaches to estimate the rate at which the Ar atom crosses the plane of the benzene molecule, namely MD and the methods developed by Rice and Ramsperger, and Kassel [50] (RRK).

For MD, we use 43ns trajectories with total energy between the energy of the transition state and the dissociation energy. Our MD results for BzAr tex2html_wrap_inline1413 use the SCM parameters only but here we present, for comparison, those obtained with OBJ as well. At regular time intervals we analyse the Cartesian coordinates for all atoms in the system to determine which side of the ring the Ar atom is on and hence count the number of times the Ar atom crosses the ring. This leads directly to the side-crossing rate which converges well on the given time-scale.

RRK theory assumes that the kinetic energy is accessible to all the vibrational normal modes in equal proportion. It can then be shown that the rate constant as a function of energy is

  equation526

where tex2html_wrap_inline2519 is the order of the point group; { tex2html_wrap_inline2521 }, the normal mode frequencies; E, the energy measured from zero at the minimum; tex2html_wrap_inline1401 , the energy of the transition state and s is the number of vibrational degrees of freedom in the system, i.e. 3 for BzAr tex2html_wrap_inline1451 with the benzene molecule considered rigid. The unlabelled quantities refer to the minimum; those labelled with tex2html_wrap_inline2531 , to the transition state.

The results are plotted in Fig. 15 and show that the RRK results are always within one order of magnitude of those calculated from MD. The RRK results are translated such that the energy is zero at dissociation. Both the behaviour at the transition state energies and close to dissociation are in quite good agreement.

   figure542

Figure 15: Comparison of BzAr tex2html_wrap_inline1451 ring-crossing rates from MD and RRK theory for both sets of Lennard-Jones parameters. All energies are measured relative to dissociation at zero and the transition state energies, tex2html_wrap_inline1401 , are plotted as vertical lines to highlight the behaviour of the rate data.

We present these results for comparison with BzAr tex2html_wrap_inline1413 for which we estimate barriers for one-sided to two-sided transitions from the MD trajectories. Calculating ring-crossing rates may not be as useful, as we have the possibility of multiple bridging atoms, and using RRK theory is tedious because of the number of transition states we would have to consider.


next up previous
Next: BzAr Up: Molecular Dynamics Results Previous: Molecular Dynamics Results

Matt Hodges & Andreas Dullweber
Fri Oct 20 09:28:06 BST 1996