Graph-based descriptors, such as bond-order matrices and adjacency matrices, offer a simple and compact way of categorizing molecular structures; furthermore, such descriptors can be readily used to catalog chemical reactions (i.e., bond-making and -breaking). As such, a number of graph-based methodologies have been developed with the goal of automating the process of generating chemical reaction network models describing the possible mechanistic chemistry in a given set of reactant species. Here, we outline the evolution of these graph-based reaction discovery schemes, with particular emphasis on more recent methods incorporating graph-based methods with semiempirical and ab initio electronic structure calculations, minimum-energy path refinements, and transition state searches. Using representative examples from homogeneous catalysis and interstellar chemistry, we highlight how these schemes increasingly act as "virtual reaction vessels" for interrogating mechanistic questions. Finally, we highlight where challenges remain, including issues of chemical accuracy and calculation speeds, as well as the inherent challenge of dealing with the vast size of accessible chemical reaction space.
The reaction path Hamiltonian (RPH) can be used to calculate chemical reaction rate constants, going beyond transition-state theory in taking account of recrossing by providing an approximation to the dynamic transmission coefficient. However, the RPH necessitates the calculation of the Hessian matrix at a number of points along the minimum energy path; the associated computational cost stands as a bottleneck in RPH calculations, especially if one is interested in using high-accuracy electronic structure methods. In this work, four different Hessian update schemes (symmetric rank-1, Powell-symmetric Broyden, Bofill, and TS-BFGS updates) are assessed to see whether or not they reliably reproduce calculated transmission coefficients for three different chemical reactions. Based on the reactions investigated, the symmetric rank-1 Hessian update was the least appropriate for RPH construction, giving different transmission coefficients from the standard analytical Hessian approach, as well as inconsistent frequencies and coupling properties. The Bofill scheme, the Powell-symmetric Broyden scheme, and the TS-BFGS scheme were the most reliable Hessian update methods, with transmission coefficients that were in good agreement with those calculated by the standard RPH calculations. The relative accuracy of the different Hessian update schemes is further rationalized by investigating the approximated Coriolis and curvature coupling terms along the reaction-path, providing insight into when these schemes would be expected to work well. Furthermore, the associated computational cost associated with the RPH calculations was substantially reduced by the tested update schemes. Together, these results provide useful rules-of-thumb for using Hessian update schemes in RPH simulations.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.