A new semiclassical theory is formulated by introducing a new distribution function to be propagated in phase space. Primarily this theory is intended to describe molecular collisions such as reactive scattering, but is also applicable to intramolecular processes including bound state problems. This new distribution function is a kind of characteristic function in the terminology of the phase-space-quantum theory developed by Weyl, Wigner, Moyal and others, and is called a dynamical characteristic function (DCF). The equation of motion and some properties of the DCF are presented. The essential feature of our theory is that the semiclassical DCF can be propagated using classical trajectories. This distribution function carried information not only on the density, but also on the phase (action integral). In this sense this theory can be regarded as a kind of phase-space-path-integral formalism. The semiclassical approximation is analyzed from the viewpoint of path integrals. The Wigner distribution function is also discussed in the context of our analysis. As a theoretical application, the semiclassical DCF is shown to reproduce the expression of the periodic orbit theory for bound states. Also, as a simple numerical application, a propagation of the DCF on a nonquadratic (Morse) potential is presented.
A new semiclassical theory is formulated by introducing a new distribution function to be propagated in phase space. Primarily this theory is intended to describe molecular collisions such as reactive scattering, but is also applicable to intramolecular processes including bound state problems. This new distribution function is a kind of characteristic function in the terminology of the phase-space-quantum theory developed by Weyl, Wigner, Moyal and others, and is called a dynamical characteristic function (DCF). The equation of motion and some properties of the DCF are presented. The essential feature of our theory is that the semiclassical DCF can be propagated using classical trajectories. This distribution function carried information not only on the density, but also on the phase (action integral). In this sense this theory can be regarded as a kind of phase-space-path-integral formalism. The semiclassical approximation is analyzed from the viewpoint of path integrals. The Wigner distribution function is also discussed in the context of our analysis. As a theoretical application, the semiclassical DCF is shown to reproduce the expression of the periodic orbit theory for bound states. Also, as a simple numerical application, a propagation of the DCF on a nonquadratic (Morse) potential is presented.
A detailed framework for a new formalism of phase-space path integral is presented, the outline of the theory of which has been reported elsewhere [K. Takatsuka, Phys. Rev. Lett. 61, 503 (1988)]. This path integral is described in terms of a phase-space distribution function, which was proposed earlier by the present author under the name of the dynamical characteristic function (DCF). The DCF is characterized by two independent wave functions and with two phase spaces. In the present paper, we extend the DCF by utilizing the Feynman kernel in place of the two independent wave functions. This is called the identity DCF. The dynamics of the DCF for general wave packets is reduced to that of the identity DCF. Some characteristics of the identity DCF are presented. For example, the quantum-mechanical time-evolution operator is represented in terms of the identity DCF and of a complete set of quantum q (coordinate) and p (momentum) operators. A semiclassical theory for the identity DCF is also developed, its primary aim being the study of heavy-particle dynamics such as molecular collisions.It is shown that the propagation of the identity DCF is represented by classical trajectories, each of which is associated with an amplitude factor and with the action integral as the quantum phase factor. This amplitude factor is free from singularity, unlike the primitive semiclassical expression of the Feynman kernel, which is well known to diverge at caustics. The physical significance of the identity DCF is analyzed in the context of chaotic dynamics. In particular, it is shown that its amplitude factor is an explicit function of the eigenvalues of a Jacobian matrix [\ensuremath{\partial}(qf,pf)/\ensuremath{\partial}(qi,pi)], where (qi,p) and (qf,pf) are the initial and final points, respectively, of a classical trajectory. These eigenvalues are well known in the study of the classical chaos of Hamilton systems and have been used to construct a measure of the linear (and local) stability of an orbit, the Liapounov characteristic exponents, and so on. Further, if the identity DCF is applied to a periodic orbit, the amplitude factor becomes equivalent to the square root of the determinant of a matrix, which is the so-called monodromy matrix minus the unit matrix. This matrix is one of the standard quantities used in bifurcation theory to investigate the stability of fixed points. In fact, the amplitude factor is closely related to the so-called residue, which was proposed by Greene [J. Math. Phys. 20, 1183 (1979)] in order to predict the occurrence of global chaos in classical systems. It is emphasized that all these classical quantities arise naturally in our representation of quantum mechanics.
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.
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.
Abstract. We establish the difference between the propagation of semiclassical Wigner functions and classical Liouville propagation. First we re-discuss the semiclassical limit for the propagator of Wigner functions, which on its own leads to their classical propagation. Then, via stationary phase evaluation of the full integral evolution equation, using the semiclassical expressions of Wigner functions, we provide the correct geometrical prescription for their semiclassical propagation. This is determined by the classical trajectories of the tips of the chords defined by the initial semiclassical Wigner function and centered on their arguments, in contrast to the Liouville propagation which is determined by the classical trajectories of the arguments themselves. 1. introduction The Wigner [11] function W is a representation of the quantum density operator ˆρ as a real function on phase space, for a non-relativistic dynamical system with classical-quantum correspondence and n degrees of freedom. Here, all the underlying geometry is euclidean and we can include the Wigner function within the general framework of Weyl representation of operators: Â → A ∈ Ck IR (IR2n), so W(x) = (2π�) −n ρ(x) , where x ≡ (p1,...pn, q1,...qn) denotes a point in phase space. In this context, the quantum evolution of the density operator ˆρt effected by an autonomous external hamiltonian
In quantum mechanics, systems can be described in phase space in terms of the Wigner function and the star-product operation. Quantum characteristics, which appear in the Heisenberg picture as the Weyl's symbols of operators of canonical coordinates and momenta, can be used to solve the evolution equations for symbols of other operators acting in the Hilbert space. To any fixed order in the Planck's constant, many-body potential scattering problem simplifies to a statistical-mechanical problem of computing an ensemble of quantum characteristics and their derivatives with respect to the initial canonical coordinates and momenta. The reduction to a system of ordinary differential equations pertains rigorously at any fixed order in ħ. We present semiclassical expansion of quantum characteristics for many-body scattering problem and provide tools for calculation of average values of time-dependent physical observables and cross sections. The method of quantum characteristics admits the consistent incorporation of specific quantum effects, such as non-locality and coherence in propagation of particles, into the semiclassical transport models. We formulate the principle of stationary action for quantum Hamilton's equations and give quantum-mechanical extensions of the Liouville theorem on conservation of the phase-space volume and the Poincaré theorem on conservation of 2p-forms. The lowest order quantum corrections to the Kepler periodic orbits are constructed. These corrections show the resonance behavior.
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.