SIMULATED ANNEALING TO LOCATE VARIOUS STATIONARY-POINTS IN SEMIEMPIRICAL METHODS

In this work, an algorithm was developed to study the potential energy surfaces in the coordinate spaces of molecules by a nonlocal way, in contrast to classic energy minimizers as the BFGS or the DFP method. This algorithm, based on the specificities of semiempirical methods, mixes simulated annealing and local searches to reduce computation costs. By this technique, the global energy minimum can be localized. Moreover, local minima that are close in energy to the global minimum are also obtained. If the search is not only for minima but for all stationary points (minima, saddle points...), then the energy is replaced by the gradient norm, which reaches its minimum values at stationary points. The annealing process is stopped before having accurately reached the global minimum and generates a list of geometries whose energies (respectively, whose gradients) are optimized by local minimizers. This list of geometries is shortened from the nearly equivalent geometries by a dynamic single-clustering analysis. The energy/gradient local minimizers act on the clustered list to produce a set of minima/stationary points. A targeted search of these points and reduction of the costs are reached by the way of several penalty functions. They eliminate - without energy calculation-most of the points generated by random walks on the potential energy surface. These penalty functions (on the total moment of inertia or on interatomic distances) are specific to the class of problem studied. They account for the nonrupture of either specific chemical bonds or rings in cyclic molecules, they assure that molecular systems are kept bonded, and they avoid the collapsing of atoms.