A periodic orbit bifurcation analysis of vibrationally excited isotopologues of sulfur dioxide and water molecules: symmetry breaking substitutions

Fundamental studies of vibrational resonance phenomena by multivalued resummation of the divergent Rayleigh-Schrödinger perturbation theory series: deciphering polyad structures of three H216O isotopologues

A fundamental quantitative study of the vibrational resonances of three H₂¹⁶O H,D-isotopologues with a quartic Watson Hamiltonian was carried out using the resummation of the high-order (∼λⁿ, n ≤ 203) divergent Rayleigh-Schrödinger perturbation theory (RSPT) series by quartic Padé-Hermite multivalued diagonal approximants PH^[m,m,m,m,m], m ≤ 40. The resonance condition between a pair of states is formulated as the existence of a common complex energy solution branch point inside the unit circle: |E(λ_j)|,  Re(λ_j)² + Im(λ_j)² ≤ 1. For the matrix formulation of the vibrational problem (VCI), the existence of common branch points is governed by the Katz theorem and they can be found as roots of discriminant polynomials. The main branches of the Padé-Hermite approximants typically reproduce VCI energies with high accuracy while alternative branches often fit nearby resonant states. The resummation of the RSPT series for H₂O and D₂O (up to the tenth polyad) revealed not only Fermi and Darling-Dennison resonances, but also unusual (0,2,-5) and (5,-2,0) resonance effects matching the (5,2,5) polyad structure, while the (3,2,1) structure was rigorously confirmed for HDO. It is demonstrated that the (5,2,5) polyad structure ensures good organization of high-energy resonating states and breaks down the classic (2,1,2) structure. The advocated methodology of quantitative description of resonance phenomena and revealing polyad structures is suitable for larger molecules and can be adapted to linear molecules and symmetric tops. Its application ensures rigorous classification of vibrational states and can be used in quantitative vibration-rotation spectroscopy.

Roaming dynamics in ion-molecule reactions: phase space reaction pathways and geometrical interpretation

A model Hamiltonian for the reaction CH4(+) -> CH3(+) + H, parametrized to exhibit either early or late inner transition states, is employed to investigate the dynamical characteristics of the roaming mechanism. Tight/loose transition states and conventional/roaming reaction pathways are identified in terms of time-invariant objects in phase space. These are dividing surfaces associated with normally hyperbolic invariant manifolds (NHIMs). For systems with two degrees of freedom NHIMS are unstable periodic orbits which, in conjunction with their stable and unstable manifolds, unambiguously define the (locally) non-recrossing dividing surfaces assumed in statistical theories of reaction rates. By constructing periodic orbit continuation/bifurcation diagrams for two values of the potential function parameter corresponding to late and early transition states, respectively, and using the total energy as another parameter, we dynamically assign different regions of phase space to reactants and products as well as to conventional and roaming reaction pathways. The classical dynamics of the system are investigated by uniformly sampling trajectory initial conditions on the dividing surfaces. Trajectories are classified into four different categories: direct reactive and non-reactive trajectories, which lead to the formation of molecular and radical products respectively, and roaming reactive and non-reactive orbiting trajectories, which represent alternative pathways to form molecular and radical products. By analysing gap time distributions at several energies, we demonstrate that the phase space structure of the roaming region, which is strongly influenced by nonlinear resonances between the two degrees of freedom, results in nonexponential (nonstatistical) decay.