Bond angles field

With some appropriate definition of nearest neighbors (for example, the Voronoi construction, described in Section III.D, provides a unique definition of nearest neighbors), a bond angle field d can be defined, where 6 is the orientation of a bond between nearest neighbors with respect to some reference direction. In the continuum limit, the bond angle field is related to the displacement field by... 

Once the bond angle field is defined, a sixfold bond orientational order parameter field can be constructed,... 

Atomistically detailed models account for all atoms. The force field contains additive contributions specified in tenns of bond lengtlis, bond angles, torsional angles and possible crosstenns. It also includes non-bonded contributions as tire sum of van der Waals interactions, often described by Lennard-Jones potentials, and Coulomb interactions. Atomistic simulations are successfully used to predict tire transport properties of small molecules in glassy polymers, to calculate elastic moduli and to study plastic defonnation and local motion in quasi-static simulations [fy7, ( ]. The atomistic models are also useful to interiDret scattering data [fyl] and NMR measurements [70] in tenns of local order. 

But the methods have not really changed. The Verlet algorithm to solve Newton s equations, introduced by Verlet in 1967 [7], and it s variants are still the most popular algorithms today, possibly because they are time-reversible and symplectic, but surely because they are simple. The force field description was then, and still is, a combination of Lennard-Jones and Coulombic terms, with (mostly) harmonic bonds and periodic dihedrals. Modern extensions have added many more parameters but only modestly more reliability. The now almost universal use of constraints for bonds (and sometimes bond angles) was already introduced in 1977 [8]. That polarisability would be necessary was realized then [9], but it is still not routinely implemented today. Long-range interactions are still troublesome, but the methods that now become popular date back to Ewald in 1921 [10] and Hockney and Eastwood in 1981 [11]. 

Fig. 1. The time evolution (top) and average cumulative difference (bottom) associated with the central dihedral angle of butane r (defined by the four carbon atoms), for trajectories differing initially in 10 , 10 , and 10 Angstoms of the Cartesian coordinates from a reference trajectory. The leap-frog/Verlet scheme at the timestep At = 1 fs is used in all cases, with an all-atom model comprised of bond-stretch, bond-angle, dihedral-angle, van der Waals, and electrostatic components, a.s specified by the AMBER force field within the INSIGHT/Discover program.

The PEF is a sum of many individual contributions, Tt can be divided into bonded (bonds, angles, and torsions) and non-bonded (electrostatic and van der Waals) contributions V, responsible for intramolecular and, in tlic case of more than one molecule, also intermoleculai interactions. Figure 7-8 shows schematically these types of interactions between atoms, which arc included in almost all force field implementations. 

The Universal Force Field, UFF, is one of the so-called whole periodic table force fields. It was developed by A. Rappe, W Goddard III, and others. It is a set of simple functional forms and parameters used to model the structure, movement, and interaction of molecules containing any combination of elements in the periodic table. The parameters are defined empirically or by combining atomic parameters based on certain rules. Force constants and geometry parameters depend on hybridization considerations rather than individual values for every combination of atoms in a bond, angle, or dihedral. The equilibrium bond lengths were derived from a combination of atomic radii. The parameters [22, 23], including metal ions [24], were published in several papers. 

The potential energy of a molecular system in a force field is the sum of individnal components of the potential, such as bond, angle, and van der Waals potentials (equation H). The energies of the individual bonding components (bonds, angles, and dihedrals) are function s of th e deviation of a molecule from a h ypo-thetical compound that has bonded in teraction s at minimum val-n es. 

The force field ec uations for M.Vf+, AMBER, BlOg and OPES are similar in the types of terms they contain bond, angle, dihedral, van der Waals. and electrostatic. There are som e differences m the form s of the etinations that can al fect your ch oice of force field. 

Three-body and higher terms are sometimes incorporated into solid-state potentials. The Axilrod-Teller term is the most obvious way to achieve this. For systems such as the alkali halides this makes a small contribution to the total energy. Other approaches involve the use of terms equivalent to the harmonic angle-bending terms in valence force fields these have the advantage of simplicity but, as we have already discussed, are only really appropriate for small deviations from the equilibrium bond angle. Nevertheless, it can make a significant difference to the quality of the results in some cases. 

In general, we know bond lengths to within an uncertainty of 0.00.5 A — 0.5 pm. Bond angles are reliably known only to one or twx) degrees, and there arc many instances of more serious angle enxirs. Tn addition to experimental uncertainties and inaccuracies due to the model (lack of coincidence between model and molecule), some models present special problems unique to their geometry. For example, some force fields calculate the ammonia molecule. Nlln to be planar when there is abundant ex p er i m en ta I evidence th at N H is a 11 i g o n a I pyramid. 

Terms in the energy expression that describe how one motion of the molecule affects another are called cross terms. A cross term commonly used is a stretch-bend term, which describes how equilibrium bond lengths tend to shift as bond angles are changed. Some force fields have no cross terms and may compensate for this by having sophisticated electrostatic functions. The MM4 force field is at the opposite extreme with nine different types of cross terms. 

In addition to these basic terms, force fields often have cross terms that combine the above interactions. For example there may be a term which causes an angle bend to interact with a bond stretch term (opening a bond angle may tend to lengthen the bonds involved). 

A molecular dynamics force field is a convenient compilation of these data (see Chapter 2). The data may be used in a much simplified fonn (e.g., in the case of metric matrix distance geometry, all data are converted into lower and upper bounds on interatomic distances, which all have the same weight). Similar to the use of energy parameters in X-ray crystallography, the parameters need not reflect the dynamic behavior of the molecule. The force constants are chosen to avoid distortions of the molecule when experimental restraints are applied. Thus, the force constants on bond angle and planarity are a factor of 10-100 higher than in standard molecular dynamics force fields. Likewise, a detailed description of electrostatic and van der Waals interactions is not necessary and may not even be beneficial in calculating NMR strucmres. 

We refer to models where we write the total potential energy in terms of chemical endties such as bond lengths, bond angles, dihedral angles and so on as valence force field models. 

---