Field Forces

Symmetry Force Constants of Phosphane. For explanation, see the text. 

Scaling of the angle-bending coordinates with a unit bond length of 1A yielded constants in mdyn/A. 

Several theoretically and empirically derived and estimated force constants of the quadratic valence force field were reported [41 to 55, 58, 63]. 

The 14 independent parameters of the cubic force field of PH3 were determined by ab initio SCF calculations [57]. 

Although quadratic expressions are reasonable approximations, as discussed above, they only works for small deviations from the natural state 

Some additional comments regarding Equation 4 are in order. The factor 143.88 converts the units to kcal/mol. There are two additional constants. The first is ks, which is the stretching force constant parameter in units of md A 1. The second constant is cs, which is the cubic term with a unitless value of 2.55. When a Morse potential is expanded in a power series, the factor 7/12 is obtained. 

In the MM3 program, the bending motions are divided into two separate descriptions in-plane and out-of-plane bending. The in-plane bending is shown in Equation 6. Depending on the geometry of the central atom in the 

The discussion thus far has focused on the forces between an array of atoms connected together through covalent bonds and their angles. Important interactions occur between atoms not directly bonded together. The theoretical explanation for attractive and repulsive forces for nonbonded atoms i and j is based on electron distributions. The motion of electrons about a nucleus creates instantaneous dipoles. The instantaneous dipoles on atom i induce dipoles of opposite polarity on atom j. The interactions between the instantaneous dipole on atom i with the induced instantaneous dipole on atom j of the two electron clouds of nonbonded atoms are responsible for attractive interactions. The attractive interactions are know as London Dispersion forces,70 which are related to r 6, where r is the distance between nonbonded atoms i and j. As the two electron clouds of nonbonded atoms i and j approach one another, they start to overlap. There is a point where electron-electron and nuclear-nuclear repulsion of like charges overwhelms the London Dispersion forces.33 The repulsive 

In MM3, the two types of nonbonded interactions which are included in the program are handled separately. The first includes attractive and repulsive forces between nonbonded atoms, their origins are described above, and the second is hydrogen bonding. 

---

The idea that unsymmetrical molecules will orient at an interface is now so well accepted that it hardly needs to be argued, but it is of interest to outline some of the history of the concept. Hardy [74] and Harkins [75] devoted a good deal of attention to the idea of force fields around molecules, more or less intense depending on the polarity and specific details of the structure. Orientation was treated in terms of a principle of least abrupt change in force fields, that is, that molecules should be oriented at an interface so as to provide the most gradual transition from one phase to the other. If we read interaction energy instead of force field, the principle could be reworded on the very reasonable basis that molecules will be oriented so that their mutual interaction energy will be a maximum. 

Hardy s explanation that the small coefficients of friction observed under boundary lubrication conditions were due to the reduction in the force fields between the surfaces as a result of adsorbed films is undoubtedly correct in a general way. The explanation leaves much to be desired, however, and it is of interest to consider more detailed proposals as to the mechanism of boundary lubrication. 

Kramers solution of the barrier crossing problem [45] is discussed at length in chapter A3.8 dealing with condensed-phase reaction dynamics. As the starting point to derive its simplest version one may use the Langevin equation, a stochastic differential equation for the time evolution of a slow variable, the reaction coordinate r, subject to a rapidly statistically fluctuating force F caused by microscopic solute-solvent interactions under the influence of an external force field generated by the PES F for the reaction... 

For a very large number of variables, the question of storing the approximate Hessian or inverse Hessian F becomes important. Wavefunction optimization problems can have a very large number of variables, a million or more. Geometry optimization at the force field level can also have thousands of degrees of freedom. In these cases, the initial inverse Hessian is always taken to be diagonal or sparse, and it is best to store the... 

Pulay P and Meyer W 1971 Ab initio calculation of the force field of ethylene J. Moi. Spectrosc. 40 59... 

The complexity of polymeric systems make tire development of an analytical model to predict tlieir stmctural and dynamical properties difficult. Therefore, numerical computer simulations of polymers are widely used to bridge tire gap between tire tlieoretical concepts and the experimental results. Computer simulations can also help tire prediction of material properties and provide detailed insights into tire behaviour of polymer systems. A simulation is based on two elements a more or less detailed model of tire polymer and a related force field which allows tire calculation of tire energy and tire motion of tire system using molecular mechanisms, molecular dynamics, or Monte Carlo teclmiques 1631. 

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. 

The present high cost of full CASSCF direct dynamics means that it is not possible to use such calculations to run large numbers of trajectories. As a result it cannot be used to build up experience of the types of effects to be found in dynamical studies of organic photochemistry, and in their interpretation. This problem can be remedied by performing calculations using the MMVB force field [63,64]. 

MMVB is a hybrid force field, which uses MM to treat the unreactive molecular framework, combined with a valence bond (VB) approach to treat the reactive part. The MM part uses the MM2 force field [58], which is well adapted for organic molecules. The VB part uses a parametrized Heisenberg spin Hamiltonian, which can be illustrated by considering a two orbital, two electron description of a sigma bond described by the VB determinants... 

The MMVB force field has also been used with Ehrenfest dynamics to propagate trajectories using mixed-state forces [84]. The motivation for this is... 

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]. 

With the danger of severe oversimplification, which unavoidably leads to improper underevaluation of important recent developments, I shall try to indicate where traditional, classical MD has brought us today, or will bring us tomorrow. This concerns the techniques rather than the applications, which cannot be reviewed in the present context. The main aspects to consider concern algorithms and force fields. 

These numbers are not negligible. At present such effects are on the average compensated by other force field terms through empirical parametrization. 

While simulations reach into larger time spans, the inaccuracies of force fields become more apparent on the one hand properties based on free energies, which were never used for parametrization, are computed more accurately and discrepancies show up on the other hand longer simulations, particularly of proteins, show more subtle discrepancies that only appear after nanoseconds. Thus force fields are under constant revision as far as their parameters are concerned, and this process will continue. Unfortunately the form of the potentials is hardly considered and the refinement leads to an increasing number of distinct atom types with a proliferating number of parameters and a severe detoriation of transferability. The increased use of quantum mechanics to derive potentials will not really improve this situation ab initio quantum mechanics is not reliable enough on the level of kT, and on-the-fly use of quantum methods to derive forces, as in the Car-Parrinello method, is not likely to be applicable to very large systems in the foreseeable future. 

This situation, despite the fact that reliability is increasing, is very undesirable. A considerable effort will be needed to revise the shape of the potential functions such that transferability is greatly enhanced and the number of atom types can be reduced. After all, there is only one type of carbon it has mass 12 and charge 6 and that is all that matters. What is obviously most needed is to incorporate essential many-body interactions in a proper way. In all present non-polarisable force fields many-body interactions are incorporated in an average way into pair-additive terms. In general, errors in one term are compensated by parameter adjustments in other terms, and the resulting force field is only valid for a limited range of environments. 

However, in many applications the essential space cannot be reduced to only one degree of freedom, and the statistics of the force fluctuation or of the spatial distribution may appear to be too poor to allow for an accurate determination of a multidimensional potential of mean force. An example is the potential of mean force between two ions in aqueous solution the momentaneous forces are two orders of magnitude larger than their average which means that an error of 1% in the average requires a simulation length of 10 times the correlation time of the fluctuating force. This is in practice prohibitive. The errors do not result from incorrect force fields, but they are of a statistical nature even an exact force field would not suffice. 

Aqvist, J., Warshel, A. Simulation of enzyme reactions using valence bond force fields and other hybrid quantum/classical approaches. Chem. Rev. 93... 

A number of issues need to be addressed before this method will become a routine tool applicable to problems as the conformational equilibrium of protein kinase. E.g. the accuracy of the force field, especially the combination of Poisson-Boltzmann forces and molecular mechanics force field, remains to be assessed. The energy surface for the opening of the two kinase domains in Pig. 2 indicates that intramolecular noncovalent energies are overestimated compared to the interaction with solvent. 

Conformational Adjustments The conformations of protein and ligand in the free state may differ from those in the complex. The conformation in the complex may be different from the most stable conformation in solution, and/or a broader range of conformations may be sampled in solution than in the complex. In the former case, the required adjustment raises the energy, in the latter it lowers the entropy in either case this effect favors the dissociated state (although exceptional instances in which the flexibility increases as a result of complex formation seem possible). With current models based on two-body potentials (but not with force fields based on polarizable atoms, currently under development), separate intra-molecular energies of protein and ligand in the complex are, in fact, definable. However, it is impossible to assign separate entropies to the two parts of the complex. 

The problems that occur when one tries to estimate affinity in terms of component terms do not arise when perturbation methods are used with simulations in order to compute potentials of mean force or free energies for molecular transformations simulations use a simple physical force field and thereby implicitly include all component terms discussed earlier. We have used the molecular transformation approach to compute binding affinities from these first principles [14]. The basic approach had been introduced in early work, in which we studied the affinity of xenon for myoglobin [11]. The procedure was to gradually decrease the interactions between xenon atom and protein, and compute the free energy change by standard perturbation methods, cf. (10). An (issential component is to impose a restraint on the... 

In our implementation of SMD, modified versions of VMD and Sigma communicate with each other using a customized, lightweight protocol. Sigma sends atomic positions resulting from each molecular dynamics time step to VMD for display. When the user specifies restraints on parts of the displayed model, VMD sends them to Sigma, where they are converted into potential-well restraints added to the force field [21]. 

---