A stable, rapidly converging conjugate gradient method for energy minimization

We apply Shanno's conjugate gradient algorithm to the problem of minimizing the potential energy function associated with molecular mechanical calculations. Shanno's algorithm is stable with respect to roundoff errors and inexact line searches and converges rapidly to a minimum. Equally important, this algorithm can improve the rate of convergence to a minimum by a factor of 5 relative to Fletcher‐Reeves or Polak‐Ribière minimizers when used within the molecular mechanics package AMBER. Comparable improvements are found for a limited number of simulations when the Polak‐Ribière direction vector is incorporated into the Shanno algorithm.

[1]  Roger Fletcher,et al.  A Rapidly Convergent Descent Method for Minimization , 1963, Comput. J..

[2]  M. J. D. Powell,et al.  An efficient method for finding the minimum of a function of several variables without calculating derivatives , 1964, Comput. J..

[3]  C. M. Reeves,et al.  Function minimization by conjugate gradients , 1964, Comput. J..

[4]  M. J. D. Powell,et al.  Recent advances in unconstrained optimization , 1971, Math. Program..

[5]  Elijah Polak,et al.  Computational methods in optimization , 1971 .

[6]  H. Crowder,et al.  Linear convergence of the conjugate gradient method , 1972 .

[7]  P K Warme,et al.  Refinement of the x-ray structure of lysozyme by complete energy minimization. , 1974, Biochemistry.

[8]  Shmuel S. Oren,et al.  On the selection of parameters in Self Scaling Variable Metric Algorithms , 1974, Math. Program..

[9]  Shmuel S. Oren,et al.  Optimal conditioning of self-scaling variable Metric algorithms , 1976, Math. Program..

[10]  M. J. D. Powell,et al.  Some convergence properties of the conjugate gradient method , 1976, Math. Program..

[11]  David A. H. Jacobs,et al.  The State of the Art in Numerical Analysis. , 1978 .

[12]  M. J. D. Powell,et al.  Restart procedures for the conjugate gradient method , 1977, Math. Program..

[13]  D. Shanno On the Convergence of a New Conjugate Gradient Algorithm , 1978 .

[14]  David F. Shanno,et al.  Conjugate Gradient Methods with Inexact Searches , 1978, Math. Oper. Res..

[15]  Avinoam Perry,et al.  Technical Note - A Modified Conjugate Gradient Algorithm , 1978, Oper. Res..

[16]  Peter A. Kollman,et al.  AMBER: Assisted model building with energy refinement. A general program for modeling molecules and their interactions , 1981 .

[17]  M Karplus,et al.  Dynamics of proteins: elements and function. , 1983, Annual review of biochemistry.

[18]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[19]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[20]  E. Padlan,et al.  Refined crystal structure of deoxyhemoglobin S. I. Restrained least-squares refinement at 3.0-A resolution. , 1982, The Journal of biological chemistry.

[21]  M. Levitt,et al.  Protein normal-mode dynamics: trypsin inhibitor, crambin, ribonuclease and lysozyme. , 1985, Journal of molecular biology.