Central trajectories

In older instruments, a velocity filter placed before the magnetic field was used to eliminate v. This was done by combining two opposing forces created by an electric and magnetic field, which would only permit ions remaining on a central trajectory to exit the filter. 

For this, in older instruments, a device situated prior to the magnetic sector, called a velocity filter was employed. This was done by submitting the ions to two opposing forces through the action of both an electric E and a magnetic field B. Only the ions remaining on the central trajectory were able to exit the filter. Because the net resulting force must be zero, then qE = qvB and therefore v = E/B. 

The force exerted upon the ions, forced to travel tangentially in the electric sector, does not modify their kinetic energy, as the force is perpendicular to the imposed central trajectory that has a radius of curvature R (Figure 16.5). Its... 

Here represents the central trajectories of the k2 states in x as found from the undistorted model and fl" represents the central estimates of the pi... 

The cylinder leading to the transition state is very obvious in this figure. There are several points to observe first, at the bottom of the figure, where the r2 arc length is 0, is the transition state. All trajectories at the transition state lie within the boundary of the perimeter family at this point. Second, as the trajectories leave the transition state, they remain clustered as they were at the transition state, sweeping out a cylinder in the phase space. Third, the cylinder moves as it evolves, following a somewhat twisted path in the phase space, predictable from the path of the central trajectory. Viewed in terms of (r, p,), the trajectories of the perimeter family oscillate around the outside of the cylinder. As r2 propagates, the trajectories then wind around the outside of the cylinder. 

Figure 11 Phase-space surface px versus r, versus r2 arc length for the trajectories in the perimeter family (dotted lines) and for the central trajectory (solid line) of the two-mode hydrazoic acid model at a total energy of 0.0610 au. The reaction cylinder is conspicuous as these trajectories move through phase space.

In general, we cannot get a closed-form solution for the central trajectory, but can use the following algorithm. 

Trajectories of projectiles impinging on a target for reactions below the barrier left), where only elastic scattering occurs and above the barrier right), where the more central trajectories are absorbed. As is seen, nuclear reactions will necessarily decrease the large angle fraction of elastically scattered projectiles... 

Optic axis In the optical as opposed to the ballistic study of particle motion in electric and magnetic fields, the behavior of particles that remain in the neighborhood of a central trajectory is studied. This central trajectory is known as the optic axis. 

An electric deflection sector is normally made from two cylindrical electrodes, with the electric field applied between the two electrodes (Figure 13.9, top left). The force exerted on the charged particle does not modify its initial kinetic energy, because it is perpendicular to the central trajectory (with curvature rj, but it induces a deviation from this trajectory on changing the applied field E. Recalling that the centripetal component of the acceleration in the system is given by Oc — v /fc together with the -dependent... 

Figure Bl.7.9. (a) Stability diagram for ions near the central axis of a quadnipole mass filter. Stable trajectories occur only if the and values lie beneath tire curve, (b) Stability diagram (now as a fiinction of U and F) for six ions with different masses. The straight line miming tlirough the apex of each set of curves is the operating line, and conesponds to values of UIVthat will produce mass resolution (reproduced with pennission of Professor R March, Trent University, Peterborough, ON, Canada).

The basic assumption here is the existence over the inelastic scattering region of a connnon classical trajectory R(t) for the relative motion under an appropriately averaged central potential y[R(t)]. The interaction V r, / (t)] between A and B may then be considered as time-dependent. The system wavefiinction therefore satisfies... 

The top part of Fig. 1 shows the time evolution of the central dihedral angle of butane, r (defined by the four carbon atoms), for trajectories... 

Fig. 1. The time evolution (top) and average cumulative difference (bottom) associated with the central dihedral angle of butane r (defined by the four carbon atoms), for trajectories differing initially in 10 , 10 , and 10 Angstoms of the Cartesian coordinates from a reference trajectory. The leap-frog/Verlet scheme at the timestep At = 1 fs is used in all cases, with an all-atom model comprised of bond-stretch, bond-angle, dihedral-angle, van der Waals, and electrostatic components, a.s specified by the AMBER force field within the INSIGHT/Discover program.

Note that Equation 25.1 shows that the field (F) has no effect along the direction of the central (z) axis of the quadmpole assembly, so, to make ions move in this direction, they must first be accelerated through a small electric potential (typically 5 V) between the ion source and the assembly. Because of the oscillatory nature of the field (F Figure 25.3), an ion trajectory as it moves through the quadmpole assembly is also oscillatory. 

Direct dynamics trajectory calculations at the MP2/6-31-FG level of theory were then used to explore the reaction dynamics of this system [63]. Sixty-four trajectories were started from the central barrier shown at A in Fig. 11, with initial conditions sampled from a 300 K Boltzmann distribution. Of the 31 trajectories that moved in the direction of products, four trajectories followed the MEP and became trapped in the hydrogen-bonded [CH3OH ... 

It was necessary periodically to generate an adiabatic trajectory in order to obtain the odd work and the time correlation functions. In calculating E (t) on a trajectory, it is essential to integrate E)(t) over the trajectory rather than use the expression for E (T(f)) given earlier. This is because is insensitive to the periodic boundary conditions, whereas j depends on whether the coordinates of the atom are confined to the central cell, or whether the itinerant coordinate is used, and problems arise in both cases when the atom leaves the central cell on a trajectory. 

Figure 5. Molecular dynamics simulation of the decay forward and backward in time of the fluctuation of the first energy moment of a Lennard-Jones fluid (the central curve is the average moment, the enveloping curves are estimated standard error, and the lines are best fits). The starting positions of the adiabatic trajectories are obtained from Monte Carlo sampling of the static probability distribution, Eq. (246). The density is 0.80, the temperature is Tq — 2, and the initial imposed thermal gradient is pj — 0.02. (From Ref. 2.)...

This chapter is organized as follows. In Section II, we briefly summarize the findings of the geometric TST for autonomous Hamiltonian systems to the extent that it is needed for the present discussion. Readers interested in a more detailed exposition are referred to Ref. 35, where the field has recently been reviewed in depth. We restrict our discussion to classical mechanics. Semiclassical extensions of geometric TST have been developed in Refs. 70-75. Section III discusses the notion of the TS trajectory in general and its incarnation in different specific settings. Section IV demonstrates how the TS trajectory allows one to carry over the central concepts of geometric TST into the time-dependent realm. 

The point q = p = 0 (or P = Q = 0) is a fixed point of the dynamics in the reactive mode. In the full-dimensional dynamics, it corresponds to all trajectories in which only the motion in the bath modes is excited. These trajectories are characterized by the property that they remain confined to the neighborhood of the saddle point for all time. They correspond to a bound state in the continuum, and thus to the transition state in the sense of Ref. 20. Because it is described by the two independent conditions q = 0 and p = 0, the set of all initial conditions that give rise to trajectories in the transition state forms a manifold of dimension 2/V — 2 in the full 2/V-dimensional phase space. It is called the central manifold of the saddle point. The central manifold is subdivided into level sets of the Hamiltonian in Eq. (5), each of which has dimension 2N — 1. These energy shells are normally hyperbolic invariant manifolds (NHIM) of the dynamical system [88]. Following Ref. 34, we use the term NHIM to refer to these objects. In the special case of the two-dimensional system, every NHIM has dimension one. It reduces to a periodic orbit and reproduces the well-known PODS [20-22]. 

Pi exponentially decreases to zero as t —> oo, whereas Q increases. Trajectories with <2i 0 therefore asymptotically approach the central manifold... 

Pi = Qi = 0 as t —y oo, or, more precisely, they approach the NHIM of the appropriate energy within the central manifold. Due to this behavior, the set of all initial conditions with Q = 0 is called the stable manifold of the NHIM. Similarly, trajectories with Pi = 0 asymptotically approach the NHIM as t > —oo. They are said to form the unstable manifold of the NHIM. 

In addition to describing the TS and the separatrices between reactive and nonreactive trajectories that are the central discovery of geometric TST, one can... 

This is exactly the autonomous linearized Hamiltonian (7), the dynamics of which was discussed in detail in Section II. One therefore finds the TS dividing surface and the full set of invariant manifolds described earlier one-dimensional stable and unstable manifolds corresponding to the dynamics of the variables A<2i and APt, respectively, and a central manifold of dimension 2N — 2 that itself decomposes into two-dimensional invariant subspaces spanned by APj and AQj. However, all these manifolds are now moving manifolds that are attached to the TS trajectory. Their actual location in phase space at any given time is obtained from their description in terms of relative coordinates by the time-dependent shift of origin, Eq. (42). 

Following Ref. 40, we study trajectories that approach the central region from the (arbitrarily chosen) reactant channel. They will be identified by their... 

For highly exothermic SN2 reactions, which have a central barrier significantly lower in energy than that of the reactants, association of the reactants may be the rate controlling step in TST.1 The SN2 rate constant can then be modeled by a capture theory9 such as VTST,10 average dipole orientation (ADO) theory,11 the statistical adiabatic channel model (SACM),12 or the trajectory capture model.13... 

Figure 4. Probability of different trajectory events versus the Cla-C-Clb angle 0, which is evaluated at the first inner turning point (ITP) for complex formation and at the central barrier for the trajectories which attain this configuration ( ), association probability (o), probability of attaining the central barrier. fre = 0.5 kcal/mol and nc-cib = 6. Trot = 0 K in (a) and 300 K in (b) (from ref. 38).

A dynamical model for SN2 nucleophilic substitution that emerges from the trajectory simulations is depicted in Figure 9. The complex formed by a collision between the reactants is an intermolecular complex CinterR. To cross the central barrier, this complex has to undergo a unimolecular transition in which energy is... 

The dynamical model described in Figure 9 indicates that the trajectories may recross the central barrier several times if the Cintra R Cintra p transition is faster... 

---