Classical trajectories direct reactions

In the context of reaction path studies. Billing has employed classical trajectories directly to simulate the dynamics of an approximation to the reaction path Hamiltonian (including rotational motion) [89,93]. Classi-... 

The Restart check box can be used in ctiii junction with the explicit editing of a IIIX file to assign completely user-specified initial velocities. This may be useful in classical trajectory analysis of chemical reactions where the initial velocities and directions of the reactants are varied to statistically determine the probability of reaction occurring, or n ot, in the process of calculating a rate con -Stan t. 

The second channel, producing CO, was first observed by Seakins and Leone [64], who estimated 40% branching to this channel. Later measurements by Lockenberg et al. [65] and Preses et al. [66] concluded the branching to CO is 18%. Note that decomposition of formaldehyde formed in reaction (26a) is not a possible source of CO due to the large barrier for formaldehyde decomposition. Marcy et al. [67] recently combined time-resolved Lourier spectroscopy experiments with direct dynamics classical trajectory calculations to examine the mechanism of the CO product channel. They observed two pathways for CO formation, neither of which involve crossing a TS. 

The interpretation of a spectrum from a dynamical point of view can also be applied to a spectrum containing a broad feature associated with direct and/or indirect dissociation reactions. From such spectra dynamics of a dissociating molecule can also be extracted via the Fourier transform of a spectrum. An application of the Fourier transform to the Hartley band of ozone by Johnson and Kinsey [3] demonstrated that a small oscillatory modulation built on a broad absorption feature contains information of the classical trajectories of the vibrational motion on PES, so-called unstable periodic orbits, at the transition state of a unimolecular dissociation. 

A perturbation-trajectory method for determining the dynamics of gas-surface collision processes was tested on the collision and subsequent surface reactions of SiH2 on a Si(lll) surface233. The predictions of an exact classical trajectory calculation234 were confirmed the sticking probabilities were unity at all temperatures, and it was found that surface SiH2 can decompose by direct elimination of H2 or by successive dissociation of Si—H bonds. 

As seen previously, the chemical reactions studied most often are the exchange ones. Those requiring several potential energy surfaces of excited states (diabatic reactions) are worth special mention, since they most certainly define a domain of application with a future for classical trajectories. An electron jump from one surface to another requires either to be given a statistical probability of occurence by the Landau Zener formula (or one of its improved versions " ) or to be described by means of complex-valued classical trajectories as a direct and gradual passage in the complex-valued extension of the potential surfaces (generalization of the classical S-matrix ). 

The key idea that supplements RRK theory is the transition state assumption. The transition state is assumed to be a point of no return. In other words, any trajectory that passes through the transition state in the forward direction will proceed to products without recrossing in the reverse direction. This assumption permits the identification of the reaction rate with the rate at which classical trajectories pass through the transition state. In combination with the ergodic approximation this means that the reaction rate coefficient can be calculated from the rate at which trajectories, sampled from a microcanonical ensemble in the reactants, cross the barrier, divided by the total number of states in the ensemble at the required energy. This quantity is conveniently formulated using the idea of phase space. 

Another advantage of the quantum calculations is that they provide a rigorous test of approximate methods for calculating dissociation rates, namely classical trajectories and statistical models. Two commonly used statistical theories are the Rice-Ramsperger-Kassel-Marcus (RRKM) theory and the statistical adiabatic channel model (SACM). The first one is thoroughly discussed in Chapter 2, while the second one is briefly reviewed in the Introduction. Moreover, the quantum mechanical approach is indispensable in analyzing the reaction mechanisms. A resonance state is characterized not only by its position, width and the distribution of product states, but also by an individual wave function. Analysis of the nodal structure of resonance wave functions gives direct access to the mechanisms of state- and mode-selectivity. 

Here again, we take a classical trajectory as a reference path in a reaction tube that passes across the transition region between two basins a and b with the flow direction b —> a. Set the time origin t = 0 at just the moment of transition. At a given time t, we take a sphere of a radius rt in 30-dimensional phase space, the center of which proceeds along the reference trajectory. Pick random points in this sphere, and let them run backward in time. Some of them will go back to the basin b if the sphere still lies inside the same tube, and the others will move to some other basins if this sphere is already out of the tube. Should the latter happen, a similar procedure is to be redone with a smaller radius r,. Repeating... 

In addition to the IRC functionality, gamess also has a direct dynamics capability, the dynamic reaction path, DRC [59]. The DRC allows one to perform dynamics on-the-fly , by performing classical trajectories at any level of theory for which analytic gradients are available. One can, for example, put an amount of energy equal to n quanta into any vibrational mode(s), in order to model mode specific chemistry. 

Generally, we expect that a classical trajectory or quantum dynamical calculation (say, by discrete variable methods) of the reaction rate for a tetraatomic system would require at least 10 -10 evaluations of the energy (or energy and gradient). This selective iterative interpolation procedure appears to produce an accurate PES with about three orders of magnitude less ab initio computation than that of the best direct method [208]. 

Reaction path methods have great promise for future progress because they can be made direct or automated. Direct ab initio dynamics is taken to mean the ab initio evaluation of the molecular energy whenever it is required in a dynamical calculation. At present this is prohibitively expensive. A more traditional approach uses electronic structure methods to provide information about a PES which is then fitted to a functional form. The dynamics employed thereafter is restricted only by the limitations of current dynamical theories. However, the process of fitting functional forms to a molecular PES is difficult, unsystematic, and extremely time consuming. The reaction path approaches that use interpolation of the reaction path data are affordable methods aimed toward direct dynamics, as they avoid the process of fitting a functional form for the PES. Such methods have been automated (programmed) and are therefore readily usable. Thus reaction path methods have been applied with various levels of ab initio theory to statistical theories of the reaction rate, to approximate quantum dynamics, and to classical trajectory studies of reactions. 

Jaffe and Anderson explored the eno gy requirements for reaction ( -88b) by looking for the products of reaction when acoeteated beams of HI were directed into a long chamber containing DI at low pressure. The mean collision i gy was thereby varied in the range 84—456 kJ mol" but in no case could HD be detected in the products. This result indicated that for reaction to proceed some vibrational - or possibly rotational - excitation of the reactants is necessary. These reactions have also been the subject of theoretical studies - using classical trajectories, although the conclusions to be drawn from these calculations are a matter of some controversy. "... 

Lack of space prevents us from describing in any detail other properties of pods. It should be mentioned though that pods may be characterized as repulsive or attractive according to the local behaviour of trajectories in their vicinity. Loosely speaking, repulsive pods repel trajectories in their vicinity, attractive pods attract them. These properties have enabled the development of a lower bound to the classical reaction probability, the construction of a new theory of direct reactions as well as new numerical methods for computation of classical product state distributions.The repulsive attractive properties provide a global picture of the classical flow of trajectories at a given energy. 

---