A Hessian-Based Interpolation Method for Locating Transition Structures

We propose a Hessian-based interpolation scheme for generating estimates of first-order saddle points on potential energy surfaces. Given reactant and product structures and associated Hessians, two molecular-specific force-field-inspired energy functions are constructed, and an estimate of the connecting first-order saddle point is calculated as a minimum on the intersection of the two diabatic energy surfaces. The minimum-energy crossing point can be located starting directly from the reactant or product geometry. At the intersection minimum, an approximate Hessian can be constructed that is suitable for initiating a search for the actual saddle point. The construction of the diabatic energy functions and interpolation procedure is computationally efficient and allows for an automated screening of thousands of potential transition structures connecting conformationally flexible reactants and products, as well as creating global reaction networks. For a benchmark set of 262 reactions, the HINT procedure can locate 90% of the transition structures in an automated workflow.