Mode modelling

The first classical trajectory study of iinimoleciilar decomposition and intramolecular motion for realistic anhannonic molecular Hamiltonians was perfonned by Bunker [12,13], Both intrinsic RRKM and non-RRKM dynamics was observed in these studies. Since this pioneering work, there have been numerous additional studies [9,k7,30,M,M, ai d from which two distinct types of intramolecular motion, chaotic and quasiperiodic [14], have been identified. Both are depicted in figure A3,12,7. Chaotic vibrational motion is not regular as predicted by tire nonnal-mode model and, instead, there is energy transfer between the modes. If all the modes of the molecule participate in the chaotic motion and energy flow is sufficiently rapid, an initial microcanonical ensemble is maintained as the molecule dissociates and RRKM behaviour is observed [9], For non-random excitation initial apparent non-RRKM behaviour is observed, but at longer times a microcanonical ensemble of states is fonned and the probability of decomposition becomes that of RRKM theory. 

Wyatt R E, lung C and Leforestier C 1992 Quantum dynamics of overtone relaxation in benzene. II. Sixteen-mode model for relaxation from CH(v = 3) J. Chem. Phys. 97 3477-86... 

It should be realized, though, that either model can be used for levels for all values of v for such stretching vibrations but the normal mode model is more practically useful at low v and the local mode model more useful at higher values of v. 

For example, in the case of H tunneling in an asymmetric 0i-H - 02 fragment the O1-O2 vibrations reduce the tunneling distance from 0.8-1.2 A to 0.4-0.7 A, and the tunneling probability increases by several orders. The expression (2.77a) is equally valid for the displacement of a harmonic oscillator and for an arbitrary Gaussian random value q. In a solid the intermolecular displacement may be contributed by various lattice motions, and the above two-mode model may not work, but once q is Gaussian, eq. (2.77a) will still hold, however complex the intermolecular motion be. 

The two-mode model has two characteristic cross-over temperatures corresponding with the freezing of each vibration. Above = hcoo/2k the dependence k(T) is Arrhenius, with activation energy equal to... 

The most important aspect of coralyne is its ability to inhibit DNA relaxation in a fashion significantly similar to the most potent antitumour alkaloid camptothecin, which is known to exert this property [242], Presumably, the most notable biological action of these alkaloids appears to be topoisomerase inhibition [238-242], which has direct relevance to their DNA intercalating property. In this context. Pilch et al. [167] described a mixed binding mode model (Fig. 16) in which the protoberberine structure constitutes portions that can intercalate or bind to the minor groove of DNA. Wang et al. [240] demonstrated that coralyne (Ci) and several of its derivatives (Ce to Ch) (Scheme 5), including the partial saturated... 

While principal components models are used mostly in an unsupervised or exploratory mode, models based on canonical variates are often applied in a supervisory way for the prediction of biological activities from chemical, physicochemical or other biological parameters. In this section we discuss briefly the methods of linear discriminant analysis (LDA) and canonical correlation analysis (CCA). Although there has been an early awareness of these methods in QSAR [7,50], they have not been widely accepted. More recently they have been superseded by the successful introduction of partial least squares analysis (PLS) in QSAR. Nevertheless, the early pattern recognition techniques have prepared the minds for the introduction of modem chemometric approaches. 

Starting from the normal mode approximation, one can introduce anharmonicity in different ways. Anharmonic perturbation theory [206] and local mode models [204] may be useful in some cases, where anharmonic effects are small or mostly diagonal. Vibrational self-consistent-field and configuration-interaction treatments [207, 208] can also be powerful and offer a hierarchy of approximation levels. Even more rigorous multidimensional treatments include variational calculations [209], diffusion quantum Monte Carlo, and time-dependent Hartree approaches [210]. 

Fleming, P. R., and Hutchinson J. S. (1988), Representation of the Hamiltonian Matrix in Non-Local Coordinates for an Acetylene Bond-Mode Model, Comp. Phys. Comm. 51, 59. 

The nuclear tunnelling factor can be accurately estimated from a 1-mode model based on the high frequency inner-sphere breathing mode (10, Tl)... 

To illustrate some of the characteristics of the two-mode model, we calculate the intervalence band of the Creutz and... 

We now show that the two-mode model can account for the observed intervalence band asymmetry of the C T ion, the weak tunneling transitions and a Ar value of 0.04 8. We fix the... 

The solid dots in Figure 5 are obtained from our two-mode model by diagonalizing a 132 x 132 matrix as discussed earlier and, hence, should be reasonably accurate over the entire range. In this calculation, the parameters kc and have been held... 

Figure 5. Mean band energy (E) vs. total reorganization energy [Eq. (18)]. The dashed curve is the equation, E = Et. The solid curves are calculated from an analytical expression (9). The solid dots are calculated from a diagonalization of the two-mode model secular determinant with c = —6.0, Ac, = 1.1, and vc, =...

If we start with the two-mode model, the problem can be formulated as follows. If the system starts in one of the minima of the lower potential surface (W of eq 12), what is the... 

Eq 17 is much more reasonable in the two-mode model when k is... 

In order to discuss various aspects of a mixed quantum-classical treatment of photoinduced nonadiabatic dynamics, we consider five different kinds of molecular models, each representing a specihc challenge for a mixed quantum-classical modeling. Here, we introduce the specifics of these models and discuss the characteristics of their nonadiabatic dynamics. The molecular parameters of the few-mode models (Model I-IV) describing intramolecular nonadiabatic dynamics are collected in Tables I-V. The parameters of Model V describing various aspects of nonadiabatic dynamics in the condensed phase will be given in the text. 

Parameters of Model I, Which Represents a Three-Mode Model of the Si — S2 Conical Intersection in Pyrazine [173] ... 

Parameters of Model II, Which Represents a Three-State Eive-Mode Model of the Ultrafast C — B — X Internal-Conversion Process in the Benzene Cation [179, 180] ... 

Parameters of Model IVa, Which Represents a One-Mode Model of Intramolecular... 

Figure 2. Diabatic (left) and adiabatic (right) population probabiUties of the C (fuU line), B (dotted line), and X (dashed line) electronic states as obtained for Model II, which represents a three-state five-mode model of the benzene cation. Shown are (A) exact quantum calculations of Ref. 180, as well as mean-field-trajectory results [(B), (E)] and surface-hopping results [(C),(D),(F),(G)]. The latter are obtained either directly from the electronic coefficients [(C),(F)] or from binned coefficients [(D),(G)].

Finally, we consider the performance of the MFT method for nonadiabatic dynamics induced by avoided crossings of the respective potential energy surfaces. We start with the discussion of the one-mode model. Model IVa, describing ultrafast intramolecular electron transfer. The comparison of the MFT method (dashed line) with the quantum-mechanical results (full line) shown in Fig. 5 demonstrates that the MFT method gives a rather good description of the short-time dynamics (up to 50 fs) for this model. For longer times, however, the dynamics is reproduced only qualitatively. Also shown is the time evolution of the diabatic electronic coherence which, too, is... 

As a first example, we again consider Model 1 describing a two-state three-mode model of the Si nn ) and 52(7171 ) states of pyrazine. Figure 11a shows the quantum-mechanical (thick line) and the SH (thin lines) results for the adiabatic population probability of the initially prepared electronic state /2). As... 

Finally, we consider Model V by describing two examples of outer-sphere electron-transfer in solution. Figures 7 and 8 display results for the diabatic electronic population for Models Va and Vb, respectively. Similar to the mean-field trajectory calculations, for Model Va the SH results are in excellent agreement with the quantum calculations, while for Model Vb the SH method only is able to describe the short-time dynamics. As for the three-mode Model IVb discussed above, the SH calculations in particular predict an incorrect long-time limit for the diabatic population. The origin of this problem will be discussed in more detail in Section VI in the context of the mapping formulation. 

Figure 25. Diabatic and adiabatic population probabilities of the C (fuU line), B (dotted hne), and X (dashed line) electronic states as obtained for a five-state 16-mode model of the benzene cation.

Because the mapping approach treats electronic and nuclear dynamics on the same dynamical footing, its classical limit can be employed to study the phase-space properties of a nonadiabatic system. With this end in mind, we adopt a onemode two-state spin-boson system (Model IVa), which is mapped on a classical system with two degrees of freedom (DoF). Studying various Poincare surfaces of section, a detailed phase-space analysis of the problem is given, showing that the model exhibits mixed classical dynamics [123]. Furthermore, a number of periodic orbits (i.e., solutions of the classical equation of motion that return to their initial conditions) of the nonadiabatic system are identified and discussed [125]. It is shown that these vibronic periodic orbits can be used to analyze the nonadiabatic quantum dynamics [126]. Finally, a three-mode model of nonadiabatic photoisomerization (Model III) is employed to demonstrate the applicability of the concept of vibronic periodic orbits to multidimensional dynamics [127]. 

Figure 32. Vibronic periodic orbits of a coupled electronic two-state system with a single vibrational mode (Model IVa). All orbits are displayed as a function of the nuclear position x and the electronic population N, where N = Aidia (left) and N = (right), respectively. As a further illustration, the three shortest orbits have been drawn as curves in between the diabatic potentials Vi and V2 (left) as well as in between the corresponding adiabatic potentials Wi and W2 (right). The shaded Gaussians schematically indicate that orbits A and C are responsible for the short-time dynamics following impulsive excitation of V2 at (xo,po) = (3,0), while orbit B and its symmetric partner determine the short-time dynamics after excitation of Vi at (xo,po) = (3, —2.45).

While the simulations of the pyrazine system discussed in the previous sections of this chapter have employed a three-mode model (Model I), the semiclassical simulations we will present here are based on two different models a four-mode model and a model including all 24 normal modes of the pyrazine molecule. Let us first consider the four-mode model of the S1-S2 conical intersection in pyrazine which was developed by Domcke and coworkers [269]. In addition to the three modes considered in Model I, it takes into account another Condon-active mode (V9a). Figure 37 shows the modulus of the autocorrelation function [cf. Eq. (24)] of this model after photoexcitation to the S2 electronic state. The exact quantum results (full line) are compared to the... 

Figure 37. Modulus of the autocorrelation function for the four-mode model of pyrazine. The full line is the quantum result, and the dotted line is the semiclassical result.

Figure 38. Absorption spectrum of pyrazine in the energy region of the S1-S2 conical intersection. Shown are (a) quantum mechanical (full line) and semiclassical (dotted line) results for the four-mode model (including a phenomenological dephasing constant of T2 = 30 fs) and (b) the experimental data [271],...

---