Mode-Tracking Based Stationary-Point Optimization

Abstract

In this work, we present a transition-state optimization protocol based on the Mode-Tracking algorithm [J. Chem. Phys. 118 (2003) 1634]. By calculating only the eigenvector of interest instead of diagonalizing the full Hessian matrix and performing an eigenvector following search based on the selectively calculated vector, we can efficiently optimize transition-state structures. The initial guess structures and eigenvectors are either chosen from a linear interpolation between the reactant and product structures, from a nudged-elastic band search, from a constrained-optimization scan, or from the minimum-energy structures. Alternatively, initial guess vectors based on chemical intuition may be defined. We then iteratively refine the selected vectors by the Davidson subspace iteration technique. This procedure accelerates finding transition states for large molecules of a few hundred atoms. It is also beneficial in cases where the starting structure is very different from the transition-state structure or where the desired vector to follow is not the one with lowest eigenvalue. Explorative studies of reaction pathways are feasible by following manually constructed molecular distortions.

28 pages, 5 figures

References in corpus (2)

• Interactive Chemical Reactivity Exploration
• Real-time Quantum Chemistry

Cited by in corpus (7)

• Exploration of Reaction Pathways and Chemical Transformation Networks
• Heuristics-Guided Exploration of Reaction Mechanisms
• Context-Driven Exploration of Complex Chemical Reaction Networks
• Systematic Error Estimation for Chemical Reaction Energies
• Minimum Energy Paths and Transition States by Curve Optimization
• Expansive Quantum Mechanical Exploration of Chemical Reaction Paths
• Molecular propensity as a driver for explorative reactivity studies