Coulomb interaction/integral force fields

Chymotrypsin, 170,171, 172, 173 Classical partition functions, 42,44,77 Classical trajectories, 78, 81 Cobalt, as cofactor for carboxypeptidase A, 204-205. See also Enzyme cofactors Condensed-phase reactions, 42-46, 215 Configuration interaction treatment, 14,30 Conformational analysis, 111-117,209 Conjugated gradient methods, 115-116. See also Energy minimization methods Consistent force field approach, 113 Coulomb integrals, 16, 27 Coulomb interactions, in macromolecules, 109, 123-126... 

The GEM force field follows exactly the SIBFA energy scheme. However, once computed, the auxiliary coefficients can be directly used to compute integrals. That way, the evaluation of the electrostatic interaction can virtually be exact for an perfect fit of the density as the three terms of the coulomb energy, namely the nucleus-nucleus repulsion, electron-nucleus attraction and electron-electron repulsion, through the use of p [2, 14-16, 58],... 

The philosophical transition from the atomic prejudice to a view of intermolecular interaction in terms of diffuse electron density has its proper computational counterpart in full quantum mechanical calculations, which, however, cannot at present provide complete intermolecular energies because of limitations in the treatment of electron correlation, a major ingredient of the intermolecular interaction recipe. In a different perspective, the classical atom-atom force-field approach is widely applicable but entirely parametric and of scarce adherence to physical principles. The need is felt for an extension to represent in a more realistic manner the effects of diffuse electron clouds. This is done in the so-called semi-classical density sums (SCDS) or briefly. Pixel approach [9], which will now be described. The Pixel method is based on numerical integrations over molecular electron densities, and allows a separation of the total intermolecular cohesion energy into coulombic, polarization, dispersion, and repulsion contributions. 

Here r and v are respectively the electron position and velocity, r = —(e2 /em)(r/r3) is the acceleration in the coulombic field of the positive ion and q = /3kBT/m. The mobility of the quasi-free electron is related to / and the relaxation time T by p = e/m/3 = et/m, so that fi = T l. In the spherically symmetrical situation, a density function n(vr, vt, t) may be defined such that n dr dvr dvt = W dr dv here, vr and vt and are respectively the radical and normal velocities. Expectation values of all dynamical variables are obtained from integration over n. Since the electron experiences only radical force (other than random interactions), it is reasonable to expect that its motion in the v space is basically a free Brownian motion only weakly coupled to r and vr by the centrifugal force. The correlations1, K(r, v,2) and fc(vr, v(2) are then neglected. Another condition, cr(r)2 (r)2, implying that the electron distribution is not too much delocalized on r, is verified a posteriori. Following Chandrasekhar (1943), the density function may now be written as an uncoupled product, n = gh, where... 

The Fock operator determines three sets of information for each electron i (1) the kinetic energy term of the electron (—1/2V ), (2) an attraction term with each nucleus, A, (—EZA/r,A), and (3) the interaction of the electron with all the other electrons in the molecule. This average force is treated by the (IJjj — Kjj) term and can be described as the potential felt by a single electron in the field of the other i — 1 electrons in the molecule. A few words about the components of this last term in the Fock operator are in order. J is called the coulomb operator and is identified as the classical repulsion between electrons. The exchange integral K is due to the quantum mechanical effect of spin correlation, an intrinsic property of the electron that keeps apart electrons of the same spin. This operator has a stabilizing effect on the energy of the system. 

---