Potential functionsΒΆ

As mentioned earlier, all interaction components are represented by some mathematical functional forms. While the precise functions and their implementations depend on the FF schemes, the force field model or the total energy of the system generally can be expressed as follows:

\[V_{total} = V_{bonded} + V_{cross-terms} + V_{non-bonded}\]

where

\[V_{bonded} = V_{bond} + V_{angle} + V_{torsion} + V_{improper} + V_{UB}\]
\[V_{cross-terms} = V_{bond1,bond2} + V_{bond,angle} + \cdots\]
\[V_{non-bonded} = V_{vdW} + V_{Coulombic}\]
../../_images/Orange_bar2.png

Below is a list of typical functional forms representing a range of interaction components.

  1. Bond interaction, \(V_{bond}\)

    Harmonic bond interaction.

    where \(b_0\) represents the equilibrium bond length and \(b\) is the bond distance (at a given moment during the MD simulation).

  2. Angle interaction, \(V_{angle}\)

    Harmonic angle interaction.

    where \(\theta_0\) represents the equilibrium angle.

  3. Dihedral interaction, \(V_{torsion}\)

    Fourier cosine dihedral interaction.

    where \(k\) is the dihedral constant, or torsion barrier, \(n\) is the periodicity and \(\delta\) is the phase shift.

  4. Cross terms, for example bond-angle, \(V_{bond, angle}\)

    Fourier cosine bond-angle interaction.
  5. Electrostatic, \(V_{Coulombic}\)

    Coulombic electrostatic interaction.
  6. Dispersive, \(V_{vdW}\)

    Lennard-Jones 12-6 van der Waals interaction.

^ GO TO TOP ^