Molecular geometry algorithms

The calculation of reaction rates has not seen as the widespread use as the calculation of molecular geometries. In recent years, it has become possible to compute reaction rates with reasonable accuracy. However, these calculations require some expertise on the part of the researcher. This is partly because of the difficulty in obtaining transition structures and partly because reaction rate algorithms have not been integrated into major computational chemistry programs and thus become automated. 

The basis for the determination of solution conformation from NMR data lies in the determination of cross relaxation rates between pairs of protons from cross peak intensities in two-dimensional nuclear Overhauser effect (NOE) experiments. In the event that pairs of protons may be assumed to be rigidly fixed in an isotopically tumbling sphere, a simple inverse sixth power relationship between interproton distances and cross relaxation rates permits the accurate determination of distances. Determination of a sufficient number of interproton distance constraints can lead to the unambiguous determination of solution conformation, as illustrated in the early work of Kuntz, et al. (25). While distance geometry algorithms remain the basis of much structural work done today (1-4), other approaches exist. For instance, those we intend to apply here represent NMR constraints as pseudoenergies for use in molecular dynamics or molecular mechanics programs (5-9). 

The example of the prototype SN2 reaction demonstrates that the algorithm of MD along the IRP applied here works correctly also in the case of a reaction with relatively large deviations of the molecular geometries at finite temperatures from the corresponding structures at the zero-temperature IRP. 

These methods can be combined with geometry optimization as well as with molecular dynamics algorithms, with forces obtained from the gradients of the total quantum energy [10]. This equally applies to all quantum methods, quoted in the following. 

This numerical problem of integration can be avoided using the ADMA technique. Within the ADMA method, the integration in Eq. (361) can be performed using the analytical expressions of macromolecular density matrices and AOs. As an option of the ADMA algorithm, the calculated ADMA Hellmann-Feynman forces can be used for macro-molecular geometry optimization and macromolecular conformational analysis. 

To optimise the geometry, the energy must be expressed as a function of atomic displacements. This yields the partial derivatives crucial to automatic minimisation algorithms. The expressions for the total energy derivatives with respect to atomic displacements are quite complex for ab initio and semi-empirical methods but trivial for empirical schemes like Molecular Mechanics (MM). Virtually all modern computer codes provide extensive, efficient facilities for determining ground state molecular geometries. 

---