Computation bond energy terms

Finally, there are a number of entirely mundane (but still very worthwhile ) steps that can be taken to reduce the total computer time required for a MD simulation. As a single example, note that any force on a particle derived from a force-field non-bonded energy term is induced by some other particle (i.e., the potential is pairwise). Newton s Third Law tells us that... 

AH°t) and the sum of bond-energy terms ( °b) of the species under consideration. Here, Nab stands for the number of equivalent bonds having E°b. The bond-energy terms are obtained from a set of reference compounds. The heat of atomization is evaluated from the difference of the heat of formation of the constituent atoms and the species under consideration (2). The procedure depends on the reliability of the heat of formation of the individual molecules. In cases where experimental values exist it is not necessary to have recourse to computational methods. However, the heats of formation of many molecules which are of interest in the context of captodative substitution are not known experimentally. 

Once the bond orders have been calculated, they are used to compute bond energies, angles, and torsions. These terms are also used in non-reactive force fields, but their use in reactive simulations requires some modifications. All covalent terms are pre-multiplied by the bond orders involved this ensures that whenever a bond is broken, all the terms involving it vanish smoothlv. Also, the equilibrium angle in covalent-angle terms depends upon the bond orders... 

In the decade following Coulson s paper, computers became sufficiently fast and numerous to allow Hartree-Fock calculations on small hydrogen bonded systems. In 1971, Morokuma defined wavefunctions associated with Coulson s energy terms (except dispersion, which non-correlated wavefunctions can not describe) and computed the energy terms. A brief description follows in which and E are, respectively, the Hartree-Fock wavefunction and energy calculated for the isolated molecule A (and similarly for B). 

Much like the RISM method, the LD approach is intermediate between a continuum model and an explicit model. In the limit of an infinite dipole density, the uniform continuum model is recovered, but with a density equivalent to, say, the density of water molecules in liquid water, some character of the explicit solvent is present as well, since the magnitude of the dipoles and their polarizability are chosen to mimic the particular solvent (Papazyan and Warshel 1997). Since the QM/MM interaction in this case is purely electrostatic, other non-bonded interaction terms must be included in order to compute, say, solvation free energies. When the same surface-tension approach as that used in many continuum models is adopted (Section 11.3.2), the resulting solvation free energies are as accurate as those from pure continuum models (Florian and Warshel 1997). Unlike atomistic models, however, the use of a fixed grid does not permit any real information about solvent structure to be obtained, and indeed the fixed grid introduces issues of how best to place the solute into the grid, where to draw the solute boundary, etc. These latter limitations have curtailed the application of the LD model. 

The first term in Eq. (8) is due to the difference in entropy between surface layer and bulk. The decrease in entropy on enrichment is balanced by the gain in enthalpy, determined by the difference between the numbers of bonds broken in the surface and bulk on enrichment. e2 and <=, are the bond energies and can be computed from the heat of sublimation. In an extensive review, Overbury et al. (12) have shown that the surface energy of a metal is proportional to the heat of sublimation of the metals. 

---