2.1.12 1D and 2D Energy vs. Torsion Plots

For studying molecular structures and their energies in detail, it is possible to calculate plots that show how energy is dependent on torsion angles of a molecule. We will get 1D or 2D plots depending on whether we use 1 or 2 torsion angles as variables.

We use 3-chloro-1-pentene as our example molecule:

example case

We have selected from the popup-menu "Render->Rendering Mode->Label Index" to make the atom index numbering visible, and also "Render->Views->Attach/Detach This View" in order to make the graphics view and the plot view visible simultaneously.

Now we need to define the bonds (or more precisely, the torsions) which we wish to rotate. For the first torsion we select atoms 2-1-0-3, and for the second torsion we select atoms 1-0-3-9. The direction in which we select the torsion is not relevant, but the sequence itself must be correct.

Now let's create a 1D plot that displays the energy vs. the first torsion, using AM1 to compute final energies. From the pop-up menu we select "Compute->Plot Energy vs. 1 Torsion Angle...". The following default command string is then displayed:

make_plot1 A B C D 36 0.0 360.0 250

The letters A-D represent the atom indices, that we have to replace using our values 2-1-0-3 determined above. If we wish to calculate the final energies with AM1, we also append the values "2 0" to the end (2 tells that the method is AM1, and 0 tells the total charge to be used in the calculation):

make_plot1  2 1 0 3  36 0.0 360.0 250  2 0

The resulting 1D energy plot is the following:

the 1D plot

If the plot is clicked using the Orbit-XY mouse tool, it will display the structure in graphics view window and print out the energy at the console window.

Now let's try the 2D plot, using both first and second torsion angles as variables. Again we select from the pop-up menu "Compute->Plot Energy vs. 2 Torsion Angle...", and type in the following changes to the default command string:

make_plot2 A B C D 36 0.0 360.0 I J K L 36 0.0 360.0 250

Since our torsions are 2-1-0-3 and 1-0-3-9 and we wish to use AM1 we will make the following modifications:

make_plot2  2 1 0 3  36 0.0 360.0  1 0 3 9  36 0.0 360.0 250  2 0

The resulting 2D energy plot is the following:

the 2D plot

In the above plot, X-axis is the first torsion angle, and the Y-axis is the second torsion angle. The colours of the plot indicate the energy; green colour means low energy and red/white colour means high energy.

Just like in the 1D plot, the structure is displayed if the the plot is clicked using the Orbit-XY mouse tool. In addition, the Zoom mouse tool changes the colour scaling of the plot, making it easier to find the desired details from the plot. The following is the same plot as above, but using colour scaling that better shows the lowest-energy conformations:

Following is one of the most probable conformations of 3-chloro-1-pentene according to the above plot:

the minimum-energy structure


Return to index