The hydrodynamic formulation of mixed quantum states involves a hierarchy of coupled equations of motion for the momentum moments of the Wigner function. In this work a closure scheme for the hierarchy is developed. The closure scheme uses information contained in the lower known moments to expand the Wigner phase-space distribution function in a Gauss-Hermite orthonormal basis. The higher moment required to terminate the hierarchy is then easily obtained from the reconstructed approximate Wigner function by a straightforward integration over the momentum space. Application of the moment closure scheme is demonstrated for the dissipative and nondissipative dynamics of two different systems: (i) double-well potential, (ii) periodic potential.
The hydrodynamic formulation of mixed quantum states involves a hierarchy of coupled equations of motion for the momentum moments of the Wigner function. In this work a closure scheme for the hierarchy is developed. The closure scheme uses information contained in the lower known moments to expand the Wigner phase-space distribution function in a Gauss-Hermite orthonormal basis. The higher moment required to terminate the hierarchy is then easily obtained from the reconstructed approximate Wigner function by a straightforward integration over the momentum space. Application of the moment closure scheme is demonstrated for the dissipative and nondissipative dynamics of two different systems: (i) double-well potential, (ii) periodic potential.
In this paper we formulate the hydrodynamic moment and cumulant expansions of the density matrix for a quantum system in which the time evolution is described by an infinite hierarchy of coupled hydrodynamic moment equations. For the case where the system interacts with a complex environment, we develop a semiclassical closure that truncates the hierarchy. Our numerical approach for solving the truncated system of equations is intimately related to the de Broglie−Bohm interpretation of quantum mechanics, in that the time evolution of each cumulant is constructed over a discrete ensemble of hydrodynamic paths. We discuss the advantages and disadvantages of this scheme and apply the method to a model problem for a particle escaping from a metastable potential well while in contact with a thermal bath.
The theoretical framework of the mixed quantum-classical description given by Burghardt and Parlant [J. Chem. Phys. 120, 3055 (2004)] is detailed. A representation in terms of partial hydrodynamic moments is developed, the dynamics of which is determined by a hierarchy of equations derived from the quantum Liouville equation. Exact equations of motion are obtained, whose quantum-classical approximants are associated with a fluid-dynamical trajectory representation which couples classical variables to quantum hydrodynamic variables. The latter evolve under a generalized hydrodynamic force which also depends upon the classical phase-space variables. The hydrodynamic moment description is shown to be closely connected to mixed quantum-classical phase-space methods.
A newly derive set of quantum fluid dynamical equations appropriate for the description of mixed state dynamics is presented. Based essentially on moments of the Wigner function, the theory presented here uses an expansion of the pseudodistribution function in conjunction with a renormalization technique to obtain a semiclassical approximation for the dynamics of the probability, probability current, and energy densities. Particularly simple equations, ideal fluid dynamics with a nonlocal potential, result from using a local Maxwellian ansatz for the Wigner function. Further, the transformed potential may be easily computed from the fields being convected. Analogies with classical kinetic theory and fluid dynamics are exploited whenever possible.
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.
We treat the non-equilibrium evolution of an open one-particle statistical system, subject to a potential and to an external “heat bath” (hb) with negligible dissipation. For the classical equilibrium Boltzmann distribution, Wc,eq, a non-equilibrium three-term hierarchy for moments fulfills Hermiticity, which allows one to justify an approximate long-time thermalization. That gives partial dynamical support to Boltzmann’s Wc,eq, out of the set of classical stationary distributions, Wc;st, also investigated here, for which neither Hermiticity nor that thermalization hold, in general. For closed classical many-particle systems without hb (by using Wc,eq), the long-time approximate thermalization for three-term hierarchies is justified and yields an approximate Lyapunov function and an arrow of time. The largest part of the work treats an open quantum one-particle system through the non-equilibrium Wigner function, W. Weq for a repulsive finite square well is reported. W’s (< 0 in various cases) are assumed to be quasi-definite functionals regarding their dependences on momentum (q). That yields orthogonal polynomials, HQ,n(q), for Weq (and for stationary Wst), non-equilibrium moments, Wn, of W and hierarchies. For the first excited state of the harmonic oscillator, its stationary Wst is a quasi-definite functional, and the orthogonal polynomials and three-term hierarchy are studied. In general, the non-equilibrium quantum hierarchies (associated with Weq) for the Wn’s are not three-term ones. As an illustration, we outline a non-equilibrium four-term hierarchy and its solution in terms of generalized operator continued fractions. Such structures also allow one to formulate long-time approximations, but make it more difficult to justify thermalization. For large thermal and de Broglie wavelengths, the dominant Weq and a non-equilibrium equation for W are reported: the non-equilibrium hierarchy could plausibly be a three-term one and possibly not far from Gaussian, and thermalization could possibly be justified.