Integrating steepest-descent reaction pathways for large molecules

Abstract

Exploring potential energy surfaces of large molecular systems can be quite challenging due to the increased number of nuclear degrees of freedom. Many techniques that are well-suited for small and moderate size systems require diagonalization of the energy second-derivative matrix. Since the cost of this step scales as O(N(atoms)(3)) (where N(atoms) is the number of atomic centers), such methods quickly become infeasible and are eventually rendered cost prohibitive. In this work, the recently developed Euler-based predictor-corrector reaction path integration method [H. P. Hratchian, M. J. Frisch, and H. B. Schlegel, J. Chem. Phys. 133, 224101 (2010)] is enhanced and proposed as a useful alternative to conventional reaction path following schemes in studies on very large systems. Because this integrator does not require Hessian diagonalization, the O(N(atoms)(3)) bottleneck afflicting other approaches is completely avoided. The effectiveness of the integrator in large system studies is demonstrated with an enzyme-catalyzed reaction employing an ONIOM (QM:MM) model chemistry and involving 5368 atomic centers.