QRC: Quick Reaction Coordinate
Cyclohexane example
J M Goodman and M A Silva
Tetrahedron Letters 2003, 44, 82338236.
Cyclohexane has two minima in its potential energy surface  a chair and
a twist boat. It also has conformational transition states. This example
examines the cyclohexane potential energy surface using RHF/631G** as
implemented in Jaguar.
The symmetrical boat conformation has one imaginary frequency. The QRC follows
this direction to the twistboat minimum. The flat hexagon of carbon atoms has
three imaginary frequencies. The lowest of these is traced by the QRC to the
chair conformation, and the second and third (which are degenerate) lead to the
twist boat.
Procedure
First, the transition states were identified and frequency calculations were
run Jaguar. The output files were:
flatfreq.out
boatfreq.out
FreqQRC was used first:
FreqQRC flatfreq.out
This creates a new file flatfreq.fdt which
contains the geometry of the molecule, the changes found from the frequencies,
and the two distorted geometries, one adding and one subtracting the distortions.
The new geometries were cut and pasted into new Jaguar input files, and minimised.
The outputs were flatfreq_add.out and flatfreq_sub.out.
Now, FindQRC was used:
FindQRC flatfreq_add.out flatfreq.fdt
This creates a new file flatfreq_add.qrc.
The first column is the energy of each structure, and the next columns are different
ways of measuring the distance between adjacent structures. The second column is
massweighted rms distance from the starting structure, and we have found this to be
the most effective. The third column gives the sum of the distances between adjacent
structures, and is similar but slightly larger. The fourth and fifth columns are
nonmassweighted versions of the second and third.
The figure above shows, in blue, the effect of plotting the energies against the
distances in the second column (massweighted distances from the initial structure).
If the corresponding calculation is done on flatfreq_sub, the same result is obtained
(shown in the figure with negative distances) because this system is symmetrical. In
general, adding and substracting the frequencies may give different results. The structure
at zerodistance is the transition state read from the flatfreq.fdt file, and
this is why this file is mentioned after the FindQRC command.
Variations
FreqQRC: Adding the whole values for the frequencies can distort the structure too
much. This occurred for the symmetrical boat structure, for which the structures generated
by FreqQRC had higher energies than the transition structures. This suggests they have
been distorted beyond the region that the harmonic frequencies describe. In this case, they
still minimise to the twist boat structures, but a smaller distortion would be preferable.
This is achieved by using the command:
FreqQRC boatfreq.out 0.2
This multiplies each of the distortions by 0.2 before adding and substracting them.
FreqQRC boatfreq.out 1.0 is equivalent to FreqQRC boatfreq.out,
as the default is to use the distortions calculated by Jaguar directly.
Some structures, including flatfreq.out have more than one imaginary frequency.
By default, FreqQRC uses the first frequency, but higher frequencies may be
selected using the f flag. For example:
FreqQRC flatfreq.out f2
generates a flatfreq.fdt file containing distortions based on the second lowest frequency
reported by Jaguar. In this example, this is also an imaginary frequency, and the results of
following the FindQRC procedure are also shown in the figure, in purple. The lowest
frequency leads to the chair conformation, and the second lowest to the twist boat.
FindQRC:
The usual form for the FindQRC command includes an .fdt file:
FindQRC flatfreq_add.out flatfreq.fdt
However, the .fdt file can be omitted, and the calculation reports distances
from the first distorted structure.
By default, FindQRC omits points for which the energy goes up as the minimisation
progresses. Occasionally it can be useful to visualise these points too, and this can be
done using the u flag:
FindQRC flatfreq_add.out flatfreq.fdt u
This has been done for the lefthand section of the yellow line in the figure. This line
would be a reflection of the righthand half if it had been generated without the u flag.
© 2003 J M Goodman, Cambridge
