Phase spacing

It is sometimes very usefiil to look at a trajectory such as the synnnetric or antisynnnetric stretch of figure Al.2.5 and figure A1.2.6 not in the physical spatial coordinates (r. . r y), but in the phase space of Hamiltonian mechanics [16, 29], which in addition to the coordinates (r. . r ) also has as additional coordinates the set of conjugate momenta. . pj. ). In phase space, a one-diniensional trajectory such as the aiitisymmetric stretch again appears as a one-diniensional curve, but now the curve closes on itself Such a trajectory is referred to in nonlinear dynamics as a periodic orbit [29]. One says that the aihiamionic nonnal modes of Moser and Weinstein are stable periodic orbits. 

Figure Al.2.7. Trajectory of two coupled stretches, obtained by integrating Hamilton s equations for motion on a PES for the two modes. The system has stable anhamionic synmretric and antisyimnetric stretch modes, like those illustrated in figrne Al.2.6. In this trajectory, semiclassically there is one quantum of energy in each mode, so the trajectory corresponds to a combination state with quantum numbers nj = [1, 1]. The woven pattern shows that the trajectory is regular rather than chaotic, corresponding to motion in phase space on an invariant torus.

Figure Al.2.12. Energy level pattern of a polyad with resonant collective modes. The top and bottom energy levels conespond to overtone motion along the two modes shown in figure Al.2.11. which have a different frequency. The spacing between adjacent levels decreases until it reaches a minimum between the third and fourth levels from the top. This minimum is the hallmark of a separatrix [29, 45] in phase space.

The question of non-classical manifestations is particularly important in view of the chaos that we have seen is present in the classical dynamics of a multimode system, such as a polyatomic molecule, with more than one resonance coupling. Chaotic classical dynamics is expected to introduce its own peculiarities into quantum spectra [29, 77]. In Fl20, we noted that chaotic regions of phase space are readily seen in the classical dynamics corresponding to the spectroscopic Flamiltonian. Flow important are the effects of chaos in the observed spectrum, and in the wavefiinctions of tire molecule In FI2O, there were some states whose wavefiinctions appeared very disordered, in the region of the... 

In the above discussion of relaxation to equilibrium, the density matrix was implicitly cast in the energy representation. However, the density operator can be cast in a variety of representations other than the energy representation. Two of the most connnonly used are the coordinate representation and the Wigner phase space representation. In addition, there is the diagonal representation of the density operator in this representation, the most general fomi of p takes the fomi... 

The set of microstates of a finite system in quantum statistical mechanics is a finite, discrete denumerable set of quantum states each characterized by an appropriate collection of quantum numbers. In classical statistical mechanics, the set of microstates fonn a continuous (and therefore infinite) set of points in f space (also called phase space). 

WIson K G 1971 Renormalization group and critical phenomena. II. Phase space cell analysis of critical behaviour P/rys. Rev. B 4 3184-205... 

A situation that arises from the intramolecular dynamics of A and completely distinct from apparent non-RRKM behaviour is intrinsic non-RRKM behaviour [9], By this, it is meant that A has a non-random P(t) even if the internal vibrational states of A are prepared randomly. This situation arises when transitions between individual molecular vibrational/rotational states are slower than transitions leading to products. As a result, the vibrational states do not have equal dissociation probabilities. In tenns of classical phase space dynamics, slow transitions between the states occur when the reactant phase space is metrically decomposable [13,14] on the timescale of the imimolecular reaction and there is at least one bottleneck [9] in the molecular phase space other than the one defining the transition state. An intrinsic non-RRKM molecule decays non-exponentially with a time-dependent unimolecular rate constant or exponentially with a rate constant different from that of RRKM theory. 

The hypersurface fomied from variations in the system s coordinates and momenta at//(p, q) = /Tis the microcanonical system s phase space, which, for a Hamiltonian with 3n coordinates, has a dimension of 6n -1. The assumption that the system s states are populated statistically means that the population density over the whole surface of the phase space is unifomi. Thus, the ratio of molecules at the dividing surface to the total molecules [dA(qi, p )/A]... 

In deriving the RRKM rate constant in section A3.12.3.1. it is assumed that the rate at which reactant molecules cross the transition state, in the direction of products, is the same rate at which the reactants fonn products. Thus, if any of the trajectories which cross the transition state in the product direction return to the reactant phase space, i.e. recross the transition state, the actual unimolecular rate constant will be smaller than that predicted by RRKM theory. This one-way crossing of the transition state, witii no recrossmg, is a fiindamental assumption of transition state theory [21]. Because it is incorporated in RRKM theory, this theory is also known as microcanonical transition state theory. 

The intennolecular Hamiltonian of the product fragments is used to calculate the sum of states of the transitional modes, when they are treated as rotations. The resulting model [28] is nearly identical to phase space theory [29],... 

The bulk of the infomiation about anhannonicity has come from classical mechanical calculations. As described above, the aidiannonic RRKM rate constant for an analytic potential energy fiinction may be detemiined from either equation (A3.12.4) [13] or equation (A3.12.24) [46] by sampling a microcanonical ensemble. This rate constant and the one calculated from the hamionic frequencies for the analytic potential give the aidiannonic correctiony j ( , J) in equation (A3.12.41). The transition state s aidiannonic classical sum of states is found from the phase space integral... 

As discussed in section A3.12.2. intrinsic non-RRKM behaviour occurs when there is at least one bottleneck for transitions between the reactant molecule s vibrational states, so drat IVR is slow and a microcanonical ensemble over the reactant s phase space is not maintained during the unimolecular reaction. The above discussion of mode-specific decomposition illustrates that there are unimolecular reactions which are intrinsically non-RRKM. Many van der Waals molecules behave in this maimer [4,82]. For example, in an initial microcanonical ensemble for the ( 211 )2 van der Waals molecule both the C2H4—C2H4 intennolecular modes and C2H4 intramolecular modes are excited with equal probabilities. However, this microcanonical ensemble is not maintained as the dimer dissociates. States with energy in the intermolecular modes react more rapidly than do those with the C2H4 intramolecular modes excited [85]. 

In other words, if we look at any phase-space volume element, the rate of incoming state points should equal the rate of outflow. This requires that be a fiinction of the constants of the motion, and especially Q=Q i). Equilibrium also implies d(/)/dt = 0 for any /. The extension of the above equations to nonequilibriiim ensembles requires a consideration of entropy production, the method of controlling energy dissipation (diennostatting) and the consequent non-Liouville nature of the time evolution [35]. 

The velocity Verlet algorithm may be derived by considering a standard approximate decomposition of the Liouville operator which preserves reversibility and is symplectic (which implies that volume in phase space is conserved). This approach [47] has had several beneficial consequences. 

Hoover W G 1985 Canonical dynamics equilibrium phase-space distributions Phys. Rev. A 31 1695-7... 

Even expression ( B3.4.31), altiiough numerically preferable, is not the end of the story as it does not fiilly account for the fact diat nearby classical trajectories (those with similar initial conditions) should be averaged over. One simple methodology for that averaging has been tln-ough the division of phase space into parts, each of which is covered by a set of Gaussians [159, 160]. This is done by recasting the initial wavefunction as... 

Hamiltonian) trajectory in the phase space of the model from which infonnation about the equilibrium dyuamics cau readily be extracted. The application to uou-equilibrium pheuomeua (e.g., the kinetics of phase separation) is, in principle, straightforward. 

It is convenient to analyse tliese rate equations from a dynamical systems point of view similar to tliat used in classical mechanics where one follows tire trajectories of particles in phase space. For tire chemical rate law (C3.6.2) tire phase space , conventionally denoted by F, is -dimensional and tire chemical concentrations, CpC2,- are taken as ortliogonal coordinates of F, ratlier tlian tire particle positions and velocities used as tire coordinates in mechanics. In analogy to classical mechanical systems, as tire concentrations evolve in time tliey will trace out a trajectory in F. Since tire velocity functions in tire system of ODEs (C3.6.2) do not depend explicitly on time, a given initial condition in F will always produce tire same trajectory. The vector R of velocity functions in (C3.6.2) defines a phase-space (or trajectory) flow and in it is often convenient to tliink of tliese ODEs as describing tire motion of a fluid in F with velocity field/ (c p). 

F after transients have decayed. This final set of phase-space points is tire attractor, and tire set of all initial conditions tliat eventually reaches tire attractor is called its basin of attraction. 

For certain parameter values tliis chemical system can exlribit fixed point, periodic or chaotic attractors in tire tliree-dimensional concentration phase space. We consider tire parameter set... 

The description of chemical reactions as trajectories in phase space requires that the concentrations of all chemical species be measured as a function of time, something that is rarely done in reaction kinetics studies. In addition, the underlying set of reaction intennediates is often unknown and the number of these may be very large. Usually, experimental data on the time variation of the concentration of a single chemical species or a small number of species are collected. (Some experiments focus on the simultaneous measurement of the concentrations of many chemical species and correlations in such data can be used to deduce the chemical mechanism [7].)... 

The trajectory description problem of chemical reactions is resolved by using phase-space reconstmction from a single time series [8] this method uses delayed data at times t, t+ip t+X2,.. ., for an -dimensional attractor,... 

---