Bifurcations of dividing surfaces in chemical reactions

Abstract

We study the dynamical behavior of the unstable periodic orbit (NHIM) associated to the non-return transition state (TS) of the H(2) + H collinear exchange reaction and their effects on the reaction probability. By means of the normal form of the Hamiltonian in the vicinity of the phase space saddle point, we obtain explicit expressions of the dynamical structures that rule the reaction. Taking advantage of the straightforward identification of the TS in normal form coordinates, we calculate the reaction probability as a function of the system energy in a more efficient way than the standard Monte Carlo method. The reaction probability values computed by both methods are not in agreement for high energies. We study by numerical continuation the bifurcations experienced by the NHIM as the energy increases. We find that the occurrence of new periodic orbits emanated from these bifurcations prevents the existence of a unique non-return TS, so that for high energies, the transition state theory cannot be longer applied to calculate the reaction probability.