Van der Waals Component. The second component of the non-bonded interactions between the QM and MM molecules is the so called van der Waals component. This component accounts for the repulsive steric and the attractive dispersion interactions between molecules. In computer simulations, van der Waals interactions between molecules are almost exclusively represented as a sum of two body potentials between the atoms within the molecules. The general behaviour of a typical van der Waals potential used in computer simulations is shown in Figure 7.5 where the interaction energy is plotted against the distance, r, between two atomic centers. The distinctive feature of the potential is the attractive well that exists at a distance Ro and with a depth Do. The steeply repulsive region of the potential at short distances, ro, models the steric or Pauli repulsion due to the overlap of electron clouds. The long range tail in the van der Waals potential at r>Ro accounts for the attractive 'London' dispersion interactions due to the instantaneous correlation of electron clouds surrounding the atoms.
Figure 7.5 Example of a van der Waals non-bonded potential between atoms.
Two common functions used to represent the van der Waals potential in empirical force fields are the Lennard-Jones 12-6 potential (Equation 7-2) and the Buckingham exponential-6 potential(Equation 7-3). These functions reproduce the general behaviour of real interactions as described above, in an empirical way. For both functions, the long ranged attractive component which models the dispersion interactions, decays with a r6 relationship. For the short-ranged steric component, the two potentials differ in their representation. The more common Lennard-Jones potentials utilizes a R12 function while the Buckingham potential uses a exponential relationship. Although, the Lennard-Jones potential is often considered more computationally efficient, its representation of the repulsive wall is sometimes too steep. Oppositely, whereas the Buckingham potential is more computationally demanding, it also more correctly models the exponential behaviour of the steric wall.
(7-2)
=12.0§ (7-3)
The van der Waals potential is characterized by the well position, Ro, and the well depth, Do. In a typical molecular mechanics force field, values are assigned to each atom type and are fit to generate results consistent with the rest of the force field. The atom type is generally defined in terms of the atomic element, its hybridization and the functional group it is contained within. Table 7.1 shows some values from the AMBER9577 molecular mechanics force field. When the van der Waals potential is required for two atoms of differing type, say between CT and O, they are generated from the two homonuclear parameters using some kind of combination rule such at that expressed in Equations 7-4 and 7-5.