The Wigner phase space formulation of quantum mechanics is a complete framework for quantum dynamic calculations that elegantly highlights connections with classical dynamics. In this series of two articles, building upon previous efforts, we derive the full hierarchy of approximate semiclassical (SC) dynamic methods for adiabatic and non-adiabatic problems in Wigner phase space. In Paper I, focusing on adiabatic single surface processes, we derive the well-known double Herman-Kluk (DHK) approximation for real-time correlation functions in Wigner phase space and connect it to the linearized SC (LSC) approximation through a stationary phase approximation. We exploit this relationship to introduce a new hybrid SC method, termed Adiabatic Hybrid Wigner Dynamics (AHWD) that allows for a few important "system" degrees of freedom (dofs) to be treated at the DHK level, while treating the rest of the dofs (the "bath") at the LSC level. AHWD is shown to accurately capture quantum interference effects in models of coupled oscillators and the decoherence of vibrational probability density of a model I2 Morse oscillator coupled to an Ohmic thermal bath. We show that AHWD significantly mitigates the sign problem and employs reduced dimensional prefactors bringing calculations of complex system-bath problems within the reach of SC methods. Paper II focuses on extending this hybrid SC dynamics to nonadiabatic processes.
The Wigner phase space formulation of quantum mechanics is a complete framework for quantum dynamic calculations that elegantly highlights connections with classical dynamics. In this series of two articles, building upon previous efforts, we derive the full hierarchy of approximate semiclassical (SC) dynamic methods for adiabatic and non-adiabatic problems in Wigner phase space. In Paper I, focusing on adiabatic single surface processes, we derive the well-known double Herman-Kluk (DHK) approximation for real-time correlation functions in Wigner phase space and connect it to the linearized SC (LSC) approximation through a stationary phase approximation. We exploit this relationship to introduce a new hybrid SC method, termed Adiabatic Hybrid Wigner Dynamics (AHWD) that allows for a few important "system" degrees of freedom (dofs) to be treated at the DHK level, while treating the rest of the dofs (the "bath") at the LSC level. AHWD is shown to accurately capture quantum interference effects in models of coupled oscillators and the decoherence of vibrational probability density of a model I2 Morse oscillator coupled to an Ohmic thermal bath. We show that AHWD significantly mitigates the sign problem and employs reduced dimensional prefactors bringing calculations of complex system-bath problems within the reach of SC methods. Paper II focuses on extending this hybrid SC dynamics to nonadiabatic processes.
ABSTRACT: The concept of open quantum system is very broad and it is related to the ability of measuring only certain degree of freedoms of a particular system. Although this idea is relatively simple, the separation between the system of interest, the degree of freedoms that are accesible experimentally, and the reservoir, the degree of freedoms that are not accesible experimentally, is not always clear. For instance, in molecular systems, electronic spectroscopy has access only to electronic degree of freedoms so that the nuclear and vibrational degree of freedom become the reservoir, this makes its description not trivial. This scenario leads to have non-trivial and structured reservoirs and to develop powerful tools to analyze them. The hallmark of non-trivial and structured reservoirs is the non-Markovian dynamics. By translating the Feynman and Vernon influence functional approach into phase-space representation, we develop two theories of semiclassical evolution of the Wigner function of the system of interest that incorporate non-Markovian dynamics and highly non-trivial quantum effects such as non-locality of quantum mechanics: (i) We translate the Caldeira-Leggett model into phase-space representation of quantum mechanics and (ii) we consider the possibility of having non-linear baths and therefore, truly quantum reservoirs. During the last forty years the study of energy loss and coherence in quantum systems has been based on the Ullersma-Caldeira-Leggett model, a model that describes the environment of quantum systems of interest as a collection of harmonic oscillators with classical evolution. We constructed this model in the Wigner-Weyl representation of quantum mechanics and discuss the classical nature of the evolution of the bath modes and the semiclassical evolution of the central system. As an application of the semiclassical Wigner propagator, the non-Markovian time evolution under the Morse potencial is analyzed. There, it is clear how decohering processes shrink the propagator to smaller regions of phase space implying that the dynamics become more local, i.e., more classical. The current level of experimentation and control of physical systems have called into question the validity of the model in, one hand, molecular systems (e.g. photosynthetic complexes immersed in solvents, chemical systems in liquid phase or gas and manipulated with intense laser pulses) and, in the other hand, solid state systems (e.g. Josephson junctions, spins in quantum dots or "spinning ice"). Given the importance of these systems in the development of new quantum technologies and in the understanding of quantum phenomena in mesoscopic systems, it has become necessary to develop new models of the environment and efficient methodologies with quantitative prediction power. However, some of them are artificial modifications of the Ullersma-Caldeira-Leggett model without solid and clean physical support. Here we formulated a general framework that allows the study of quantum correlations in quantum systems in the presence of non-harmonic thermal baths (e.g. baths formed by strongly coupled diatomic molecules). This formulation will allow a more precise and quantitative description of processes such as the transport of excitons in photosynthetic complexes, the transfer of heat in solid state devices, among others. Results clearly show the non-classical time dynamics of the bath modes. The implementation of this particular theory remains, however, as a challenge.
We present a comprehensive study of semiclassical phase-space propagation in the Wigner representation, emphasizing numerical applications, in particular as an initial-value representation. Two semiclassical approximation schemes are discussed. The propagator of the Wigner function based on van Vleck's approximation replaces the Liouville propagator by a quantum spot with an oscillatory pattern reflecting the interference between pairs of classical trajectories. Employing phase-space path integration instead, caustics in the quantum spot are resolved in terms of Airy functions. We apply both to two benchmark models of nonlinear molecular potentials, the Morse oscillator and the quartic double well, to test them in standard tasks such as computing autocorrelation functions and propagating coherent states. The performance of semiclassical Wigner propagation is very good even in the presence of marked quantum effects, e.g., in coherent tunneling and in propagating Schrodinger cat states, and of classical chaos in four-dimensional phase space. We suggest options for an effective numerical implementation of our method and for integrating it in Monte-Carlo-Metropolis algorithms suitable for high-dimensional systems.
We formulate various semiclassical propagators for the Wigner phase space representation from a unified point of view. As is shown in several studies, the Moyal equation, which is an equation of motion for the Wigner distribution function, can be regarded as the Schrödinger equation of an extended Hamiltonian system where its "position" and "momentum" correspond to the middle point of two points of the original phase space and the difference between them, respectively. Then we show that various phase-space semiclassical propagators can be formulated just by applying existing semiclassical propagators to the extended system. As a result, a phase space version of the Van Vleck propagator, the initial-value Van Vleck propagator, the Herman-Kluk propagator, and the thawed Gaussian approximation are obtained. In addition, we numerically compare the initial-value phase-space Van Vleck propagator, the phase-space Herman-Kluk propagator, and the classical mechanical propagation as approximation methods for the time propagation of the Wigner distribution function in terms of both accuracy and convergence speed. As a result, we find that the convergence speed of the Van Vleck propagator is far slower than others as is the case of the Hilbert space, and the Herman-Kluk propagator keeps its accuracy for a long period compared with the classical mechanical propagation while the convergence speed of the latter is faster than the former.
The semiclassical Herman–Kluk initial value approximation of the quantum-mechanical propagator, which is heavily used in dynamical studies of atomic and molecular systems, takes into account nonlinearity of the underlying classical dynamics by using multiple classical initial value solutions. In applications to the propagation of Gaussian wavepackets non-Gaussian distortions in the course of time and therefore a realistic description of the quantum dynamics are possible. This is in contrast to the single-trajectory approximations investigated in a recent paper by Baranger et al , one of which was erroneously termed Herman–Kluk approximation.
We present a new approach to calculate real-time quantum dynamics in complex systems. The formalism is based on the partitioning of a system's environment into "core" and "reservoir" modes with the former to be treated quantum mechanically and the latter classically. The presented method only requires the calculation of the system's reduced density matrix averaged over the quantum core degrees of freedom which is then coupled to a classically evolved reservoir to treat the remaining modes. We demonstrate our approach by applying it to the spin-boson problem using the noninteracting blip approximation to treat the system and core, and Ehrenfest dynamics to treat the reservoir. The resulting hybrid methodology is accurate for both fast and slow baths, since it naturally reduces to its composite methods in their respective regimes of validity. In addition, our combined method is shown to yield good results in intermediate regimes where neither approximation alone is accurate and to perform equally well for both strong and weak system-bath coupling. Our approach therefore provides an accurate and efficient methodology for calculating quantum dynamics in complex systems.
In this work we provide a complete model of semiclassical theories by including back-reaction and correlation into the picture. We specially aim at the interaction between light and a two-level atom, and we also illustrate it via the coupling of two harmonic oscillators. Quantum and classical systems are treated on the same grounds via the Wigner-Weyl phase-space correspondence of the quantum theory. We show that this model provides a suitable mixture of the quantum and classical degrees of freedom, including the fact that the evolution transfers nonclassical features to the classical subsystem and nonquantum behavior to the quantum subsystem. In that sense, we can no longer distinguish between classical and quantum variables and we need to talk about a hybrid model.