We develop Wigner's approach to a dynamical transition state theory in phase space in both the classical and quantum mechanical settings. The key to our development is the construction of a normal form for describing the dynamics in the neighborhood of a specific type of saddle point that governs the evolution from reactants to products in high dimensional systems. In the classical case this is the standard Poincaré-Birkhoff normal form. In the quantum case we develop a normal form based on the Weyl calculus and an explicit algorithm for computing this quantum normal form. The classical normal form allows us to discover and compute the phase space structures that govern classical reaction dynamics. From this knowledge we are able to provide a direct construction of an energy dependent dividing surface in phase space having the properties that trajectories do not locally "re-cross" the surface and the directional flux across the surface is minimal. Using this, we are able to give a formula for the directional flux through the dividing surface that goes beyond the harmonic approximation. We relate this construction to the fluxflux autocorrelation function which is a standard ingredient in the expression for the reaction rate in the chemistry community. We also give a classical mechanical interpretation of the activated complex as a normally hyperbolic invariant manifold (NHIM), and further describe the structure of the NHIM. The quantum normal form provides us with an efficient algorithm to compute quantum reaction rates and we relate this algorithm to the quantum version of the flux-flux autocorrelation function formalism. The significance of the classical phase space structures for the quantum mechanics of reactions is elucidated by studying the phase space distribution of scattering states. The quantum normal form also provides an efficient way of computing Gamov-Siegert resonances. We relate these resonances to the lifetimes of the quantum activated complex. We consider several one, two, and three degree-of-freedom systems and show explicitly how calculations of the above quantities can be carried out. Our theoretical framework is valid for Hamiltonian systems with an arbitrary number of degrees of freedom and we demonstrate that in several situations it gives rise to algorithms that are computationally more efficient than existing methods.
The fundamental assumption of conventional transition state theory is the existence of a dividing surface having the property that trajectories originating in reactants must cross the surface only once and then proceed to products. Recently it has been shown [1,2] how to construct a dividing surface in phase space for Hamiltonian systems with an arbitrary (but finite) number of degrees of freedom having the property that trajectories only cross once locally. In this letter we provide an argument showing that the flux across this dividing surface is a minimum with respect to certain types of variations of the dividing surface. In this sense, conventional transition state theory is shown to be equivalent to variational transition state theory. PACS: 82.20.Db, 82.20.Nk,
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.
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 Շ 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.
The exact and semiclassical quantum mechanics of the elliptic billiard is investigated. The classical system is integrable and exhibits a separatrix, dividing the phase space into regions of oscillatory and rotational motion. The classical separability carries over to quantum mechanics, and the Schrödinger equation is shown to be equivalent to the spheroidal wave equation. The quantum eigenvalues show a clear pattern when transformed into the classical action space. The implication of the separatrix on the wave functions is illustrated. A uniform WKB quantization taking into account complex orbits is shown to be adequate for the semiclassical quantization in the presence of a separatrix. The pattern of states in classical action space is nicely explained by this quantization procedure. We extract an effective Maslov phase varying smoothly on the energy surface, which is used to modify the Berry-Tabor trace formula, resulting in a summation over non-periodic orbits. This modified trace formula produces the correct number of states, even close to the separatrix. The Fourier transform of the density of states is explained in terms of classical orbits, and the amplitude and form of the different kinds of peaks is analytically calculated.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.