Defect in the Joint Spectrum of Hydrogen due to Monodromy.

In addition to the well-known case of spherical coordinates, the Schrödinger equation of the hydrogen atom separates in three further coordinate systems. Separating in a particular coordinate system defines a system of three commuting operators. We show that the joint spectrum of the Hamilton operator, the z component of the angular momentum, and an operator involving the z component of the quantum Laplace-Runge-Lenz vector obtained from separation in prolate spheroidal coordinates has quantum monodromy for energies sufficiently close to the ionization threshold. The precise value of the energy above which monodromy is observed depends on the distance of the focus points of the spheroidal coordinates. The presence of monodromy means that one cannot globally assign quantum numbers to the joint spectrum. Whereas the principal quantum number n and the magnetic quantum number m correspond to the Bohr-Sommerfeld quantization of globally defined classical actions a third quantum number cannot be globally defined because the third action is globally multivalued.

The phase space geometry underlying roaming reaction dynamics.

Recent studies have found an unusual way of dissociation in formaldehyde. It can be characterized by a hydrogen atom that separates from the molecule, but instead of dissociating immediately it roams around the molecule for a considerable amount of time and extracts another hydrogen atom from the molecule prior to dissociation. This phenomenon has been coined roaming and has since been reported in the dissociation of a number of other molecules. In this paper we investigate roaming in Chesnavich's CH 4 + model. During dissociation the free hydrogen must pass through three phase space bottleneck for the classical motion, that can be shown to exist due to unstable periodic orbits. None of these orbits is associated with saddle points of the potential energy surface and hence related to transition states in the usual sense. We explain how the intricate phase space geometry influences the shape and intersections of invariant manifolds that form separatrices, and establish the impact of these phase space structures on residence times and rotation numbers. Ultimately we use this knowledge to attribute the roaming phenomenon to particular heteroclinic intersections.

Reaction dynamics through kinetic transition states.

The transformation of a system from one state to another is often mediated by a bottleneck in the system's phase space. In chemistry, these bottlenecks are known as transition states through which the system has to pass in order to evolve from reactants to products. The chemical reactions are usually associated with configurational changes where the reactants and products states correspond, e.g., to two different isomers or the undissociated and dissociated state of a molecule or cluster. In this Letter, we report on a new type of bottleneck which mediates kinetic rather than configurational changes. The phase space structures associated with such kinetic transition states and their dynamical implications are discussed for the rotational vibrational motion of a triatomic molecule. An outline of more general related phase space structures with important dynamical implications is given.

The flux-flux correlation function for anharmonic barriers.

The flux-flux correlation function formalism is a standard and widely used approach for the computation of reaction rates. In this paper we introduce a method to compute the classical and quantum flux-flux correlation functions for anharmonic barriers essentially analytically through the use of the classical and quantum normal forms. In the quantum case we show that for a general f degree-of-freedom system having an index one saddle the quantum normal form reduces the computation of the flux-flux correlation function to that of an effective one-dimensional anharmonic barrier. The example of the computation of the quantum flux-flux correlation function for a fourth order anharmonic barrier is worked out in detail, and we present an analytical expression for the quantum mechanical microcanonical flux-flux correlation function. We then give a discussion of the short-time and harmonic limits.

The quantum normal form approach to reactive scattering: The cumulative reaction probability for collinear exchange reactions.

The quantum normal form approach to quantum transition state theory is used to compute the cumulative reaction probability for collinear exchange reactions. It is shown that for heavy-atom systems such as the nitrogen-exchange reaction, the quantum normal form approach gives excellent results and has major computational benefits over full reactive scattering approaches. For light atom systems such as the hydrogen-exchange reaction however, the quantum normal approach is shown to give only poor results. This failure is attributed to the importance of tunneling trajectories in light atom reactions that are not captured by the quantum normal form as indicated by the only very slow convergence of the quantum normal form for such systems.

Microcanonical rates, gap times, and phase space dividing surfaces.

The general approach to classical unimolecular reaction rates due to Thiele is revisited in light of recent advances in the phase space formulation of transition state theory for multidimensional systems. Key concepts, such as the phase space dividing surface separating reactants from products, the average gap time, and the volume of phase space associated with reactive trajectories, are both rigorously defined and readily computed within the phase space approach. We analyze in detail the gap time distribution and associated reactant lifetime distribution for the isomerization reaction HCN <==> CNH, previously studied using the methods of phase space transition state theory. Both algebraic (power law) and exponential decay regimes have been identified. Statistical estimates of the isomerization rate are compared with the numerically determined decay rate. Correcting the RRKM estimate to account for the measure of the reactant phase space region occupied by trapped trajectories results in a drastic overestimate of the isomerization rate. Compensating but as yet not fully understood trapping mechanisms in the reactant region serve to slow the escape rate sufficiently that the uncorrected RRKM estimate turns out to be reasonably accurate, at least at the particular energy studied. Examination of the decay properties of subensembles of trajectories that exit the HCN well through either of two available symmetry related product channels shows that the complete trajectory ensemble effectively attains the full symmetry of the system phase space on a short time scale t approximately < 0.5 ps, after which the product branching ratio is 1:1, the "statistical" value. At intermediate times, this statistical product ratio is accompanied by nonexponential (nonstatistical) decay. We point out close parallels between the dynamical behavior inferred from the gap time distribution for HCN and nonstatistical behavior recently identified in reactions of some organic molecules.

Nonuniqueness of the phase shift in central scattering due to monodromy.

Scattering at a central potential is completely characterized by the phase shifts which are the differences in phase between outgoing scattered and unscattered partial waves. In this Letter, it is shown that, for 2D scattering at a repulsive central potential, the phase shift cannot be uniquely defined due to a topological obstruction which is similar to monodromy in bound systems.

Efficient computation of transition state resonances and reaction rates from a quantum normal form.

A quantum version of a recent formulation of transition state theory in phase space is presented. The theory developed provides an algorithm to compute quantum reaction rates and the associated Gamov-Siegert resonances with very high accuracy. The algorithm is especially efficient for multi-degree-of-freedom systems where other approaches are no longer feasible.

Efficient procedure to compute the microcanonical volume of initial conditions that lead to escape trajectories from a multidimensional potential well.

A procedure is presented for computing the phase space volume of initial conditions for trajectories that escape or "react" from a multidimensional potential well. The procedure combines a phase space transition state theory, which allows one to construct dividing surfaces that are free of local recrossing and that minimize the directional flux, and a classical spectral theorem. The procedure gives the volume of reactive initial conditions in terms of a sum over each entrance channel of the well of the product of the phase space flux across the dividing surface associated with the channel and the mean residence time in the well of trajectories which enter through the channel. This approach is illustrated for HCN isomerization in three dimensions, for which the method is several orders of magnitude more efficient than standard Monte Carlo sampling.

Phase space conduits for reaction in multidimensional systems: HCN isomerization in three dimensions.

The three-dimensional hydrogen cyanide/isocyanide isomerization problem is taken as an example to present a general theory for computing the phase space structures which govern classical reaction dynamics in systems with an arbitrary (finite) number of degrees of freedom. The theory, which is algorithmic in nature, comprises the construction of a dividing surface of minimal flux which is locally a "surface of no return." The theory also allows for the computation of the global phase space transition pathways that trajectories must follow in order to react. The latter are enclosed by the stable and unstable manifolds of a so-called normally hyperbolic invariant manifold (NHIM). A detailed description of the geometrical structures and the resulting constraints on reaction dynamics is given, with particular emphasis on the three degrees of freedom case. A procedure is given which uses these structures to compute orbits homoclinic to, and heteroclinic between, NHIMs. The role of homoclinic and heteroclinic orbits in global recrossings of dividing surfaces and transport in complex systems is explained. The complete description provided here is inherently one within phase space; it cannot be inferred from a configuration space picture. A complexification of the classical phase space structures to incorporate quantum effects is also discussed. The results presented here call into question certain assumptions routinely made on the global dynamics; this paper provides methods that enable one to understand and quantify the phase space dynamics of reactions without making such assumptions.