Efficient Basin-Hopping Sampling of Reaction Intermediates through Molecular Fragmentation and Graph Theory.

Basin-hopping sampling has been widely used for searching local minima on a potential energy surface. Reaction intermediates including reactants and products are also local minima composed of a reaction path, but their brute-force sampling is too demanding because of large degrees of freedom. We developed an efficient Monte Carlo basin-hopping method to sample reaction intermediates through the fragmentation of molecules and a postanalysis scheme using the graph theory with a matrix representation of molecular structures. The former greatly reduces the dimension of a given potential energy surface, while the latter offers not only the effective screening of resulting local minima toward desirable intermediates but also their automatic ordering along a reaction path. We combined it with the density functional tight binding method for rapid calculations and tested its performance for organic reactions.

[1]  L. Pogliani.  From molecular connectivity indices to semiempirical connectivity terms: recent trends in graph theoretical descriptors. , 2000, Chemical reviews.

[2]  Sándor Suhai,et al.  Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties , 1998 .

[3]  W. R. Wadt,et al.  Ab initio effective core potentials for molecular calculations , 1984 .

[4]  A. Schug,et al.  Basin hopping simulations for all-atom protein folding. , 2006, The Journal of chemical physics.

[5]  Yanhui Shi,et al.  A Titanium-Catalyzed Three-Component Coupling To Generate α,β-Unsaturated β-Iminoamines , 2003 .

[6]  S. Bulusu,et al.  A global search of highly stable gold-covered bimetallic clusters M@Au(n) (n=8-17): endohedral gold clusters. , 2006, Chemphyschem : a European journal of chemical physics and physical chemistry.

[7]  Yanhui Shi,et al.  Titanium-catalyzed iminohydrazination of alkynes , 2005 .

[8]  A. Odom.  New C-N and C-C bond forming reactions catalyzed by titanium complexes. , 2005, Dalton transactions.

[9]  Craig A. Shelley,et al.  Heuristic approach for displaying chemical structures , 1983, J. Chem. Inf. Comput. Sci..

[10]  J. Moussa.  Comment on "Fast and accurate modeling of molecular atomization energies with machine learning". , 2012, Physical review letters.

[11]  Wing-Ki Liu,et al.  Multicanonical basin hopping: a new global optimization method for complex systems. , 2004, The Journal of chemical physics.

[12]  K. Morokuma,et al.  Finding Reaction Pathways of Type A + B → X: Toward Systematic Prediction of Reaction Mechanisms. , 2011, Journal of chemical theory and computation.

[13]  Klaus-Robert Müller,et al.  Assessment and Validation of Machine Learning Methods for Predicting Molecular Atomization Energies. , 2013, Journal of chemical theory and computation.

[14]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[15]  John P. Wolfe,et al.  Rational Development of Practical Catalysts for Aromatic Carbon−Nitrogen Bond Formation , 1998 .

[16]  K. Lam,et al.  DFT Studies on the Effect of the Nature of the Aryl Halide Y−C6H4−X on the Mechanism of Its Oxidative Addition to Pd0L versus Pd0L2 , 2007 .

[17]  Jijun Zhao,et al.  B80 and B101-103 clusters: remarkable stability of the core-shell structures established by validated density functionals. , 2012, The Journal of chemical physics.

[18]  Richard F. Heck,et al.  The Reaction of Cobalt Hydrotetracarbonyl with Olefins , 1961 .

[19]  J. Doye,et al.  Thermodynamics and the Global Optimization of Lennard-Jones clusters , 1998, cond-mat/9806020.

[20]  S. Woodley,et al.  Crystal structure prediction from first principles. , 2008, Nature materials.

[21]  Danail Bonchev,et al.  Application of Graph Theory to Chemical Kinetics, 3. Topological Specificity of Multiroute Reaction Mechanisms , 1995, J. Chem. Inf. Comput. Sci..

[22]  Benjamin D. Ward,et al.  Insertions into azatitanacyclobutenes: new insights into three-component coupling reactions involving imidotitanium intermediates , 2008 .

[23]  Beatriz Cordero,et al.  Covalent radii revisited. , 2008, Dalton transactions.

[24]  P J Hsu,et al.  Structures of bimetallic clusters. , 2006, The Journal of chemical physics.

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

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

[27]  Zhi-Pan Liu,et al.  Constrained Broyden Dimer Method with Bias Potential for Exploring Potential Energy Surface of Multistep Reaction Process. , 2012, Journal of chemical theory and computation.

[28]  R. García-Domenech,et al.  Some new trends in chemical graph theory. , 2008, Chemical reviews.

[29]  S. Yoo,et al.  Global geometry optimization of silicon clusters described by three empirical potentials , 2003 .

[30]  Milan Randic,et al.  Resolution of ambiguities in structure-property studies by use of orthogonal descriptors , 1991, J. Chem. Inf. Comput. Sci..

[31]  P. Bandyopadhyay.  Two-surface Monte Carlo with basin hopping: quantum mechanical trajectory and multiple stationary points of water cluster. , 2008, The Journal of chemical physics.

[32]  Michael K. Gilson,et al.  Tork: Conformational analysis method for molecules and complexes , 2003, J. Comput. Chem..

[33]  Travis E. Oliphant,et al.  Python for Scientific Computing , 2007, Computing in Science & Engineering.

[34]  H. Scheraga,et al.  Global optimization of clusters, crystals, and biomolecules. , 1999, Science.

[35]  David J. Wales,et al.  Global minima of protonated water clusters , 2000 .

[36]  J. Doye,et al.  Tetrahedral global minimum for the 98-atom Lennard-Jones cluster. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[37]  Hyuck-Mo Lee,et al.  New algorithm in the basin hopping Monte Carlo to find the global minimum structure of unary and binary metallic nanoclusters. , 2008, The Journal of chemical physics.

[38]  David J Wales,et al.  Potential energy and free energy landscapes. , 2006, The journal of physical chemistry. B.

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

[40]  A. Kerber,et al.  MOLGEN+, a generator of connectivity isomers and stereoisomers for molecular structure elucidation , 1995 .

[41]  Jinlan Wang,et al.  Optimally stuffed fullerene structures of silicon nanoclusters , 2005 .

[42]  T. Walsh,et al.  Structural isomers and reactivity for Rh6 and Rh6+. , 2006, The journal of physical chemistry. B.

[43]  Gábor Náray-Szabó,et al.  Multicoordinate driven method for approximating enzymatic reaction paths: automatic definition of the reaction coordinate using a subset of chemical coordinates. , 2006, The journal of physical chemistry. A.

[44]  Chang-Yu Hu,et al.  Principles for structure generation of organic isomers from molecular formula , 1994 .

[45]  Satoshi Maeda,et al.  Global mapping of equilibrium and transition structures on potential energy surfaces by the scaled hypersphere search method: applications to ab initio surfaces of formaldehyde and propyne molecules. , 2005, The journal of physical chemistry. A.

[46]  Satoshi Maeda,et al.  Toward Predicting Full Catalytic Cycle Using Automatic Reaction Path Search Method: A Case Study on HCo(CO)3-Catalyzed Hydroformylation. , 2012, Journal of chemical theory and computation.

[47]  Zhi-Pan Liu,et al.  Stochastic Surface Walking Method for Structure Prediction and Pathway Searching. , 2013, Journal of chemical theory and computation.

[48]  Thomas Frauenheim,et al.  Parameter Calibration of Transition-Metal Elements for the Spin-Polarized Self-Consistent-Charge Density-Functional Tight-Binding (DFTB) Method:  Sc, Ti, Fe, Co, and Ni. , 2007, Journal of chemical theory and computation.

[49]  P. Jurs,et al.  Chemical Applications of Graph Theory: Part I. Fundamentals and Topological Indices. , 1988 .

[50]  Halim Kusumaatmaja,et al.  A Local Rigid Body Framework for Global Optimization of Biomolecules. , 2012, Journal of chemical theory and computation.

[51]  D. Wales.  Energy landscapes and properties of biomolecules , 2005, Physical biology.

[52]  A. Baranger,et al.  Variable regiochemistry in the stoichiometric and catalytic hydroamination of alkynes by imidozirconium complexes caused by an unusual dependence of the rate law on alkyne structure and temperature , 1993 .

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

[54]  P. Jurs,et al.  Chemical applications of graph theory. Part II. Isomer enumeration , 1988 .

[55]  Thanh N. Truong,et al.  Application of Chemical Graph Theory for Automated Mechanism Generation , 2003, J. Chem. Inf. Comput. Sci..

[56]  Giulia Rossi,et al.  Electronic and structural shell closure in AgCu and AuCu nanoclusters. , 2006, The journal of physical chemistry. B.

[57]  Masao Iwamatsu,et al.  Basin hopping with occasional jumping , 2004 .

[58]  T. Frauenheim,et al.  DFTB+, a sparse matrix-based implementation of the DFTB method. , 2007, The journal of physical chemistry. A.

[59]  M. Randic.  Characterization of molecular branching , 1975 .

[60]  Satoshi Maeda,et al.  Global reaction route mapping on potential energy surfaces of formaldehyde, formic acid, and their metal-substituted analogues. , 2006, The journal of physical chemistry. A.

[61]  S. Lago,et al.  Structure and stability of small TiO2 nanoparticles. , 2005, The journal of physical chemistry. B.

[62]  Woo Youn Kim,et al.  Noncovalent Interactions of DNA Bases with Naphthalene and Graphene. , 2013, Journal of chemical theory and computation.

[63]  O. Sǐnanoğlu.  Theory of chemical reaction networks. All possible mechanisms or synthetic pathways with given number of reaction steps or species , 1975 .

[64]  Alexandru T. Balaban,et al.  Applications of graph theory in chemistry , 1985, J. Chem. Inf. Comput. Sci..

[65]  Howard R. Mayne,et al.  An investigation of two approaches to basin hopping minimization for atomic and molecular clusters , 1998 .

[66]  Lionello Pogliani,et al.  Plot methods in quantitative structure–property studies , 2004 .

[67]  Zoe E. Brain,et al.  Optimization of a genetic algorithm for searching molecular conformer space. , 2011, The Journal of chemical physics.

[68]  Wing-Ki Liu,et al.  Asynchronous multicanonical basin hopping method and its application to cobalt nanoclusters. , 2005, The Journal of chemical physics.

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

[70]  Nicholas A Besley,et al.  Structural optimization of molecular clusters with density functional theory combined with basin hopping. , 2012, The Journal of chemical physics.

[71]  Timm Lankau,et al.  A constrained reduced-dimensionality search algorithm to follow chemical reactions on potential energy surfaces. , 2013, The Journal of chemical physics.

[72]  David J. Wales,et al.  Global minima of water clusters (H2O)n, n≤21, described by an empirical potential , 1998 .