Quasi-continuous Interpolation Scheme for Pathways between Distant Configurations

Abstract

A quasi-continuous interpolation (QCI) scheme is introduced for characterizing physically realistic initial pathways from which to initiate transition state searches and construct kinetic transition networks. Applications are presented for peptides, proteins, and a morphological transformation in an atomic cluster. The first step in each case involves end point alignment, and we describe the use of a shortest augmenting path algorithm for optimizing permutational isomers. The QCI procedure then employs an interpolating potential, which preserves the covalent bonding framework for the biomolecules and includes repulsive terms between unconstrained atoms. This potential is used to identify an interpolating path by minimizing contributions from a connected set of images, including terms corresponding to minima in the interatomic distances between them. This procedure detects unphysical geometries in the line segments between images. The most difficult cases, where linear interpolation would involve chain crossings, are treated by growing the structure an atom at a time using the interpolating potential. To test the QCI procedure, we carry through a series of benchmark calculations where the initial interpolation is coupled to explicit transition state searches to produce complete pathways between specified local minima.This work was supported by the Engineering and Physical Sciences Research Council [grant number EP/H042660/1]This document is the unedited Author’s version of a Submitted Work that was subsequently accepted for publication in the Journal of Chemical Theory and Computation, copyright © American Chemical Society after peer review. To access the final edited and published work see http://dx.doi.org/10.1021/ct300483