A Valley Following Method

We present a procedure to follow the "path along the valley floor" of a hypersurface. The aim is either to find minima, or to go from a minimum to a saddle point of index one, if the saddle is at the top of the valley floor. The motivation is that of taking into account local nonconvexity of the hypersurface and possibly to determine valleys. The method uses a projector technique where the projector is built by the tangent of the valley floor line. The projector is applied to the gradient and Hessian matrix of a given function, and it is used for predictor and corrector steps in path following. The resulting path is the "valley floor gradient extremal" which corresponds to the smallest (absolute) eigenvalue of the Hessian. Convergence properties are analysed.

[1]  Michael Hirsch,et al.  Improved RGF method to find saddle points , 2002, J. Comput. Chem..

[2]  H. H. Rosenbrock,et al.  An Automatic Method for Finding the Greatest or Least Value of a Function , 1960, Comput. J..

[3]  Pedro J. Zufiria,et al.  On the role of singularities in Branin's method from dynamic and continuation perspectives , 2002, Appl. Math. Comput..

[4]  Eugene L. Allgower,et al.  Numerical continuation methods - an introduction , 1990, Springer series in computational mathematics.

[5]  P. Deuflhard,et al.  Affine Invariant Convergence Theorems for Newton’s Method and Extensions to Related Methods , 1979 .

[6]  Klaus Ruedenberg,et al.  Gradient fields of potential energy surfaces , 1994 .

[7]  Wolfgang Quapp.  Comment On the quadratic reaction path evaluated in a reduced potential energy surface model and the problem to locate transition states [by J. M. Anglada, E. Besalú, J. M. Bofill, and R. Crehuet, J Comput Chem 2001, 22, 4, 387-406] , 2001, J. Comput. Chem..

[8]  A note on Branin's method for binding the critical points of smooth functions , 1987 .

[9]  Jun-qiang Sun,et al.  Gradient extremals and steepest descent lines on potential energy surfaces , 1993 .

[10]  Wolfgang Quapp,et al.  Searching minima of an N-dimensional surface: A robust valley following method , 2001 .

[11]  A. G. Shamov,et al.  The local definition of the Optimum ascent path on a multi-dimensional potential energy surface and its practical application for the location of saddle points , 1981 .

[12]  P. Landsberg.  On the Theory of Reaction Rates , 1953 .

[13]  Gradient Extremals and Their Relation to the Minimum Energy Path , 1995 .

[14]  Michael Hirsch,et al.  Bifurcation of reaction pathways: the set of valley ridge inflection points of a simple three-dimensional potential energy surface , 1998 .

[15]  C. Floudas,et al.  Primal-relaxed dual global optimization approach , 1993 .

[16]  Wolfgang Quapp,et al.  Gradient extremals and valley floor bifurcations on potential energy surfaces , 1989 .

[17]  David J Rowe,et al.  Valleys and fall lines on a Riemannian manifold , 1982 .

[18]  F. Jensen.  Locating transition structures by mode following: A comparison of six methods on the Ar8 Lennard‐Jones potential , 1995 .

[19]  Dietmar Heidrich,et al.  Analysis of the concept of minimum energy path on the potential energy surface of chemically reacting systems , 1984 .

[20]  Emili Besalú,et al.  On the quadratic reaction path evaluated in a reduced potential energy surface model and the problem to locate transition states * , 2001 .

[21]  H. E. Avery.  Theory of Reaction Rates , 1974 .

[22]  Ahmed H. Zewail,et al.  Femtochemistry: Atomic-Scale Dynamics of the Chemical Bond† , 2000 .

[23]  Robert Schaback,et al.  An extended continuous Newton method , 1990 .

[24]  Michael Hirsch,et al.  Following the streambed reaction on potential-energy surfaces: a new robust method , 2000 .

[25]  Michael Hirsch,et al.  Searching for saddle points of potential energy surfaces by following a reduced gradient , 1998, J. Comput. Chem..

[26]  Hubertus Th. Jongen,et al.  Nonlinear Optimization in Finite Dimensions - Morse Theory, Chebyshev Approximation, Transversality, Flows, Parametric Aspects , 2000 .

[27]  David K. Hoffman,et al.  Gradient extremals , 1986 .

[28]  Farhad Arbab,et al.  Distributed Evolutionary Optimization, in Manifold: Rosenbrock's Function Case Study , 2000, Inf. Sci..

[29]  F. H. Branin.  Widely convergent method for finding multiple solutions of simultaneous nonlinear equations , 1972 .