Gaussian wavepacket

Figure Al.6.1. Gaussian wavepacket in a hannonic oscillator. Note tliat the average position and momentum change according to the classical equations of motion (adapted from [6]).

It is not difficult to show that, for a constant potential, equation (A3.11.218) and equation (A3.11.219) can be solved to give the free particle wavepacket in equation (A3.11.7). More generally, one can solve equation (A3.11.218) and equation (A3.11.219) numerically for any potential, even potentials that are not quadratic, but the solution obtained will be exact only for potentials that are constant, linear or quadratic. The deviation between the exact and Gaussian wavepacket solutions for other potentials depends on how close they are to bemg locally quadratic, which means... 

Gaussian wavepackets are very special fiinctions which, in a sense, bridge the gap between classical and quantum... 

Figure B3.4.17. When a wavepacket comes to a crossing point, it will split into two parts (schematic Gaussians). One will remain on the same adiabat (difFerent diabat) and the other will hop to the other adiabat (same diabat). The adiabatic curves are shown by fidl lines and denoted by ground and excited die diabatic curves are shown by dashed lines and denoted 1, 2.

To remedy this diflSculty, several approaches have been developed. In some metliods, the phase of the wavefunction is specified after hopping [178]. In other approaches, one expands the nuclear wavefunction in temis of a limited number of basis-set fiinctions and works out the quantum dynamical probability for jumping. For example, the quantum dynamical basis fiinctions could be a set of Gaussian wavepackets which move forward in time [147]. This approach is very powerfLil for short and intemiediate time processes, where the number of required Gaussians is not too large. 

Figure C3.5.7. Possible modes of vibrational wavepacket (smootli Gaussian curve) motion for a highly vibrationally excited diatomic molecule produced by photodissociation of a linear triatomic such as Hglj, from [8].

In this seiniclassical calculation, we use only one wavepacket (the classical path limit), that is, a Gaussian wavepacket, rather than the general expansion of the total wave function. Equation (39) then takes the following form ... 

By substituting these expressions into Eq. (55), one can see after some algebra that ln,g(x, t) can be identified with lnx (t) + P t) shown in Section III.C.4. Moreover, In (f) = 0. It can be verified, numerically or algebraically, that the log-modulus and phase of In X-(t) obey the reciprocal relations (9) and (10). In more realistic cases (i.e., with several Gaussians), Eq. (56-58) do not hold. It still may be due that the analytical properties of the wavepacket remain valid and so do relations (9) and (10). If so, then these can be thought of as providing numerical checks on the accuracy of approximate wavepackets. 

A different approach is to represent the wavepacket by one or more Gaussian functions. When using a local harmonic approximation to the trae PES, that is, expanding the PES to second-order around the center of the function, the parameters for the Gaussians are found to evolve using classical equations of motion [22-26], Detailed reviews of Gaussian wavepacket methods are found in [27-29]. 

To add non-adiabatic effects to semiclassical methods, it is necessary to allow the trajectories to sample the different surfaces in a way that simulates the population transfer between electronic states. This sampling is most commonly done by using surface hopping techniques or Ehrenfest dynamics. Recent reviews of these methods are found in [30-32]. Gaussian wavepacket methods have also been extended to include non-adiabatic effects [33,34]. Of particular interest here is the spawning method of Martinez, Ben-Nun, and Levine [35,36], which has been used already in a number of direct dynamics studies. 

The Gaussian wavepacket based spawning method, mentioned above, has also been used in direct dynamics where it is called ab initio multiple spawning... 

Finally, Gaussian wavepacket methods are described in which the nuclear wavepacket is described by one or more Gaussian functions. Again the equations of motion to be solved have the fomi of classical trajectories in phase space. Now, however, each trajectory has a quantum character due to its spread in coordinate space. 

The fundamental method [22,24] represents a multidimensional nuclear wavepacket by a multivariate Gaussian with time-dependent width niaUix, A center position vector, R, momentum vector, p and phase, y,... 

The big advantage of the Gaussian wavepacket method over the swarm of trajectory approach is that a wave function is being used, which can be easily manipulated to obtain quantum mechanical information such as the spechum, or reaction cross-sections. The initial Gaussian wave packet is chosen so that it... 

One drawback is that, as a result of the time-dependent potential due to the LHA, the energy is not conserved. Approaches to correct for this approximation, which is valid when the Gaussian wavepacket is narrow with respect to the width of the potential, include that of Coalson and Karplus [149], who use a variational principle to derive the equations of motion. This results in replacing the function values and derivatives at the central point, V, V, and V" in Eq. (41), by values averaged over the wavepacket. 

The method will, however, fail badly if the Gaussian form is not a good approximation. For example, looking at the dynamics shown in Figure 2, a problem arises when a barrier causes the wavepacket to bifurcate. Under these circumstances it is necessary to use a superposition of functions. As will be seen later, this is always the case when non-adiabatic effects are present. 

Sawada et al. [26] made a detailed study of the methodology and numerical properties of the method. They paid paiticular attention to the problem of using a superposition of Gaussian wavepackets... 

The standard semiclassical methods are surface hopping and Ehrenfest dynamics (also known as the classical path (CP) method [197]), and they will be outlined below. More details and comparisons can be found in [30-32]. The multiple spawning method, based on Gaussian wavepacket propagation, is also outlined below. See [1] for further infomiation on both quantum and semiclassical non-adiabatic dynamics methods. 

Both the BO dynamics and Gaussian wavepacket methods described above in Section n separate the nuclear and electronic motion at the outset, and use the concept of potential energy surfaces. In what is generally known as the Ehrenfest dynamics method, the picture is still of semiclassical nuclei and quantum mechanical electrons, but in a fundamentally different approach the electronic wave function is propagated at the same time as the pseudoparticles. These are driven by standard classical equations of motion, with the force provided by an instantaneous potential energy function... 

A more comprehensive Gaussian wavepacket based method has been introduced by Martinez et al. [35,36,218]. Called the multiple spawning method, it has already been used in direct dynamics studies (see Section V.B), and shows much promise. It has also been applied to adiabatic problems in which tunneling plays a role [219], as well as the interaction of a... 

---