Systematic exploration of the mechanism of chemical reactions: the global reaction route mapping (GRRM) strategy using the ADDF and AFIR methods.

Global reaction route mapping (GRRM), a fully-automated search for all important reaction pathways relevant to a given purpose, on the basis of quantum chemical calculations enables systematic elucidation of complex chemical reaction mechanisms. However, GRRM had previously been limited to very simple systems. This is mainly because such calculations are highly demanding even in small systems when a brute-force sampling is considered. Hence, we have developed two independent but complementary methods: anharmonic downward distortion following (ADDF) and artificial force induced reaction (AFIR) methods. ADDF can follow reaction pathways starting from local minima on the potential energy surface (PES) toward transition structures (TSs) and dissociation channels. AFIR can find pathways starting from two or more reactants toward TSs for their associative reactions. In other words, ADDF searches for A → X type isomerization and A → X + Y type dissociation pathways, whereas AFIR finds A + B → X (+ Y) type associative pathways. Both follow special paths called the ADDF path and the AFIR path, and these tend to pass through near TSs of corresponding reaction pathways, giving approximate TSs. Such approximate TSs can easily be re-optimized to corresponding true TSs by standard geometry optimizations. On the basis of these two methods, we have proposed practical strategies of GRRM. The GRRM strategies have been applied to a variety of chemical systems ranging from thermal- and photochemical-reactions in small systems to organometallic- and enzyme-catalysis, on the basis of quantum chemical calculations. In this perspective, we present an overview of the GRRM strategies and some results of applications. Their practical usage for systematic prediction is also discussed.

[1]  K. Ohno,et al.  D–L Conversion Pathways between Optical Isomers of Alanine: Applications of the Scaled Hypersphere Search Method to Explore Unknown Reaction Routes in a Chiral System , 2006 .

[2]  Satoshi Maeda,et al.  A new method for constructing multidimensional potential energy surfaces by a polar coordinate interpolation technique , 2003 .

[3]  K. Ohno,et al.  Ab initio Studies on Synthetic Routes of Glycine from Simple Molecules via Ammonolysis of Acetolactone: Applications of the Scaled Hypersphere Search Method , 2004 .

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

[5]  J. Kästner.  Umbrella sampling , 2011 .

[6]  K. Dill,et al.  The protein folding problem. , 1993, Annual review of biophysics.

[7]  M. Frisch,et al.  Using redundant internal coordinates to optimize equilibrium geometries and transition states , 1996, J. Comput. Chem..

[8]  Hajime Hirao,et al.  Transition-state optimization by the free energy gradient method : Application to aqueous-phase Menshutkin reaction between ammonia and methyl chloride , 2001 .

[9]  M. Tachikawa,et al.  Theoretical study on mechanisms of structural rearrangement and ionic dissociation in the HCl(H2O)4 cluster with path-integral molecular dynamics simulations , 2011 .

[10]  Frank Jensen,et al.  Gradient extremal bifurcation and turning points: An application to the H2CO potential energy surface , 1996 .

[11]  Paul N. Mortenson,et al.  Energy landscapes: from clusters to biomolecules , 2007 .

[12]  M. Polanyi,et al.  Further considerations on the thermodynamics of chemical equilibria and reaction rates , 1936 .

[13]  Keiji Morokuma,et al.  Theoretical studies of carbonyl photochemistry. I. ab initio potential energy surfaces for the photodissociation H2CO*→H + HCO , 1972 .

[14]  C. J. Tsai,et al.  Use of an eigenmode method to locate the stationary points on the potential energy surfaces of selected argon and water clusters , 1993 .

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

[16]  Karl K. Irikura,et al.  Predicting Unexpected Chemical Reactions by Isopotential Searching , 2000 .

[17]  Keiji Morokuma,et al.  Ab initio Molecular Orbital Studies of Catalytic Elementary Reactions and Catalytic Cycles of Transition-Metal Complexes , 1991 .

[18]  Satoshi Maeda,et al.  An Automated and Systematic Transition Structure Explorer in Large Flexible Molecular Systems Based on Combined Global Reaction Route Mapping and Microiteration Methods. , 2009, Journal of chemical theory and computation.

[19]  J. Doye,et al.  Surveying a potential energy surface by eigenvector-following , 1997 .

[20]  George S. Hammond,et al.  A Correlation of Reaction Rates , 1955 .

[21]  Thom Vreven,et al.  Geometry optimization with QM/MM, ONIOM, and other combined methods. I. Microiterations and constraints , 2003, J. Comput. Chem..

[22]  Kazuhiro Ishida,et al.  The intrinsic reaction coordinate. An ab initio calculation for HNC→HCN and H−+CH4→CH4+H− , 1977 .

[23]  Peter Pulay,et al.  The calculation of ab initio molecular geometries: efficient optimization by natural internal coordinates and empirical correction by offset forces , 1992 .

[24]  H. Schlegel,et al.  Optimization of equilibrium geometries and transition structures , 1982 .

[25]  James W. McIver,et al.  Structure of transition states in organic reactions. General theory and an application to the cyclobutene-butadiene isomerization using a semiempirical molecular orbital method , 1972 .

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

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

[28]  David J. Wales,et al.  Exploring potential energy surfaces with transition state calculations , 1990 .

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

[30]  K. Ohno,et al.  Conversion pathways between a fullerene and a ring among C(20) clusters by a sphere contracting walk method: Remarkable difference in local potential energy landscapes around the fullerene and the ring. , 2006, The Journal of chemical physics.

[31]  W. Lipscomb,et al.  The synchronous-transit method for determining reaction pathways and locating molecular transition states , 1977 .

[32]  K. Morokuma,et al.  No Straight Path: Roaming in Both Ground- and Excited-State Photolytic Channels of NO3 → NO + O2 , 2012, Science.

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

[34]  Xin Li,et al.  The ONIOM method: its foundation and applications to metalloenzymes and photobiology , 2012 .

[35]  Satoshi Maeda,et al.  Automated global mapping of minimal energy points on seams of crossing by the anharmonic downward distortion following method: a case study of H2CO. , 2009, The journal of physical chemistry. A.

[36]  A. Chakraborty,et al.  A growing string method for determining transition states: comparison to the nudged elastic band and string methods. , 2004, The Journal of chemical physics.

[37]  Jaroslav Koča,et al.  VADER: New Software for Exploring Interconversions on Potential Energy Surfaces , 1999, J. Chem. Inf. Comput. Sci..

[38]  K. Morokuma,et al.  Automated Exploration of Photolytic Channels of HCOOH: Conformational Memory via Excited-State Roaming. , 2012, The journal of physical chemistry letters.

[39]  H. Schlegel,et al.  Using Hessian Updating To Increase the Efficiency of a Hessian Based Predictor-Corrector Reaction Path Following Method. , 2005, Journal of chemical theory and computation.

[40]  Satoshi Maeda,et al.  A systematic study on the RuHCl-BINAP-catalyzed asymmetric hydrogenation mechanism by the global reaction route mapping method , 2010 .

[41]  K. Murata,et al.  Applications of homogeneous water-gas shift reaction. III. A further study of the hydrocarbonylation. A highly selective formation of diethyl ketone from ethene, CO, and H2O. , 1981 .

[42]  E. Carter,et al.  Ridge method for finding saddle points on potential energy surfaces , 1993 .

[43]  Satoshi Maeda,et al.  A scaled hypersphere search method for the topography of reaction pathways on the potential energy surface , 2004 .

[44]  Peter Pulay,et al.  Systematic AB Initio Gradient Calculation of Molecular Geometries, Force Constants, and Dipole Moment Derivatives , 1979 .

[45]  D. Heidrich,et al.  Searching for saddle points of potential energy surfaces by following a reduced gradient , 1998 .

[46]  Shigehiko Hayashi,et al.  QM/MM Reweighting Free Energy SCF for Geometry Optimization on Extensive Free Energy Surface of Enzymatic Reaction. , 2012, Journal of chemical theory and computation.

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

[48]  M. Karplus,et al.  Protein Folding: A Perspective from Theory and Experiment. , 1998, Angewandte Chemie.

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

[50]  K. Müller,et al.  Location of saddle points and minimum energy paths by a constrained simplex optimization procedure , 1979 .

[51]  W. Heisenberg,et al.  Zur Quantentheorie der Molekeln , 1924 .

[52]  Satoshi Maeda,et al.  Automated exploration of reaction channels , 2008 .

[53]  F. Jensen,et al.  Searching Peptide Conformational Space. , 2011, Journal of chemical theory and computation.

[54]  Satoshi Maeda,et al.  Finding reaction pathways for multicomponent reactions: the Passerini reaction is a four-component reaction. , 2011, Angewandte Chemie.

[55]  Fumio Hirata,et al.  Reference interaction site model self-consistent field study for solvation effect on carbonyl compounds in aqueous solution , 1994 .

[56]  Satoshi Maeda,et al.  Insight into global reaction mechanism of [C2, H4, O] system from ab initio calculations by the scaled hypersphere search method. , 2007, The journal of physical chemistry. A.

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

[58]  K. Ohno,et al.  Ab initio anharmonic calculations of vibrational frequencies of benzene by means of efficient construction of potential energy functions , 2011 .

[59]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[60]  Satoshi Maeda,et al.  Exploring Multiple Potential Energy Surfaces: Photochemistry of Small Carbonyl Compounds , 2012 .

[61]  Josep Maria Bofill,et al.  Finding transition states using reduced potential-energy surfaces , 2001 .

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

[63]  K. Fukui,et al.  Differential geometry of chemically reacting systems , 1978 .

[64]  K. Morokuma,et al.  Excited-State Roaming Dynamics in Photolysis of a Nitrate Radical , 2011 .

[65]  Jochen Autschbach,et al.  Theoretical methods of potential use for studies of inorganic reaction mechanisms. , 2005, Chemical reviews.

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

[67]  S. Niu,et al.  Theoretical studies on reactions of transition-metal complexes. , 2000, Chemical reviews.

[68]  Fumio Hirata,et al.  Analytical energy gradient for the reference interaction site model multiconfigurational self‐consistent‐field method: Application to 1,2‐difluoroethylene in aqueous solution , 1996 .

[69]  K. Ohno,et al.  A new approach for finding a transition state connecting a reactant and a product without initial guess: applications of the scaled hypersphere search method to isomerization reactions of HCN, (H2O)2, and alanine dipeptide , 2005 .

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

[71]  Satoshi Maeda,et al.  Communications: A systematic method for locating transition structures of A+B-->X type reactions. , 2010, The Journal of chemical physics.

[72]  Paul Ha-Yeon Cheong,et al.  Computational prediction of small-molecule catalysts , 2008, Nature.

[73]  Satoshi Maeda,et al.  A new global reaction route map on the potential energy surface of H2CO with unrestricted level , 2008 .

[74]  H. Bernhard Schlegel,et al.  Exploring potential energy surfaces for chemical reactions: An overview of some practical methods , 2003, J. Comput. Chem..

[75]  Thom Vreven,et al.  Geometry optimization with QM/MM methods II: Explicit quadratic coupling , 2006 .

[76]  Josep Maria Bofill,et al.  Updated Hessian matrix and the restricted step method for locating transition structures , 1994, J. Comput. Chem..

[77]  Satoshi Maeda,et al.  Structures of water octamers (H2O)8: exploration on ab initio potential energy surfaces by the scaled hypersphere search method. , 2007, The journal of physical chemistry. A.

[78]  Tetsuya Taketsugu,et al.  Exploring Potential Energy Surfaces of Large Systems with Artificial Force Induced Reaction Method in Combination with ONIOM and Microiteration. , 2012, Journal of chemical theory and computation.

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

[80]  E. Vanden-Eijnden,et al.  String method for the study of rare events , 2002, cond-mat/0205527.

[81]  W. Quapp,et al.  Reaction channels of the potential energy surface: application of Newton trajectories , 2004 .

[82]  Ajit Banerjee,et al.  Search for stationary points on surfaces , 1985 .

[83]  H. Bernhard Schlegel,et al.  Following gradient extremal paths , 1992 .

[84]  D. Sánchez-Portal,et al.  The SIESTA method for ab initio order-N materials simulation , 2001, cond-mat/0111138.

[85]  K. Ohno,et al.  Microsolvation of hydrogen sulfide: exploration of H2S.(H2O)n and SH-.H3O+.(H2O)n-1 (n = 5-7) cluster structures on ab initio potential energy surfaces by the scaled hypersphere search method. , 2008, The journal of physical chemistry. A.

[86]  G. Frenking,et al.  Theoretical Studies of Some Transition-Metal-Mediated Reactions of Industrial and Synthetic Importance. , 2000, Chemical reviews.

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

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

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

[90]  R. Bell,et al.  The Theory of Reactions Involving Proton Transfers , 1936 .

[91]  E. Corey,et al.  Robert Robinson Lecture. Retrosynthetic thinking—essentials and examples , 1988 .

[92]  J Andrew McCammon,et al.  Studying functional dynamics in bio-molecules using accelerated molecular dynamics. , 2011, Physical chemistry chemical physics : PCCP.

[93]  P. Culot,et al.  A quasi-Newton algorithm for first-order saddle-point location , 1992 .

[94]  K. Ohno,et al.  No activation barrier synthetic route of glycine from simple molecules (NH3, CH2, and CO2) via carboxylation of ammonium ylide: a theoretical study by the scaled hypersphere search method , 2004 .

[95]  K. Morokuma,et al.  Finding Minimum Structures on the Seam of Crossing in Reactions of Type A + B → X: Exploration of Nonadiabatic Ignition Pathways of Unsaturated Hydrocarbons. , 2011, The journal of physical chemistry letters.

[96]  Christodoulos A. Floudas,et al.  Locating all transition states and studying the reaction pathways of potential energy surfaces , 1999 .

[97]  Bruce A. Murtagh,et al.  Computational Experience with Quadratically Convergent Minimisation Methods , 1970, Comput. J..

[98]  Keiji Morokuma,et al.  Photodissociation of formaldehyde: Potential energy surfaces for H2CO → H2 + CO , 1974 .

[99]  A. Suits,et al.  Roaming Radical Reactions , 2011 .

[100]  Yi Luo,et al.  Automated exploration of stable isomers of H+(H2O)n (n = 5–7) via ab initio calculations: An application of the anharmonic downward distortion following algorithm , 2009, J. Comput. Chem..

[101]  Satoshi Maeda,et al.  Global reaction route mapping on potential energy surfaces of C2H7+ and C3H9+ , 2007 .

[102]  Paul W Ayers,et al.  Dual Grid Methods for Finding the Reaction Path on Reduced Potential Energy Surfaces. , 2010, Journal of chemical theory and computation.

[103]  FRANCESCO AQUILANTE,et al.  MOLCAS 7: The Next Generation , 2010, J. Comput. Chem..

[104]  Satoshi Maeda,et al.  Automated exploration of adsorption structures of an organic molecule on RuH2-BINAP by the ONIOM method and the scaled hypersphere search method. , 2007, The journal of physical chemistry. A.

[105]  Y. Abashkin,et al.  Transition state structures and reaction profiles from constrained optimization procedure. Implementation in the framework of density functional theory , 1994 .

[106]  Kazuhiro Ishida,et al.  Efficient determination and characterization of transition states using ab-initio methods , 1977 .

[107]  M. Shiga,et al.  Path-integral molecular dynamics simulations of hydrated hydrogen chloride cluster HCl(H2O)4 on a semiempirical potential energy surface , 2009 .

[108]  H. Grubmüller,et al.  Predicting unimolecular chemical reactions: Chemical flooding , 2002 .

[109]  Hao Hu,et al.  Free energies of chemical reactions in solution and in enzymes with ab initio quantum mechanics/molecular mechanics methods. , 2008, Annual review of physical chemistry.

[110]  Michael Hirsch,et al.  Determination of energy minima and saddle points using multireference configuration interaction methods in combination with reduced gradient following: The S0 surface of H2CO and the T1 and T2 surfaces of acetylene , 2002, J. Comput. Chem..

[111]  Satoshi Maeda,et al.  Lowest transition state for the chirality-determining step in Ru((R)-BINAP)-catalyzed asymmetric hydrogenation of methyl-3-oxobutanoate. , 2008, Journal of the American Chemical Society.

[112]  Ron Elber,et al.  A method for determining reaction paths in large molecules: application to myoglobin , 1987 .

[113]  K. Morokuma,et al.  Experimental and theoretical investigations of isomerization reactions of ionized acetone and its dimer. , 2012, Physical chemistry chemical physics : PCCP.

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

[115]  Josep Maria Bofill,et al.  A new look at the reduced-gradient-following path , 2002 .

[116]  Satoshi Maeda,et al.  Global investigation on the potential energy surface of CH3CN: application of the scaled hypersphere search method. , 2005, The journal of physical chemistry. A.

[117]  K. Ohno,et al.  Theoretical investigation of the reaction pathway of O atom on Si(001)-(2×1) , 2010 .

[118]  Poul Jørgensen,et al.  Proper characterization of MC SCF stationary points , 1983 .

[119]  M. V. Basilevsky.  The topography of potential energy surfaces , 1982 .

[120]  J. Pancíř.  Calculation of the least energy path on the energy hypersurface , 1975 .

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

[122]  K. Ohno,et al.  Finding important anharmonic terms in the sixth-order potential energy function by the scaled hypersphere search method: an application to vibrational analyses of molecules and clusters. , 2008, The Journal of chemical physics.

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

[124]  Peter Pulay,et al.  Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules , 1969 .

[125]  K. Ohno,et al.  Generation Mechanisms of Amino Acids in Interstellar Space via Reactions between Closed-Shell Species: Significance of Higher Energy Isomers in Molecular Evolution , 2006 .

[126]  P. Jørgensen,et al.  A gradient extremal walking algorithm , 1988 .

[127]  H. Bernhard Schlegel,et al.  An improved algorithm for reaction path following , 1989 .

[128]  K. Ohno,et al.  Intramolecular vibrational frequencies of water clusters (H2O)n (n=2-5): Anharmonic analyses using potential functions based on the scaled hypersphere search method. , 2008, The Journal of chemical physics.

[129]  Feliu Maseras,et al.  IMOMM: A new integrated ab initio + molecular mechanics geometry optimization scheme of equilibrium structures and transition states , 1995, J. Comput. Chem..

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

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

[132]  Michael A. Collins,et al.  Molecular potential-energy surfaces for chemical reaction dynamics , 2002 .

[133]  T. Helgaker.  Transition-state optimizations by trust-region image minimization , 1991 .

[134]  H. Bernhard Schlegel,et al.  Methods for optimizing large molecules. II. Quadratic search , 1999 .

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

[136]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.