Semiclassical calculation of chemical reaction dynamics via wavepacket correlation functions

Computing Quantum Mean Values in the Deep Chaotic Regime

We study the time evolution of mean values of quantum operators in a regime plagued by two difficulties: the smallness of ℏ and the presence of strong and ubiquitous classical chaos. While numerics become too computationally expensive for purely quantum calculations as ℏ→0, methods that take advantage of the smallness of ℏ-that is, semiclassical methods-suffer from both conceptual and practical difficulties in the deep chaotic regime. We implement an approach which addresses these conceptual problems, leading to a deeper understanding of the origin of the interference contributions to the operator's mean value. We show that in the deep chaotic regime our approach is capable of unprecedented accuracy, while a standard semiclassical method (the Herman-Kluk propagator) produces only numerical noise. Our work paves the way to the development and employment of more efficient and accurate methods for quantum simulations of systems with strongly chaotic classical limits.

Semiclassical initial value representation: From Møller to Miller. II

As shown by W. H. Miller in a seminal article [J. Chem. Phys. 53, 3578 (1970)], the most convenient and accurate semiclassical (SC) theory of molecular scattering in action-angle coordinates is based on the initial value representation (IVR) and the use of shifted angles, which are different from the natural angles usually used in the quantum and classical treatments. Here, we show for an inelastic molecular collision that the initial and final shifted angles define three-segment classical paths that are exactly those involved in the classical-limit of Tannor-Weeks quantum scattering theory [J. Chem. Phys. 98, 3884 (1993)], provided that the translational wave packets |g₊⟩ and |g_-⟩ into play in this theory are both taken at |0⟩. Assuming this to be the case, using van Vleck propagators, and applying the stationary phase approximation, Miller's SCIVR expression of S-matrix elements is found, with an additional cut-off factor canceling the energetically forbidden transition probabilities. This factor, however, is close to unity in most practical cases. Furthermore, these developments show that the Møller operators underlie Miller's formulation, thus confirming, for molecular collisions, the results recently established in the simpler case of light-induced rotational transitions [L. Bonnet, J. Chem. Phys. 153, 174102 (2020)]. Last but not least, we show, based on the previous results, that for processes involving long-range anisotropic forces, implementing the Skinner-Miller method [Chem. Phys. Lett. 300, 20 (1999)] in shifted coordinates makes its predictions both easier and more accurate than in natural coordinates.

A Semiclassical Framework for Mixed Quantum Classical Dynamics

Semiclassical (SC) approximations for quantum dynamic simulations in complex chemical systems range from rigorously accurate methods that are computationally expensive to methods that exhibit near-classical scaling with system size but are limited in their ability to describe quantum effects. In practical studies of high-dimensional reactions, neither extreme is the best choice: frequently a high-level quantum mechanical description is only required for a handful of modes, while the majority of environment modes that do not play a key role in the reactive event of interest are well served with a lower level of theory. In this feature, we introduce modified Filinov filtration as a powerful tool to construct mixed quantum-classical SC theories where different subsystems can be quantized to different extents without introducing ad hoc intersubsystem interaction terms. We demonstrate that these Filinov-based SC methods can systematically tune between quantum and classical limit SC behavior, offering a practical way forward to accurate and computationally efficient simulations of high-dimensional quantum processes.

Semiclassical descriptions of rotational transitions in natural and shifted angles: Analysis of unexpected results

In the semiclassical theory of rotational transitions, S-matrix elements are expressed as integrals over initial and final angles of probability amplitudes calculated along the classical paths joining these angles, before final passage to an initial value representation [W. H. Miller, J. Phys. Chem. A 105, 2942 (2001)]. These angles can be either natural angles fixing the orientation of the rotor or angles shifted with respect to the previous ones so as to vary only within the interaction region causing the transitions. The two approaches, however, were recently shown to lead to different predictions. While the theory in natural angles lacks precision and exhibits unphysical behavior, the theory in shifted angles is much more accurate and physically well behaved [L. Bonnet, J. Chem. Phys. 153, 174102 (2020)]. The present work is devoted to the analysis of this unexpected finding.

When classical trajectories get to quantum accuracy: II. The scattering of rotationally excited H2 on Pd(111)

The classical trajectory method in a quantum spirit assigns statistical weights to classical paths on the basis of two semiclassical corrections: Gaussian binning and the adiabaticity correction.

Theoretical Study of Barrierless Chemical Reactions Involving Nearly Elastic Rebound: The Case of S(1D) + X2, X = H, D

For some values of the total angular momentum consistent with reaction, the title processes involve nonreactive trajectories proceeding through a single rebound mechanism during which the internal motion of the reagent diatom is nearly unperturbed. When such paths are in a significant amount, the classical reaction probability is found to be markedly lower than the quantum mechanical one. This finding was recently attributed to an unusual quantum effect called diffraction-mediated trapping, and a semiclassical correction was proposed in order to take into account this effect in the classical trajectory method. In the present work, we apply the resulting approach to the calculation of opacity functions as well as total and state-resolved integral cross sections (ICSs) and compare the values obtained with exact quantum ones, most of which are new. As the title reactions proceed through a deep insertion well, mean potential statistical calculations are also presented. Seven values of the collision energy, ranging from 30 to 1127 K, are considered. Two remarkable facts stand out: (i) The corrected classical treatment strongly improves the accuracy of the opacity function as compared to the usual classical treatment. When the entrance transition state is tight, however, those trajectories crossing it with a bending vibrational energy below the zero point energy must be discarded. (ii) The quantum opacity function, particularly its cutoff, is finely reproduced by the statistical approach. Consequently, the total ICS is also very well described by the two previous approximate methods. These, however, do not predict state-resolved ICSs with the same accuracy, proving thereby that (i) one or several genuine quantum effects involved in the dynamics are missed by the corrected classical treatment and (ii) the dynamics are not fully statistical.

Time-dependent wave packet simulations of transport through Aharanov-Bohm rings with an embedded quantum dot

We have performed time-dependent wave packet simulations of realistic Aharonov-Bohm (AB) devices with a quantum dot embedded in one of the arms of the interferometer. The AB ring can function as a measurement device for the intrinsic transmission phase through the quantum dot, however, care has to be taken in analyzing the influence of scattering processes in the junctions of the interferometer arms. We consider a harmonic quantum dot and show how the Darwin-Fock spectrum emerges as a unique pattern in the interference fringes of the AB oscillations.

Communication: importance sampling including path correlation in semiclassical initial value representation calculations for time correlation functions

Full semiclassical (SC) initial value representation (IVR) for time correlation functions involves a double phase space average over a set of two phase points, each of which evolves along a classical path. Conventionally, the two initial phase points are sampled independently for all degrees of freedom (DOF) in the Monte Carlo procedure. Here, we present an efficient importance sampling scheme by including the path correlation between the two initial phase points for the bath DOF, which greatly improves the performance of the SC-IVR calculations for large molecular systems. Satisfactory convergence in the study of quantum coherence in vibrational relaxation has been achieved for a benchmark system-bath model with up to 21 DOF.

SEMICLASSICAL FROZEN GAUSSIAN PROPAGATION METHOD FOR ELECTRONICALLY NONADIABATIC CHEMICAL DYNAMICS: MØLLER OPERATOR FORMULATION AND INCORPORATION OF THE ZHU-NAKAMURA THEORY

The title theory is developed by combining the Herman–Kluk semiclassical theory for adiabatic propagation on single potential energy surface and the Zhu–Nakamura theory for nonadiabatic transition. A fairly simple expression for the propagator based on classical trajectories is derived using the Møller operator formulation of scattering theory. The theory takes into account almost all quantum effects that occur during nonadiabatic transition, especially at low energies and is expected to be applicable to general chemical dynamics of high dimensions. Application to a two-dimensional model system shows that the theory works well for the propagation duration of several molecular vibrational periods and wide wave packet energy range.

A semiclassical correction for quantum mechanical energy levels

We propose a semiclassical method for correcting molecular energy levels obtained from a quantum mechanical variational calculation. A variational calculation gives the energy level (i.e., eigenvalue) as the expectation value of the molecular Hamiltonian <phi/H/phi>, where /phi> is the trial wave function. The true (i.e., exact) eigenvalue E can thus be expressed as this variational result plus a correction, i.e., E=<phi/H/phi>+DeltaE, the correction being due to the lack of exactness of the trial wave function. A formally exact expression for DeltaE is usually given (via Löwdin partitioning methodology) in terms of the Greens function of the Hamiltonian projected onto the orthogonal complement of /phi>. Formal treatment of this expression (using Brillouin-Wigner perturbation theory to infinite order) leads to an expression for DeltaE that involves matrix elements of the Greens function for the unprojected, i.e., full molecular Hamiltonian, which can then be approximated semiclassically. (Specifically, the Greens function is expressed as the Fourier transform of the quantum mechanical time evolution operator, e(-iHt/variant Planck's over 2pi), which in turn is approximated by using an initial value representation of semiclassical theory.) Calculations for several test problems (a one dimensional quartic potential, and vibrational energy levels of H(2)O and H(2)CO) clearly support our proposition that the error in the total eigenvalue E arises solely due to the semiclassical error in approximating DeltaE, which is usually a small fraction of the total energy E itself.

Gaussian approximation for the structure function in semiclassical forward-backward initial value representations of time correlation functions

Initial value representations (IVRs) of semiclassical (SC) theory provide a general approach for adding quantum mechanical effects to classical molecular dynamics simulations of large molecular systems. Of the various versions of SC-IVR methodology for evaluating time correlation functions, the Fourier transform forward-backward (FB) approach is the simplest one that is able to describe true quantum coherence effects, so it is of considerable importance to find efficient and systematic ways for implementing it. It is shown in this paper that a Gaussian approximation for the "structure function"-the dependence of the correlation function on the (typically) momentum jump parameter-provides an efficient and accurate way for doing so. The approach is illustrated by an application to the time-dependent radial distribution function of I(2) (after photoexcitation) in a cluster of (up to 16) argon atoms.

Continuum limit semiclassical initial value representation for dissipative systems

In this paper, we consider a dissipative system in which the system is coupled linearly to a harmonic bath. In the continuum limit, the bath is defined via a spectral density and the classical system dynamics is given in terms of a generalized Langevin equation. Using the path integral formulation and factorized initial conditions, it is well known that one can integrate out the harmonic bath, leaving only a path integral over the system degrees of freedom. However, the semiclassical initial value representation treatment of dissipative systems has usually been limited to a discretized treatment of the bath in terms of a finite number of bath oscillators. In this paper, the continuum limit of the semiclassical initial value representation is derived for dissipative systems. As in the path integral, the action is modified with an added nonlocal term, which expresses the influence of the bath on the dynamics. The first order correction term to the semiclassical initial value approximation is also derived in the continuum limit.

A Hybrid Hydrodynamic−Liouvillian Approach to Mixed Quantum−Classical Dynamics: Application to Tunneling in a Double Well

The hybrid quantum-classical approach of Burghardt and Parlant [Burghardt, I.; Parlant, G. J. Chem. Phys. 2004, 120, 3055], referred to here as the quantum-classical moment (QCM) approach, is demonstrated for the dynamics of a quantum double well coupled to a classical harmonic coordinate. The approach combines the quantum hydrodynamic and classical Liouvillian representations by the construction of a particular type of moments (that is, partial hydrodynamic moments) whose evolution is determined by a hierarchy of coupled equations. For pure states, which are at the center of the present study, this hierarchy terminates at the first order. In the Lagrangian picture, the deterministic trajectories result in dynamics which is Hamiltonian in the classical subspace, while the projection onto the quantum subspace evolves under a generalized hydrodynamic force. Importantly, this force also depends upon the classical (Q, P) variables. The present application demonstrates the tunneling dynamics in both the Eulerian and Lagrangian representations. The method is exact if the classical subspace is harmonic, as is the case for the systems studied here.

New coherent state representation for the imaginary time propagator with applications to forward-backward semiclassical initial value representations of correlation functions

There have been quite a few attempts in recent years to provide an initial value coherent state representation for the imaginary time propagator exp(-betaH). The most notable is the recent time evolving Gaussian approximation of Frantsuzov and Mandelshtam [J. Chem. Phys. 121, 9247 (2004)] which may be considered as an expansion of the imaginary time propagator in terms of coherent states whose momentum is zero. In this paper, a similar but different expression is developed in which exp(-betaH) is represented in a series whose terms are weighted phase space averages of coherent states. Such a representation allows for the formulation of a new and simplified forward-backward semiclassical initial value representation expression for thermal correlation functions.

Semiclassical initial value series solution of the spin boson problem

A numerical solution for the quantum dynamics of the spin boson problem is obtained using the semiclassical initial value series representation approach to the quantum dynamics. The zeroth order term of the series is computed using the new forward-backward representation for correlation functions presented in the preceding adjacent paper. This leads to a rapid convergence of the Monte Carlo sampling, as compared to previous attempts. The zeroth order results are already quite accurate. The first order term of the series is small, demonstrating the rapid convergence of the semiclassical initial value representation series. This is the first time that the first order term in the semiclassical initial value representation series has been converged for systems with the order of 50 degrees of freedom.

Frozen Gaussian Wavepacket Study of the Ground State of the He Atom

The Rayleigh-Ritz functional is used in conjunction with an approximate time evolution to improve ab initio estimates of ground-state energies. The improvement is due in part to the introduction of a novel variational "normalization function" for the approximate propagator. An additional variational parameter was introduced in the form of a constant shift energy of the Hamiltonian. The approximate propagator used was the frozen Gaussian propagator; however, the trajectories evolved on the coherent-state averaged Hamiltonian (Q representation). For Coulombic forces, this removes the singularity, easing the computation. An additional variational parameter was the width parameter used for the coherent states appearing in the frozen Gaussian propagator. Using an initial combination of nine Gaussian functions for He, with an initial energy of -2.5115 au, the variational method, with a very short time interval of integration, led to an improved energy of -2.81 ± 0.04 au.

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.

A prefactor free semiclassical initial value series representation of the propagator

A new class of prefactor free semiclassical initial value representations (SCIVR) of the quantum propagator is presented. The derivation is based on the physically motivated demand, that on the average in phase space and in time, the propagator obey the exact quantum equation of motion. The resulting SCIVR series representation of the exact quantum propagator is also free of prefactors. When using a constant width parameter, the prefactor free SCIVR propagator is identical to the frozen Gaussian propagator of Heller [J. Chem. Phys. 75, 2923 (1981)]. A numerical study of the prefactor free SCIVR series is presented for scattering through a double slit potential, a system studied extensively previously by Gelabert et al. [J. Chem. Phys. 114, 2572 (2001)]. As a basis for comparison, the SCIVR series is also computed using the optimized Herman-Kluk SCIVR. We find that the sum of the zeroth order and the first order terms in the series suffice for an accurate determination of the diffraction pattern. The same exercise, but using the prefactor free propagator series needs also the second order term in the series, however the numerical effort is not greater than that needed for the Herman-Kluk propagator, since one does not need to compute the monodromy matrix elements at each point in time. The numerical advantage of the prefactor free propagator grows with increasing dimensionality of the problem.

Forward-backward semiclassical initial value series representation of quantum correlation functions

The forward-backward (FB) approximation as applied to semiclassical initial value representations (SCIVR's) has enabled the practical application of the SCIVR methodology to systems with many degrees of freedom. However, to date a systematic representation of the exact quantum dynamics in terms of the FB-SCIVR has proven elusive. In this paper, we provide a new derivation of a forward-backward phase space SCIVR expression (FBPS-SCIVR) derived previously by Thompson and Makri [Phys. Rev. E 59, R4729 (1999)]. This enables us to represent quantum correlation functions exactly in terms of a series whose leading order term is the FBPS-SCIVR expression. Numerical examples for systems with over 50 degrees of freedom are presented for the spin boson problem. Comparison of the FBPS-SCIVR with the numerically exact results of Wang [J. Chem. Phys. 113, 9948 (2000)] obtained using a multiconfigurational time dependent method shows that the leading order FBPS-SCIVR term already provides an excellent approximation.

A new time evolving Gaussian series representation of the imaginary time propagator

Frantsuzov and Mandelshtam [J. Chem. Phys. 121, 9247 (2004)] have recently demonstrated that a time evolving Gaussian approximation (TEGA) to the imaginary time propagator exp(-betaH) is useful for numerical computations of anharmonically coupled systems with many degrees of freedom. In this paper we derive a new exact series representation for the imaginary time propagator whose leading order term is the TEGA. One can thus use the TEGA not only as an approximation but also to obtain the exact imaginary time propagator. We also show how the TEGA may be generalized to provide a family of TEGA's. Finally, we find that the equations of motion governing the evolution of the center and width of the Gaussian may be thought of as introducing a quantum friction term to the classical evolution equations.