A reduced two dimensional model is used to study Ketene isomerization reaction. In light of recent results by Ulusoy textit{et al.} [J. Phys. Chem. A {bf 117}, 7553 (2013)], the present work focuses on the generalization of the roaming mechanism to the Ketene isomerization reaction by applying our phase space approach previously used to elucidate the roaming phenomenon in ion-molecule reactions. Roaming is again found be associated with the trapping of trajectories in a phase space region between two dividing surfaces; trajectories are classified as reactive or nonreactive, and are further naturally classified as direct or non-direct (roaming). The latter long-lived trajectories are trapped in the region of non-linear mechanical resonances, which in turn define alternative reaction pathways in phase space. It is demonstrated that resonances associated with periodic orbits provide a dynamical explanation of the quantum mechanical resonances found in the isomerization rate constant calculations by Gezelter and Miller [J. Chem. Phys. {bf 103}, 7868-7876 (1995)]. Evidence of the trapping of trajectories by `sticky resonant periodic orbits is provided by plotting Poincare surfaces of section, and a gap time analysis is carried out in order to investigate the statistical assumption inherent in transition state theory for Ketene isomerization.
We study reaction dynamics on a model potential energy surface exhibiting post-transition state bifurcation in the vicinity of a valley ridge inflection point. We compute fractional yields of products reached after the VRI region is traversed, both with and without dissipation. It is found that apparently minor variations in the potential lead to significant changes in the reaction dynamics. Moreover, when dissipative effects are incorporated, the product ratio depends in a complicated and highly non-monotonic fashion on the dissipation parameter. Dynamics in the vicinity of the VRI point itself play essentially no role in determining the product ratio, except in the highly dissipative regime.
We provide a dynamical interpretation of the recently identified `roaming mechanism for molecular dissociation reactions in terms of geometrical structures in phase space. These are NHIMs (Normally Hyperbolic Invariant Manifolds) and their stable/unstable manifolds that define transition states for ion-molecule association or dissociation reactions. The associated dividing surfaces rigorously define a roaming region of phase space, in which both reactive and nonreactive trajectories can be trapped for arbitrarily long times.
