Bond stretching

As you know, a useful function for modelling typical potential energy curves for bond stretching is the Morse potential. This is usually not used in molecular mechanics (it involves 3 parameters for every bond, and is also time-consuming to calculate). Instead, the simpler harmonic form (above) is used, requiring two parameters (the force constant, and the reference bond length l0). With the harmonic form, the energy of bond stretching can be calculated quickly, and the behaviour close to the minimum is given reasonably well. However, modelling bonds as springs means that they can't dissociate - this is one of the reasons that standard molecular mechanics methods cannot be used to model chemical reactions. Typical values of force constants are of the order of several hundred kcal mol-1Å-2. 

 

Angle bending

The energy bending of valence angles is also calculated using a harmonic term. Less energy is needed to distort an angle away from its equilibrium value, so force constants for angle bending are smaller.

 

Rotation around bonds: torsion angle energy

Torsion terms are very important in determining the structure of 'organic' type molecules. The energy for rotation of torsions are represented by series of cosine functions, i.e. the energy for a torsion angle w is given by:

                   (2)

 

Typically 1-3 cosine terms will be used for every torsion. Vn is often called the barrier height, but this is not a very good description - it can even take negative values - and each torsion may be associated with several Vn. Also, other terms (e.g. non-bonded 1-4 interactions) will contribute to the barrier height. n is the multiplicity, giving the number of minima as the torsion is rotated through 360º. g is the phase factor - this gives where the cosine function has its minimum.

 

'Improper' torsions

These are needed, for example, to keep the 3 atoms bonded to an sp2 carbon in the same plane. Molecular mechanics potential function do not include the effects of p-bonding, for example, and so to model a carbonyl group correctly, an additional improper torsion term is needed. This is described as an improper torsion because the atoms involved are not bonded together in the normal A-B-C-D order.

For example, for the carbonyl fragment shown, an improper torsion could be defined as 1-4-3-2.

 

A cosine potential is then used to restrain the improper torsion angle to 0º or 180º:

 

                         (3)

 

Cross terms

In reality, individual bond stretches, angle bends and torsion twists are not independent of one another. For example, as a valence angle bends, the bonds involved are likely to change at the same time. Another example is shown in this figure:

With rotation around the bond, the bond should lengthen, (stretching of the bond leads to a reduction in the barrier, facilitating rotation around the bond). This could be represented by a stretch-torsion cross term (for the bond length, l (e.g. B-C in the figure), and the torsion angle w (e.g. A-B-C-D in the figure):

 

                      

To represent the coupling between angle bending and bond stretching, Urey-Bradley terms may be used, effectively applying a harmonic restraint on the distance, r, between atoms 1,3 with respect to each other.

 

                           

 

Non-bonded Interactions

Van der Waals and electrostatic forces play vital roles in determining the structures of molecules, and their interactions. Electrostatic forces act over large distances. Non-bonded interactions play an important role not only in the interaction of separate molecules, but also within molecules.

 

van der Waals interactions

Deviations from ideal gas behaviour (e.g. in noble gases) were studied in detail by van der Waals, and so the interactions giving rise to these deviations (not due to the effects of electrostatics) are referred to as van der Waals interactions. These can be thought of as being made up of an attractive (longer range) and a (short range) repulsive component. The attractive contribution is due to dispersion forces, created by instantaneous fluctuations in electronic distribution, giving rise to an instantaneous dipole in a molecule or atom, which in turn creates a dipole in neighbouring molecules or atoms. Dispersion interactions were analysed in the early days of molecular quantum mechanics by London, and so are sometimes referred to as London forces.

 

Next: more on non-bonded interactions