Minima hopping guided path search: an efficient method for finding complex chemical reaction pathways.

The Minima Hopping global optimization method uses physically realizable molecular dynamics moves in combination with an energy feedback that guarantees the escape from any potential energy funnel. For the purpose of finding reaction pathways, we argue that Minima Hopping is particularly suitable as a guide through the potential energy landscape and as a generator for pairs of minima that can be used as input structures for methods capable of finding transition states between two minima. For Lennard-Jones benchmark systems we compared this Minima Hopping guided path search method to a known approach for the exploration of potential energy landscapes that is based on deterministic mode-following. Although we used a stabilized mode-following technique that reliably allows to follow distinct directions when escaping from a local minimum, we observed that Minima Hopping guided path search is far superior in finding lowest-barrier reaction pathways. We, therefore, suggest that Minima Hopping guided path search can be used as a simple and efficient way to identify energetically low-lying chemical reaction pathways. Finally, we applied the Minima Hopping guided path search approach to 75-atom and 102-atom Lennard-Jones systems. For the 75-atom system we found pathways whose highest energies are significantly lower than the highest energy along the previously published lowest-barrier pathway. Furthermore, many of these pathways contain a smaller number of intermediate transition states than the previously publish lowest-barrier pathway. In case of the 102-atom system Minima Hopping guided path search found a previously unknown and energetically low-lying funnel.

[1]  Stefan Goedecker,et al.  Crystal structure prediction using the minima hopping method. , 2010, The Journal of chemical physics.

[2]  John H. Holland,et al.  Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence , 1992 .

[3]  C. Dellago,et al.  Precision shooting: Sampling long transition pathways. , 2008, The Journal of chemical physics.

[4]  A. Krall.  Applied Analysis , 1986 .

[5]  William H. Press,et al.  Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition , 1987 .

[6]  Edsger W. Dijkstra,et al.  A note on two problems in connexion with graphs , 1959, Numerische Mathematik.

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

[8]  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.

[9]  Hannes Jonsson,et al.  Reversible work transition state theory: application to dissociative adsorption of hydrogen , 1995 .

[10]  Robert T. Downs,et al.  Morphology-tuned wurtzite-type ZnS nanobelts , 2005, Nature materials.

[11]  Janet E. Jones,et al.  On the Calculation of Certain Crystal Potential Constants, and on the Cubic Crystal of Least Potential Energy , 1925 .

[12]  Linda Broadbelt,et al.  Computational tools for guided discovery and engineering of metabolic pathways. , 2013, Methods in molecular biology.

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

[14]  Stefan Goedecker,et al.  An enhanced splined saddle method. , 2011, The Journal of chemical physics.

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

[16]  Normand Mousseau,et al.  Evolution of the potential-energy surface of amorphous silicon. , 2010, Physical review letters.

[17]  David J. Wales,et al.  Energy landscapes: calculating pathways and rates , 2006 .

[18]  J. Doye,et al.  Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms , 1997, cond-mat/9803344.

[19]  N. Mousseau,et al.  The complex folding pathways of protein A suggest a multiple-funnelled energy landscape. , 2008, The Journal of chemical physics.

[20]  G. Henkelman,et al.  Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points , 2000 .

[21]  David J Wales,et al.  Folding pathways and rates for the three-stranded beta-sheet peptide Beta3s using discrete path sampling. , 2008, The journal of physical chemistry. B.

[22]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

[23]  D. Wales.  Discrete path sampling , 2002 .

[24]  Christoph Dellago,et al.  Efficient transition path sampling: Application to Lennard-Jones cluster rearrangements , 1998 .

[25]  K. Fukui.  Formulation of the reaction coordinate , 1970 .

[26]  David J Wales,et al.  Comparison of double-ended transition state search methods. , 2007, The Journal of chemical physics.

[27]  J. Doye,et al.  THE DOUBLE-FUNNEL ENERGY LANDSCAPE OF THE 38-ATOM LENNARD-JONES CLUSTER , 1998, cond-mat/9808265.

[28]  Thomas F. Miller,et al.  Sampling diffusive transition paths. , 2006, The Journal of chemical physics.

[29]  Jean-Francois Joly,et al.  The Activation-Relaxation Technique: ART Nouveau and Kinetic ART , 2012 .

[30]  G. Barkema,et al.  Traveling through potential energy landscapes of disordered materials: The activation-relaxation technique , 1997, cond-mat/9710023.

[31]  Mark A. Miller,et al.  Archetypal energy landscapes , 1998, Nature.

[32]  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.

[33]  A. Oganov,et al.  Crystal structure prediction using ab initio evolutionary techniques: principles and applications. , 2006, The Journal of chemical physics.

[34]  A. Châtelain,et al.  The European Physical Journal D , 1999 .

[35]  Normand Mousseau,et al.  Energy landscape of small clusters of self-interstitial dumbbells in iron , 2011 .

[36]  J. Pintér.  Global optimization : scientific and engineering case studies , 2006 .

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

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

[39]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[40]  Stefan Goedecker,et al.  Bell-Evans-Polanyi principle for molecular dynamics trajectories and its implications for global optimization. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[41]  Normand Mousseau,et al.  Understanding the EF-hand closing pathway using non-biased interatomic potentials. , 2012, The Journal of chemical physics.

[42]  Normand Mousseau,et al.  Numerical characterization of the Ga interstitial self-diffusion mechanisms in GaAs , 2008 .

[43]  J. Doye,et al.  Evolution of the Potential Energy Surface with Size for Lennard-Jones Clusters , 1999, cond-mat/9903305.

[44]  J. Gale,et al.  The prediction of inorganic crystal structures using a genetic algorithm and energy minimisation , 1999 .

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

[46]  Paul M Zimmerman,et al.  Efficient exploration of reaction paths via a freezing string method. , 2011, The Journal of chemical physics.

[47]  R. Rosenfeld.  Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[48]  Roi Baer,et al.  A spline for your saddle. , 2008, The Journal of chemical physics.

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

[50]  Mills,et al.  Quantum and thermal effects in H2 dissociative adsorption: Evaluation of free energy barriers in multidimensional quantum systems. , 1994, Physical review letters.

[51]  N. Mousseau,et al.  Dynamics of lennard-jones clusters: A characterization of the activation-relaxation technique , 2000 .

[52]  Normand Mousseau,et al.  Exploring the energy landscape of proteins: A characterization of the activation-relaxation technique , 2002 .

[53]  Stefan Goedecker,et al.  Structure of large gold clusters obtained by global optimization using the minima hopping method , 2009 .

[54]  C. Q. Lee,et al.  The Computer Journal , 1958, Nature.

[55]  Pathways for conformational change in nitrogen regulatory protein C from discrete path sampling. , 2008, The journal of physical chemistry. B.

[56]  N. Mousseau,et al.  Large loop conformation sampling using the activation relaxation technique, ART‐nouveau method , 2012, Proteins.

[57]  Lindsey J. Munro,et al.  DEFECT MIGRATION IN CRYSTALLINE SILICON , 1999 .

[58]  C. Dellago,et al.  Transition path sampling and the calculation of rate constants , 1998 .

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

[60]  D. Anderson,et al.  Algorithms for minimization without derivatives , 1974 .

[61]  V. Bazterra,et al.  Modified genetic algorithm to model crystal structures. I. Benzene, naphthalene and anthracene , 2002 .

[62]  S. Goedecker.  Global Optimization with the Minima Hopping Method , 2010 .

[63]  David J. Wales,et al.  Rearrangements of 55‐atom Lennard‐Jones and (C60)55 clusters , 1994 .

[64]  David J. Wales,et al.  Some further applications of discrete path sampling to cluster isomerization , 2004 .

[65]  N. Mousseau,et al.  Exploring energy landscapes of protein folding and aggregation. , 2008, Frontiers in bioscience : a journal and virtual library.

[66]  Thierry Deutsch,et al.  Optimized energy landscape exploration using the ab initio based activation-relaxation technique. , 2011, The Journal of chemical physics.

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

[68]  Christoph Dellago,et al.  Nucleation and growth in structural transformations of nanocrystals. , 2009, Nano letters.

[69]  Artem R. Oganov,et al.  Modern methods of crystal structure prediction , 2011 .

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

[71]  C. Dellago,et al.  Transition Path Sampling , 2005 .

[72]  Christoph Dellago,et al.  UvA-DARE ( Digital Academic Repository ) Role of the prestructured surface cloud in crystal nucleation , 2011 .

[73]  Stefan Goedecker,et al.  The performance of minima hopping and evolutionary algorithms for cluster structure prediction. , 2009, The Journal of chemical physics.

[74]  U. Landman,et al.  Genetic Algorithms for Structural Cluster Optimization , 1998 .

[75]  Janet E. Jones.  On the determination of molecular fields. —II. From the equation of state of a gas , 1924 .

[76]  Barkema,et al.  Event-Based Relaxation of Continuous Disordered Systems. , 1996, Physical review letters.

[77]  S. Goedecker,et al.  Questioning the existence of a unique ground-state structure for Si clusters , 2007 .

[78]  N. Mousseau,et al.  Energy landscapes of the monomer and dimer of the Alzheimer's peptide Abeta(1-28). , 2008, The Journal of chemical physics.

[79]  N. Mousseau,et al.  Ab initio study of the diffusion mechanisms of gallium in a silicon matrix , 2007 .

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

[81]  D. Wales,et al.  Theoretical study of the water pentamer , 1996 .

[82]  Peter Steenkiste,et al.  On path selection for traffic with bandwidth guarantees , 1997, Proceedings 1997 International Conference on Network Protocols.

[83]  Nancy R. Sottos,et al.  Biasing reaction pathways with mechanical force , 2007, Nature.

[84]  M. Karplus,et al.  The topology of multidimensional potential energy surfaces: Theory and application to peptide structure and kinetics , 1997 .

[85]  David Chandler,et al.  Transition path sampling: throwing ropes over rough mountain passes, in the dark. , 2002, Annual review of physical chemistry.