Stabilized quasi-Newton optimization of noisy potential energy surfaces.

Optimizations of atomic positions belong to the most commonly performed tasks in electronic structure calculations. Many simulations like global minimum searches or characterizations of chemical reactions require performing hundreds or thousands of minimizations or saddle computations. To automatize these tasks, optimization algorithms must not only be efficient but also very reliable. Unfortunately, computational noise in forces and energies is inherent to electronic structure codes. This computational noise poses a severe problem to the stability of efficient optimization methods like the limited-memory Broyden-Fletcher-Goldfarb-Shanno algorithm. We here present a technique that allows obtaining significant curvature information of noisy potential energy surfaces. We use this technique to construct both, a stabilized quasi-Newton minimization method and a stabilized quasi-Newton saddle finding approach. We demonstrate with the help of benchmarks that both the minimizer and the saddle finding approach are superior to comparable existing methods.

[1]  Stefan Goedecker,et al.  Global minimum determination of the Born-Oppenheimer surface within density functional theory. , 2005, Physical review letters.

[2]  R. Fletcher,et al.  A New Approach to Variable Metric Algorithms , 1970, Comput. J..

[3]  W. Miller,et al.  ON FINDING TRANSITION STATES , 1981 .

[4]  Siegfried Schmauder,et al.  Comput. Mater. Sci. , 1998 .

[5]  Pekka Koskinen,et al.  Structural relaxation made simple. , 2006, Physical review letters.

[6]  H. Eyring.  The Activated Complex in Chemical Reactions , 1935 .

[7]  S. Goedecker.  Minima hopping: an efficient search method for the global minimum of the potential energy surface of complex molecular systems. , 2004, The Journal of chemical physics.

[8]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[9]  F. Jensen.  Introduction to Computational Chemistry , 1998 .

[10]  G. Kresse,et al.  Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set , 1996 .

[11]  D. Shanno.  Conditioning of Quasi-Newton Methods for Function Minimization , 1970 .

[12]  J. Nocedal.  Updating Quasi-Newton Matrices With Limited Storage , 1980 .

[13]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[14]  G. Henkelman,et al.  Comparison of methods for finding saddle points without knowledge of the final states. , 2004, The Journal of chemical physics.

[15]  Arthur F. Voter,et al.  Highly optimized empirical potential model of silicon , 2000 .

[16]  Stefan Goedecker,et al.  Linear scaling relaxation of the atomic positions in nanostructures , 2001 .

[17]  K. Varga,et al.  Stochastic Variational Approach to Quantum-Mechanical Few-Body Problems , 1998 .

[18]  Paul Sherwood,et al.  Superlinearly converging dimer method for transition state search. , 2008, The Journal of chemical physics.

[19]  K. Tamura,et al.  Metabolic engineering of plant alkaloid biosynthesis. Proc Natl Acad Sci U S A , 2001 .

[20]  D. Wales.  Locating stationary points for clusters in cartesian coordinates , 1993 .

[21]  Car,et al.  Acceleration schemes for ab initio molecular-dynamics simulations and electronic-structure calculations. , 1994, Physical review. B, Condensed matter.

[22]  Reinhold Schneider,et al.  Daubechies wavelets as a basis set for density functional pseudopotential calculations. , 2008, The Journal of chemical physics.

[23]  P. Pulay.  Improved SCF convergence acceleration , 1982 .

[24]  Giovanni Ciccotti,et al.  Book Review: Classical and Quantum Dynamics in Condensed Phase Simulations , 1998 .

[25]  Julius Jellinek,et al.  Energy Landscapes: With Applications to Clusters, Biomolecules and Glasses , 2005 .

[26]  Hafner,et al.  Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. , 1994, Physical review. B, Condensed matter.

[27]  D. Goldfarb.  A family of variable-metric methods derived by variational means , 1970 .

[28]  S. Goedecker,et al.  Metrics for measuring distances in configuration spaces. , 2013, The Journal of chemical physics.

[29]  Stefan Goedecker,et al.  Minima hopping guided path search: an efficient method for finding complex chemical reaction pathways. , 2014, The Journal of chemical physics.

[30]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[31]  G. Kresse,et al.  From ultrasoft pseudopotentials to the projector augmented-wave method , 1999 .

[32]  Stefan Goedecker,et al.  Daubechies wavelets for linear scaling density functional theory. , 2014, The Journal of chemical physics.

[33]  G. Henkelman,et al.  A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives , 1999 .

[34]  I. Mayer.  Simple theorems, proofs, and derivations in quantum chemistry , 2003 .

[35]  Graeme Henkelman,et al.  EON: software for long time simulations of atomic scale systems , 2014 .

[36]  D H Weinstein,et al.  Modified Ritz Method. , 1934, Proceedings of the National Academy of Sciences of the United States of America.

[37]  A. Bell,et al.  Efficient methods for finding transition states in chemical reactions: comparison of improved dimer method and partitioned rational function optimization method. , 2005, The Journal of chemical physics.

[38]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[39]  Stefan Goedecker.  Optimization and parallelization of a force field for silicon using OpenMP ? ? This program can be d , 2002 .

[40]  Michael Page,et al.  On evaluating the reaction path Hamiltonian , 1988 .

[41]  C. G. Broyden.  The Convergence of a Class of Double-rank Minimization Algorithms 1. General Considerations , 1970 .

[42]  M. Tuckerman,et al.  IN CLASSICAL AND QUANTUM DYNAMICS IN CONDENSED PHASE SIMULATIONS , 1998 .

[43]  Matt Probert,et al.  Improved algorithm for geometry optimisation using damped molecular dynamics , 2003 .

[44]  G. Kresse,et al.  Ab initio molecular dynamics for liquid metals. , 1993 .