Interaction bonded

The bonded interactions are needed to model the energetic penalty that will result fi-om two covalently bonded atoms moving too close or too far away from one another. The most common functional form that is used to model the bond bending interactions is that of a harmonic term  

Another more quantitative strategy that has been used to determine whether a pair of atoms is bonded or not is based on the bond critical point properties of the electron density distribution. According to Bader (1998), two atoms are indicated to be bonded if 

---

There is a fair amount of work reported with films at the mercury-air interface. Rice and co-workers [107] used grazing incidence x-ray diffraction to determine that a crystalline stearic acid monolayer induces order in the Hg substrate. Quinone derivatives spread at the mercury-n-hexane interface form crystalline structures governed primarily by hydrogen bonding interactions [108]. 

As the tip is brought towards the surface, there are several forces acting on it. Firstly, there is the spring force due to die cantilever, F, which is given by = -Icz. Secondly, there are the sample forces, which, in the case of AFM, may comprise any number of interactions including (generally attractive) van der Waals forces, chemical bonding interactions, meniscus forces or Bom ( hard-sphere ) repulsion forces. The total force... 

In LN, the bonded interactions are treated by the approximate linearization, and the local nonbonded interactions, as well as the nonlocal interactions, are treated by constant extrapolation over longer intervals Atm and At, respectively). We define the integers fci,fc2 > 1 by their relation to the different timesteps as Atm — At and At = 2 Atm- This extrapolation as used in LN contrasts the modern impulse MTS methods which only add the contribution of the slow forces at the time of their evaluation. The impulse treatment makes the methods symplectic, but limits the outermost timestep due to resonance (see figures comparing LN to impulse-MTS behavior as the outer timestep is increased in [88]). In fact, the early versions of MTS methods for MD relied on extrapolation and were abandoned because of a notable energy drift. This drift is avoided by the phenomenological, stochastic terms in LN. 

In addition, the non-bonded forces can be divided into several regions according to pair distances. The near region is normally more important than the distant region because the non-bonded forces decay with distance. Since most of the CPU time in a MD simulation is spent in the calculation of these non-bonded interactions, the separation in pair distance results in valuable speedups. Using a 3-fold distance split, the non-bonded forces are separated in 3 regions near, medium, and fax distance zones. Thus, the Liouville operator can be express as a sum of five terms... 

In the next step, the suggested models are translated into 3D space by subsequently combining the templates. Again, each model is assessed and ranked according to various structural criteria, such as the geometric fit of the 3D templates and non-bonding interactions (steric clashes). If none of the solu-... 

Figure 2-99. Elimination of non-bonded interactions (close contacts).

To calculate the bonded interaction of two atoms, a Morse function is often used. It has the form described by Eq, (19). 

Figure 7-9. Variation of the potential energy of the bonded interaction of two atoms with the distance between them. The solid line comes close to the experimental situation by using a Morse function the broken line represents the approximation by a harmonic potential.

Vn is often called the barrier of rotation. This is intuitive but misleading, because the exact energetic barrier of a particular rotation is the sum of all V components and other non-bonding interactions with the atoms under consideration. The multiplicity n gives the number of minima of the function during a 360° rotation of the dihedral angle o). The phase y defines the exact position of the minima. 

The origin of a torsional barrier can be studied best in simple cases like ethane. Here, rotation about the central carbon-carbon bond results in three staggered and three eclipsed stationary points on the potential energy surface, at least when symmetry considerations are not taken into account. Quantum mechanically, the barrier of rotation is explained by anti-bonding interactions between the hydrogens attached to different carbon atoms. These interactions are small when the conformation of ethane is staggered, and reach a maximum value when the molecule approaches an eclipsed geometry. 

It is noteworthy that it is not obligatory to use a torsional potential within a PEF. Depending on the parameterization, it is also possible to represent the torsional barrier by non-bonding interactions between the atoms separated by three bonds. In fact, torsional potentials and non-bonding 1,4-interactions are in a close relationship. This is one reason why force fields like AMBER downscale the 1,4-non-bonded Coulomb and van der Waals interactions. 

N is the number of point charges within the molecule and Sq is the dielectric permittivity of the vacuum. This form is used especially in force fields like AMBER and CHARMM for proteins. As already mentioned, Coulombic 1,4-non-bonded interactions interfere with 1,4-torsional potentials and are therefore scaled (e.g., by 1 1.2 in AMBER). Please be aware that Coulombic interactions, unlike the bonded contributions to the PEF presented above, are not limited to a single molecule. If the system under consideration contains more than one molecule (like a peptide in a box of water), non-bonded interactions have to be calculated between the molecules, too. This principle also holds for the non-bonded van der Waals interactions, which are discussed in Section 7.2.3.6. 

The raie gas atoms reveal through their deviation from ideal gas behavior that electrostatics alone cannot account for all non-bonded interactions, because all multipole moments are zero. Therefore, no dipole-dipole or dipole-induced dipole interactions are possible. Van der Waals first described the forces that give rise to such deviations from the expected behavior. This type of interaction between two atoms can be formulated by a Lennaid-Jones [12-6] function Eq. (27)). 

A second idea to save computational time addresses the fact that hydrogen atoms, when involved in a chemical bond, show the fastest motions in a molecule. If they have to be reproduced by the simulation, the necessary integration time step At has to be at least 1 fs or even less. This is a problem especially for calculations including explicit solvent molecules, because in the case of water they do not only increase the number of non-bonded interactions, they also increase the number of fast-moving hydrogen atoms. This particular situation is taken into account... 

Figure 7-14. All-atom and united-atom representation of the amino acid isoleucine. In this example, 13 atoms, which are able to form explicit non-bonding interactions, are reduced to only four pseudo-atoms,...

Ways to circumvent the above-mentioned problems have been to simply increase the cutoff distance to larger values, to use more than one cutoff value with different update frequencies, or to define more sophisticated cutoff schemes. In the last case, a truncation of the non-bonded interactions was replaced by shifting the interaction energies to zero or by additionally applying a switched sigmoidal func-... 

Although in teraetion s between vicinal I 4 atom s arc n om in ally treated as non bonded interactions, triost of the force fields treat these somewhat differently from normal 1 5 and greater non-bonded interactions. HyperCbern allows each of these 1 4 non-bonded interactions to be scaled down by a scale factor < 1.0 with AMBHR or OPI-S. bor HIO+ the electrostatic may be scaled and different param eters rn ay be ti sed for I 4 van dcr Waals interactions, fh e. AMBHR force field, for exam pie, n orrn a lly uses a seal in g factor of 0.5 for both van der Waals an d electrostatic interactions. 

OPTS (Optim i/.ed Potentials for Liquid Simulations) is based on a force field developed by the research group of Bill Jorgensen now at Yale University and previously at Purdue University. Like AMBER, the OPLS force field is designed for calculations on proteins an d nucleic acids. It in troduces non bonded in leraclion parameters that have been carefully developed from extensive Monte Carlo liquid sim u lation s of small molecules. These n on-bonded interactions have been added to the bonding interactions of AMBER to produce a new force field that is expected to be better than AMBER at describing simulations w here the solvent isexplic-... 

A typical force field model for propane contains ten bond-stretching terms, eighteen angle-bending terms, eighteen torsional terms and 27 non-bonded interactions. 

In a Urey-Bradley force field, angle bending is achieved using 1,3 non-bonded interaction rather than an explicit angle-bending potential. The stretch-bond term in such a forci field would be modelled by a harmonic function of the distance between the 1,3 atoms ... 

---