Identifying Reaction Pathways via Asymptotic Trajectories
Perturbation theories have proven to be useful for identifying the geometric structure of reactions, but are often limited to regions close to the zeroth order reference. Numerical analysis may also breakdown with increasing dimensionality. In this paper, we revisit the concepts of the reactivity map and the reactivity bands to address these limitations. We introduce a reformulated metric, called the asymptotic trajectory indicator, and an efficient algorithm to obtain reactivity boundaries. We demonstrate that this method has sufficient accuracy to reproduce phase space structures such as turnstiles for a 1D model of the isomerization of ketene in an external field. The asymptotic trajectory indicator can be applied to higher dimensional systems coupled to Langevin baths as we demonstrate for a 3D model of the isomerization of ketene.