Nuclear dynamics

Vos M H, Jones M R, Hunter C N, Breton J and Martin J-L 1994 Coherent nuclear dynamics at room temperature in bacterial reaction centers Proc. Natl Acad. Sci. USA 91 12 701-5... 

The discussion in the previous sections assumed that the electron dynamics is adiabatic, i.e. the electronic wavefiinction follows the nuclear dynamics and at every nuclear configuration only the lowest energy (or more generally, for excited states, a single) electronic wavefiinction is relevant. This is the Bom-Oppenlieimer approxunation which allows the separation of nuclear and electronic coordinates in the Schrodinger equation. 

Diz A, Deumens E and Ohm Y 1990 Quantum electron-nuclear dynamics Chem. Phys. Lett. 166 203... 

We now consider the connection between the preceding equations and the theory of Aharonov et al. [18] [see Eqs. (51)-(60)]. The tempting similarity between the structures of Eqs. (56) and (90), hides a fundamental difference in the roles of the vector operator A in Eq. (56) and the vector potential a in Eq. (90). The fomrer is defined, in the adiabatic partitioning scheme, as a stiictly off-diagonal operator, with elements (m A n) = (m P n), thereby ensuring that (P — A) is diagonal. By contiast, the Mead-Truhlar vector potential a arises from the influence of nonzero diagonal elements, (n P /i) on the nuclear equation for v), an aspect of the problem not addressed by Arahonov et al. [18]. Suppose, however, that Eq. (56) was contracted between (n and n) v) in order to handle the adiabatic nuclear dynamics within the Aharonov scheme. The result becomes... 

Some final comments on the relevance of non-adiabatic coupling matrix elements to the nature of the vector potential a are in order. The above analysis of the implications of the Aharonov coupling scheme for the single-surface nuclear dynamics shows that the off-diagonal operator A provides nonzero contiibutions only via the term (n A n). There are therefore no necessary contributions to a from the non-adiabatic coupling. However, as discussed earlier, in Section IV [see Eqs. (34)-(36)] in the context of the x e Jahn-Teller model, the phase choice t / = —4>/2 coupled with the identity... 

V is the derivative with respect to R.) We stress that in this formalism, I and J denote the complete adiabatic electronic state, and not a component thereof. Both /) and y) contain the nuclear coordinates, designated by R, as parameters. The above line integral was used and elaborated in calculations of nuclear dynamics on potential surfaces by several authors [273,283,288-301]. (For an extended discussion of this and related matters the reviews of Sidis [48] and Pacher et al. [49] are especially infonnative.)... 

C, General Electron Nuclear Dynamics TV. Molecular Processes... 

Obviously, the BO or the adiabatic states only serve as a basis, albeit a useful basis if they are determined accurately, for such evolving states, and one may ask whether another, less costly, basis could be Just as useful. The electron nuclear dynamics (END) theory [1-4] treats the simultaneous dynamics of electrons and nuclei and may be characterized as a time-dependent, fully nonadiabatic approach to direct dynamics. The END equations that approximate the time-dependent Schrddinger equation are derived by employing the time-dependent variational principle (TDVP). 

A major drawback to the approaches that use preconstructed PESs is that there are many more interesting systems undergoing reactive dynamical processes than there are available PESs. It would be much better if the electronic structure part of the problem, which provides the forces for the nuclear dynamics, could be performed simultaneously with the dynamics part, thus making possible the treatment of systems for which preconstructed PESs... 

Election nuclear dynamics theory is a direct nonadiababc dynamics approach to molecular processes and uses an electi onic basis of atomic orbitals attached to dynamical centers, whose positions and momenta are dynamical variables. Although computationally intensive, this approach is general and has a systematic hierarchy of approximations when applied in an ab initio fashion. It can also be applied with semiempirical treatment of electronic degrees of freedom [4]. It is important to recognize that the reactants in this approach are not forced to follow a certain reaction path but for a given set of initial conditions the entire system evolves in time in a completely dynamical manner dictated by the inteiparbcle interactions. 

Knowledge of the underlying nuclear dynamics is essential for the classification and description of photochemical processes. For the study of complicated systems, molecular dynamics (MD) simulations are an essential tool, providing information on the channels open for decay or relaxation, the relative populations of these channels, and the timescales of system evolution. Simulations are particularly important in cases where the Bom-Oppenheimer (BO) approximation breaks down, and a system is able to evolve non-adiabatically, that is, in more than one electronic state. 

In this section, the basic theory of molecular dynamics is presented. Starting from the BO approximation to the nuclear Schrddinger equation, the picture of nuclear dynamics is that of an evolving wavepacket. As this picture may be unusual to readers used to thinking about nuclei as classical particles, a few prototypical examples are shown. 

Full quantum wavepacket studies on large molecules are impossible. This is not only due to the scaling of the method (exponential with the number of degrees of freedom), but also due to the difficulties of obtaining accurate functions of the coupled PES, which are required as analytic functions. Direct dynamics studies of photochemical systems bypass this latter problem by calculating the PES on-the-fly as it is required, and only where it is required. This is an exciting new field, which requires a synthesis of two existing branches of theoretical chemistry—electronic structure theory (quantum chemistiy) and mixed nuclear dynamics methods (quantum-semiclassical). 

Quantum chemical methods, exemplified by CASSCF and other MCSCF methods, have now evolved to an extent where it is possible to routinely treat accurately the excited electronic states of molecules containing a number of atoms. Mixed nuclear dynamics, such as swarm of trajectory based surface hopping or Ehrenfest dynamics, or the Gaussian wavepacket based multiple spawning method, use an approximate representation of the nuclear wavepacket based on classical trajectories. They are thus able to use the infoiination from quantum chemistry calculations required for the propagation of the nuclei in the form of forces. These methods seem able to reproduce, at least qualitatively, the dynamics of non-adiabatic systems. Test calculations have now been run using duect dynamics, and these show that even a small number of trajectories is able to produce useful mechanistic infomiation about the photochemistry of a system. In some cases it is even possible to extract some quantitative information. 

In this section, the spin-orbit interaction is treated in the Breit-Pauli [13,24—26] approximation and incoi porated into the Hamiltonian using quasidegenerate perturbation theory [27]. This approach, which is described in [8], is commonly used in nuclear dynamics and is adequate for molecules containing only atoms with atomic numbers no larger than that of Kr. 

The full quantum mechanical study of nuclear dynamics in molecules has received considerable attention in recent years. An important example of such developments is the work carried out on the prototypical systems H3 [1-5] and its isotopic variant HD2 [5-8], Li3 [9-12], Na3 [13,14], and HO2 [15-18], In particular, for the alkali metal trimers, the possibility of a conical intersection between the two lowest doublet potential energy surfaces introduces a complication that makes their theoretical study fairly challenging. Thus, alkali metal trimers have recently emerged as ideal systems to study molecular vibronic dynamics, especially the so-called geometric phase (GP) effect [13,19,20] (often referred to as the molecular Aharonov-Bohm effect [19] or Berry s phase effect [21]) for further discussion on this topic see [22-25], and references cited therein. The same features also turn out to be present in the case of HO2, and their exact treatment assumes even further complexity [18],... 

In this chapter, we discussed the permutational symmetry properties of the total molecular wave function and its various components under the exchange of identical particles. We started by noting that most nuclear dynamics treatments carried out so far neglect the interactions between the nuclear spin and the other nuclear and electronic degrees of freedom in the system Hamiltonian. Due to... 

---