A semiclasical justification for the use of non-spreading wavepackets in dynamics calculations

A justification is given for the use of non-spreading or frozen gaussian packets in dynamics calculations. In this work an initial wavefunction or quantum density operator is expanded in a complete set of grussian wavepackets. It is demonstrated that the time evolution of this wavepacket expansion for the quantum wavefunction or density is correctly given within the approximations employed by the classical propagation of the avarage position and momentum of each gaussian packet, holding the shape of these individual gaussians fixed. The semiclassical approximation is employed for the quantum propagator and the stationary phase approximation for certain integrals is utilized in this derivation. This analysis demonstrates that the divergence of the classical trajectories associated with the individual gaussian packets accounts for the changes in shape of the quantum wavefunction or density, as has been suggested on intuitive grounds by Heller. The method should be exact for quadratic potentials and this is verified by explicitly applying it for the harmonic oscillator example.     

Optimized Monte Carlo sampling in forward-backward semiclassical dynamics

Forward-backward semiclassical dynamics (FBSD) provides a rigorous and powerful methodology for calculating time correlation functions in condensed phase systems characterized by substantial quantum mechanical effects associated with zero-point motion, quantum dispersion, or identical particle exchange symmetries. The efficiency of these simulations arises from the use of classical trajectories to capture all dynamical information. However, full quantization of the density operator makes these calculations rather expensive compared to fully classical molecular dynamics simulations. This article discusses the convergence properties of various correlation functions and introduces an optimal Monte Carlo sampling scheme that leads to a significant reduction of statistical error. A simple and efficient procedure for normalizing the FBSD results is also discussed. Illustrative examples on model systems are presented.     

Interference and quantization in semiclassical response functions

Application of the Herman-Kluk semiclassical propagator to the calculation of spectroscopic response functions for anharmonic oscillators has demonstrated the quantitative accuracy of these approximate dynamics. In this approach, spectroscopic response functions are expressed as multiple phase-space integrals over pairs of classical trajectories and their associated stability matrices. Here we analyze the Herman-Kluk semiclassical approximation to a linear response function and determine the origin of the capacity of this method to reproduce quantum effects in a response function from classical dynamical information. Our analysis identifies those classical trajectories that contribute most significantly to the response function on different time scales. This finding motivates a procedure for computing the linear response function in which the interference between pairs of classical trajectories is treated approximately, resulting in an integral over a single average trajectory, as in a purely classical calculation.     

Finite temperature application of the corrected propagator method to reactive dynamics in a condensed-phase environment

The recently proposed mixed quantum-classical method is extended to applications at finite temperatures. The method is designed to treat complex systems consisting of a low-dimensional quantum part (the primary system) coupled to a dissipative bath described classically. The method is based on a formalism showing how to systematically correct the approximate zeroth-order evolution rule. The corrections are defined in terms of the total quantum Hamiltonian and are taken to the classical limit by introducing the frozen Gaussian approximation for the bath degrees of freedom. The evolution of the primary system is governed by the corrected propagator yielding the exact quantum dynamics. The method has been tested on a standard model system describing proton transfer in a condensed-phase environment: a symmetric double-well potential bilinearly coupled to a bath of harmonic oscillators. Flux correlation functions and thermal rate constants have been calculated at two different temperatures for a range of coupling strengths. The results have been compared to the fully quantum simulations of Topaler and Makri [J. Chem. Phys. 101, 7500 (1994)] with the real path integral method.     

Quantum-instanton evaluation of the kinetic isotope effects

A general quantum-mechanical method for computing kinetic isotope effects is presented. The method is based on the quantum-instanton approximation for the rate constant and on the path-integral Metropolis-Monte Carlo evaluation of the Boltzmann operator matrix elements. It computes the kinetic isotope effect directly, using a thermodynamic integration with respect to the mass of the isotope, thus avoiding the more computationally expensive process of computing the individual rate constants. The method should be more accurate than variational transition-state theories or the semiclassical instanton method since it does not assume a single tunneling path and does not use a semiclassical approximation of the Boltzmann operator. While the general Monte Carlo implementation makes the method accessible to systems with a large number of atoms, we present numerical results for the Eckart barrier and for the collinear and full three-dimensional isotope variants of the hydrogen exchange reaction H + H2 --> H2 + H. In all seven test cases, for temperatures between 250 and 600 K, the error of the quantum instanton approximation for the kinetic isotope effects is less than approximately 10%.     

Semiclassical calculation of the vibrational echo

The infrared echo measurement probes the time scales of the molecular motions that couple to a vibrational transition. Computation of the echo observable within rigorous quantum mechanics is problematic for systems with many degrees of freedom, motivating the development of semiclassical approximations to the nonlinear optical response. We present a semiclassical approximation to the echo observable, based on the Herman-Kluk propagator. This calculation requires averaging over a quantity generated by two pairs of classical trajectories and associated stability matrices, connected by a pair of phase-space jumps. Quantum, classical, and semiclassical echo calculations are compared for a thermal ensemble of noninteracting anharmonic oscillators. The semiclassical approach uses input from classical mechanics to reproduce the significant features of a complete, quantum mechanical calculation of the nonlinear response.     

Using the thermal Gaussian approximation for the Boltzmann operator in semiclassical initial value time correlation functions

The thermal Gaussian approximation (TGA) recently developed by Frantsuzov et al. [Chem. Phys. Lett. 381, 117 (2003)] has been demonstrated to be a practical way for approximating the Boltzmann operator exp(-betaH) for multidimensional systems. In this paper the TGA is combined with semiclassical (SC) initial value representations (IVRs) for thermal time correlation functions. Specifically, it is used with the linearized SC-IVR (LSC-IVR, equivalent to the classical Wigner model), and the "forward-backward semiclassical dynamics" approximation developed by Shao and Makri [J. Phys. Chem. A 103, 7753 (1999); 103, 9749 (1999)]. Use of the TGA with both of these approximate SC-IVRs allows the oscillatory part of the IVR to be integrated out explicitly, providing an extremely simple result that is readily applicable to large molecular systems. Calculation of the force-force autocorrelation for a strongly anharmonic oscillator demonstrates its accuracy, and calculation of the velocity autocorrelation function (and thus the diffusion coefficient) of liquid neon demonstrates its applicability.     

Semiclassical initial value series representation in the continuum limit: application to vibrational relaxation

A recently formulated continuum limit semiclassical initial value series representation (SCIVR) of the quantum dynamics of dissipative systems is applied to the study of vibrational relaxation of model harmonic and anharmonic oscillator systems. As is well known, the classical dynamics of dissipative systems may be described in terms of a generalized Langevin equation. The continuum limit SCIVR uses the Langevin trajectories as input, albeit with a quantum noise rather than a classical noise. Combining this development with the forward-backward form of the prefactor-free propagator leads to a tractable scheme for computing quantum thermal correlation functions. Here we present the first implementation of this continuum limit SCIVR series method to study two model problems of vibrational relaxation. Simulations of the dissipative harmonic oscillator system over a wide range of parameters demonstrate that at most only the first two terms in the SCIVR series are needed for convergence of the correlation function. The methodology is then applied to the vibrational relaxation of a dissipative Morse oscillator. Here, too, the SCIVR series converges rapidly as the first two terms are sufficient to provide the quantum mechanical relaxation with an estimated accuracy on the order of a few percent. The results in this case are compared with computations obtained using the classical Wigner approximation for the relaxation dynamics.     

Nonadiabatic electron wavepacket dynamics of molecules in an intense optical field: an ab initio electronic state study

A theory of quantum electron wavepacket dynamics that nonadiabatically couples with classical nuclear motions in intense optical fields is studied. The formalism is intended to track the laser-driven electron wavepackets in terms of the linear combination of configuration-state functions generated with ab initio molecular orbitals. Beginning with the total quantum Hamiltonian for electrons and nuclei in the vector potential of classical electromagnetic field, we reduce the Hamiltonian into a mixed quantum-classical representation by replacing the quantum nuclear momentum operators with the classical counterparts. This framework gives equations of motion for electron wavepackets in an intense laser field through the time dependent variational principle. On the other hand, a generalization of the Newtonian equations provides a matrix form of forces acting on the nuclei for nonadiabatic dynamics. A mean-field approximation to the force matrix reduces this higher order formalism to the semiclassical Ehrenfest theory in intense optical fields. To bring these theories into a practical quantum chemical package for general molecules, we have implemented the relevant ab initio algorithms in it. Some numerical results in the level of the semiclassical Ehrenfest-type theory with explicit use of the nuclear kinematic (derivative) coupling and the velocity form for the optical interaction are presented.     

Generalization of classical mechanics for nuclear motions on nonadiabatically coupled potential energy surfaces in chemical reactions

Classical trajectory study of nuclear motion on the Born-Oppenheimer potential energy surfaces is now one of the standard methods of chemical dynamics. In particular, this approach is inevitable in the studies of large molecular systems. However, as soon as more than a single potential energy surface is involved due to nonadiabatic coupling, such a naive application of classical mechanics loses its theoretical foundation. This is a classic and fundamental issue in the foundation of chemistry. To cope with this problem, we propose a generalization of classical mechanics that provides a path even in cases where multiple potential energy surfaces are involved in a single event and the Born-Oppenheimer approximation breaks down. This generalization is made by diagonalization of the matrix representation of nuclear forces in nonadiabatic dynamics, which is derived from a mixed quantum-classical representation of the electron-nucleus entangled Hamiltonian [Takatsuka, K. J. Chem. Phys. 2006, 124, 064111]. A manifestation of quantum fluctuation on a classical subsystem that directly contacts with a quantum subsystem is discussed. We also show that the Hamiltonian thus represented gives a theoretical foundation to examine the validity of the so-called semiclassical Ehrenfest theory (or mean-field theory) for electron quantum wavepacket dynamics, and indeed, it is pointed out that the electronic Hamiltonian to be used in this theory should be slightly modified.     

Renormalization of the frozen Gaussian approximation to the quantum propagator

The frozen Gaussian approximation to the quantum propagator may be a viable method for obtaining "on the fly" quantum dynamical information on systems with many degrees of freedom. However, it has two severe limitations, it rapidly loses normalization and one needs to know the Gaussian averaged potential, hence it is not a purely local theory in the force field. These limitations are in principle remedied by using the Herman-Kluk (HK) form for the semiclassical propagator. The HK propagator approximately conserves unitarity for relatively long times and depends only locally on the bare potential and its second derivatives. However, the HK propagator involves a much more expensive computation due to the need for evaluating the monodromy matrix elements. In this paper, we (a) derive a new formula for the normalization integral based on a prefactor free HK propagator which is amenable to "on the fly" computations; (b) show that a frozen Gaussian version of the normalization integral is not readily computable "on the fly"; (c) provide a new insight into how the HK prefactor leads to approximate unitarity; and (d) how one may construct a prefactor free approximation which combines the advantages of the frozen Gaussian and the HK propagators. The theoretical developments are backed by numerical examples on a Morse oscillator and a quartic double well potential.     

Matching-pursuit split-operator Fourier-transform simulations of excited-state intramolecular proton transfer in 2-(2'-hydroxyphenyl)-oxazole

The excited-state intramolecular proton-transfer dynamics associated with the keto-enolic tautomerization reaction in 2-(2(')-hydroxyphenyl)-oxazole is simulated according to a numerically exact quantum-dynamics propagation method and a full-dimensional excited-state potential energy surface, based on an ab initio reaction surface Hamiltonian. The reported simulations involve the propagation of 35-dimensional wave packets according to the recently developed matching-pursuit/split-operator-Fourier-transform (MP/SOFT) method by Wu and Batista. The underlying propagation scheme recursively applies the time-evolution operator as defined by the Trotter expansion to second order accuracy in dynamically adaptive coherent-state expansions. Computations of time-dependent survival amplitudes, photoabsorption cross sections, and time-dependent reactant(product) populations are compared to the corresponding calculations based on semiclassical approaches, including the Herman-Kluk semiclassical initial value representation method. The reported results demonstrate the capabilities of the MP/SOFT method as a valuble computational tool to study ultrafast reaction dynamics in polyatomic systems as well as to validate semiclassical simulations of complex (nonintegrable) quantum dynamics in multidimensional model systems.     

Long-time and unitary properties of semiclassical initial value representations

We numerically compare the semiclassical "frozen Gaussian" Herman-Kluk propagator [Chem. Phys. 91, 27 (1984)] and the "thawed Gaussian" propagator put forward recently by Baranger et al. [J. Phys. A 34, 7227 (2001)] by studying the quantum dynamics in some nonlinear one-dimensional potentials. The reasons for the lack of long-time accuracy and norm conservation in the latter method are uncovered. We amend the thawed Gaussian propagator with a global harmonic approximation for the stability of the trajectories and demonstrate that this revised propagator is a true alternative to the Herman-Kluk propagator with similar accuracy.     

Independent trajectory implementation of the semiclassical Liouville method: application to multidimensional reaction dynamics

We describe an independent trajectory implementation of semiclassical Liouville method for simulating quantum processes using classical trajectories. In this approach, a single ensemble of trajectories describes all semiclassical density matrix elements of a coupled electronic state problem, with the ensemble evolving classically under a single reference Hamiltonian chosen on the basis of physical grounds. In this paper, we introduce an additional uncoupled trajectory approximation, allowing the members of the ensemble to evolve independently of one another and eliminating the major computational costs of our previous coupled trajectory implementation. The accuracy of the method is demonstrated for model one-dimensional problems. In addition, the approach is applied to the chemical reaction dynamics of a collinear triatomic system, yielding excellent agreement with exact calculations. This method allows molecular dynamics involving coupled electronic surfaces to be modeled with essentially the same effort as classical molecular dynamics and ensemble averaging.     

Complete dynamics of a two-level quantum system in interaction with a radiation field

Using the general theory of systems of linear differential equations with periodic coefficients, we derive a complete set of solutions of the equations of motion of a two-level quantum system in interaction with a classical radiation field. The structure of the solutions is investigated and a method of obtaining approximate solutions is discussed. These solutions are compared quantitatively with those in the rotating-wave approximation in dependence of field amplitude, detuning from resonance, time of interaction between two-level system and field, and initial state of the two-level system. In an appendix it is shown that the semiclassical dynamics of a harmonically driven two-level system may be derived from an associated fully quantum-mechanical motion by an asymptotic limit.     

A coherent state approach to semiclassical nonadiabatic dynamics

A semiclassical (SC) approximation to the quantum mechanical propagator for nonadiabatic systems is derived. Our derivation starts with an exact path integral expression that uses canonical coherent states for the nuclear degrees of freedom and spin coherent states for the electronic degrees of freedom. A stationary path approximation (SPA) is then applied to the path integral to obtain the SC approximation. The SPA results in complex classical trajectories of both nuclear and electronic degrees of freedom and a double ended boundary condition. The root search problem is solved using the previously proposed "real trajectory local search" algorithm. The SC approximation is tested on three simple one dimensional two-state systems proposed by Tully [J. Chem. Phys. 93, 1061 (1990)], and the SC results are compared to Ehrenfest and surface hopping predictions. Excellent agreement with quantum results is reached when the SC trajectory is far away from caustics. We discuss the origin of caustics in this SC formalism and the strengths and weaknesses of this approach.     

Efficient evaluation of the accuracy of molecular quantum dynamics on an approximate analytical or interpolated ab initio potential energy surface

Ab initio electronic structure methods have reached a satisfactory accuracy for the calculation of static properties, but remain too expensive for quantum dynamical calculations. Recently, an efficient semiclassical method was proposed to evaluate the accuracy of quantum dynamics on an approximate potential without having to perform the expensive quantum dynamics on the accurate potential. Here, this method is applied for the first time to evaluate the accuracy of quantum dynamics on an approximate analytical or interpolated potential in comparison to the quantum dynamics on an accurate potential obtained by an ab initio electronic structure method. Specifically, the vibrational dynamics of H₂ on a Morse potential is compared with that on the full CI potential, and the photodissociation dynamics of CO₂ on a LEPS potential with that on the excited ¹Π surface computed at the EOM‐CCSD/aug‐cc‐pVDZ level of theory. Finally, the effect of discretization of a potential energy surface on the quantum dynamics is evaluated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2426–2435, 2010