Quantum corrections dynamical

Classical and semiclassical moment expressions. The expressions for the spectral moments can be made classical by substituting the classical distribution function, g(R) = exp (— V(R)/kT), for the quantum expressions. Wigner-Kirkwood corrections are known which account to lowest order for the static quantum corrections, Eq. 5.44 [177, 292]. For the second and higher moments, dynamic quantum corrections must also be made [177]. As was mentioned in the previous Chapter, such semiclassical corrections are useful in supplementing quantum computations of the spectral moments at large separations where the quantum effects are small the computational effort of quantum calculations, which is substantial at large separations, may thus be avoided. 

Equation (C3.5.3) shows tire VER lifetime can be detennined if tire quantum mechanical force-correlation Emotion is computed. However, it is at present impossible to compute tliis Emotion accurately for complex systems. It is straightforward to compute tire classical force-correlation Emotion using classical molecular dynamics (MD) simulations. Witli tire classical force-correlation function, a quantum correction factor Q is needed 5,... 

Behrens, P.H., Mackay, D.H.J., White, G.M., Wilson, K.R. Thermodynamics and quantum corrections from molecular dynamics for liquid water. J. Chem. Phys. 79 (1983) 2375-2389. 

Fig. 7.12 Experimental and calculated infrared spectra for liquid water. The black dots are the experimental values. The thick curve is the classical profile produced by the molecular dynamics simulation. The thin curve is obtained by applying quantum corrections. (Figure redrawn from Guilbt B 1991. A Molecular Dynamics Study of the Infrared Spectrum of Water. Journal of Chemical Physics 95 1543-1551.)...

The approach to the evaluation of vibrational spectra described above is based on classical simulations for which quantum corrections are possible. The incorporation of quantum effects directly in simulations of large molecular systems is one of the most challenging areas in theoretical chemistry today. The development of quantum simulation methods is particularly important in the area of molecular spectroscopy for which quantum effects can be important and where the goal is to use simulations to help understand the structural and dynamical origins of changes in spectral lineshapes with environmental variables such as the temperature. The direct evaluation of quantum time- correlation functions for anharmonic systems is extremely difficult. Our initial approach to the evaluation of finite temperature anharmonic effects on vibrational lineshapes is derived from the fact that the moments of the vibrational lineshape spectrum can be expressed as functions of expectation values of positional and momentum operators. These expectation values can be evaluated using extremely efficient quantum Monte-Carlo techniques. The main points are summarized below. 

After an overview of the main papers devoted to chaos in lasers (Section I.A) and in nonlinear optical processes (Section I.B), we present a more detailed analysis of dynamics in a process of second-harmonic generation of light (Section II) as well as in Kerr oscillators (Section III). The last case we consider particularly in the context of coupled nonlinear systems. Finally, we present a cumulant approach to the problem of quantum corrections to the classical dynamics in second-harmonic generation and Kerr processes (Section IV). 

This is the semi-classical, second, binary moment. While the zeroth and first moments require only static quantum corrections of the Wigner-Kirkwood type, the second and all higher moments require also dynamical corrections involving yC4 and higher moments. 

Induced dipole autocorrelation functions of three-body systems have not yet been computed from first principles. Such work involves the solution of Schrodinger s equation of three interacting atoms. However, classical and semi-classical methods, especially molecular dynamics calculations, exist which offer some insight into three-body dynamics and interactions. Very useful expressions exist for the three-body spectral moments, with the lowest-order Wigner-Kirkwood quantum corrections which were discussed above. 

Higher-order classical moments have also been reported. We mention the classical expressions for the translational spectral moments M , with n = 0, 2, 4, and 6, for pairs of linear molecules given in an appendix of [204]. Spectral moments of spherical top molecules have been similarly considered [163, 205], We note that for n > 1, spectral moments show dynamic as well as static quantum correction, which become more important as the order n of the spectral moments is increased. The discussions on pp. 219, and Table 5.1, suggest that, even for the near-classical systems, quantum corrections may be substantial and can rarely be ignored. 

Because of equipartition all frequencies should be equally weighted (i.e. fV((o) = 1 for all ). This clearly contradicts our understanding of dynamics from quantum mechanics (QM) theory. Since QM is essential to the treatment of certain phenomena, quantum corrections must be introduced, regardless of whether structural or dynamic simulations are required. This discrepancy is illustrated if we calculate the total kinetic energy of the system quantum-mechanically ... 

For temperatures below the Debye temperature (9d), quantum corrections must be applied to the temperature and thermal conductivity obtained from molecular dynamics. These quantum corrections are negligible for T 0, where the system behaves classically. A quantum correction for the temperature can be estimated by equating the ensemble s total energy to the phonons total energy [10, 57, 61] as ... 

Figure 2. Quantum corrected temperature (right axis), and ratio of quantum corrected to molecular dynamics predicted thermal conductivity k/kMo) (left axis), as a function of the temperature of the MD simulations for silicon.

Computer simulation of molecular dynamics is concerned with solving numerically the simultaneous equations of motion for a few hundred atoms or molecules that interact via specified potentials. One thus obtains the coordinates and velocities of the ensemble as a function of time that describe the structure and correlations of the sample. If a model of the induced polarizabilities is adopted, the spectral lineshapes can be obtained, often with certain quantum corrections [425,426]. One primary concern is, of course, to account as accurately as possible for the pairwise interactions so that by carefully comparing the calculated with the measured band shapes, new information concerning the effects of irreducible contributions of inter-molecular potential and cluster polarizabilities can be identified eventually. Pioneering work has pointed out significant effects of irreducible long-range forces of the Axilrod-Teller triple-dipole type [10]. Very recently, on the basis of combined computer simulation and experimental CILS studies, claims have been made that irreducible three-body contributions are observable, for example, in dense krypton [221]. 

F. Barocchi, M. Zoppi, and M. Neumann. First-order quantum corrections to depolarized interaction induced light scattering spectral moments Molecular dynamics calculation. Phys. Rev. A, 27 1587-1593 (1983). 

Inclusion of dynamical effects allows calculation of corrections to simple Transition State Theory, often described by a transmission coefficient k to be multiplied with the TST rate constant (Section 12.1), or used in connection with variational TST (Section 12.3). Classical dynamics allow corrections due to recrossing to be calculated, while a quantum treatment is necessary for including tunnelling effects. Owing to the stringent... 

Dynamical Quantum Effects. When dynamical quantum effects become important, corrections must be made. At the simplest level, quantum corrections can be made following Hirschfelder et al 51 allowing corrections for... 

Region I The dynamics is over the barrier (as described by TST) or just below the barrier (small quantum corrections). 

It is important to clarify here that the description of PT processes by curve crossing formulations is not a new approach nor does it provide new dynamical insight. That is, the view of PT in solutions and proteins as a curve crossing process has been formulated in early realistic simulation studies [1, 2, 42] with and without quantum corrections and the phenomenological formulation of such models has already been introduced even earlier by Kuznetsov and others [47]. Furthermore, the fact that the fluctuations of the environment in enzymes and solution modulate the activation barriers of PT reactions has been demonstrated in realistic microscopic simulations of Warshel and coworkers [1, 2]. However, as clarified in these works, the time dependence of these fluctuations does not provide a useful way to determine the rate constant. That is, the electrostatic fluctuations of the environment are determined by the corresponding Boltzmann probability and do not represent a dynamical effect. In other words, the rate constant is determined by the inverse of the time it takes the system to produce a reactive trajectory, multiplied by the time it takes such trajectories to move to the TS. The time needed for generation of a reactive trajectory is determined by the corresponding Boltzmann probability, and the actual time it takes the reactive trajectory to reach the transition state (of the order of picoseconds), is more or less constant in different systems. 

---