The mixed quantum-classical Liouville equation (QCLE) provides an approximate perturbative framework for describing the dynamics of systems with coupled quantum and classical degrees of freedom of disparate thermal wavelengths. The evolution governed by the Liouville operator preserves many properties of full quantum dynamics, including the conservation of total population, energy, and purity, and has shown quantitative agreement with exact quantum results for the expectation values of many observables where direct comparisons are feasible. However, since the QCLE density matrix operator is obtained from the partial Wigner transform of the full quantum density matrix, its matrix elements can have negative values, implying that the diagonal matrix elements behave as pseudo-densities rather than densities of classical phase space. Here, we compare phase-space distributions generated by exact quantum dynamics with those produced by QCLE evolution from pure quantum initial states. We show that resonance effects in the off-diagonal matrix elements differ qualitatively, particularly for low-energy states. Furthermore, numerical and analytical results for low-dimensional models reveal that the QCLE can violate the positivity of marginal phase-space densities, a property that should hold at all times for any physical system. A perturbative analysis of a model system confirms that such violations arise generically. We also show that the violations of positivity of the marginal densities vanish as the initial energy of the system increases relative to the energy gap between subsystem states. These findings suggest that a negativity index, quantifying deviations from positivity, may provide a useful metric for assessing the validity of mixed quantum\textendash{}classical descriptions.The mixed quantum-classical Liouville equation (QCLE) provides an approximate perturbative framework for describing the dynamics of systems with coupled quantum and classical degrees of freedom of disparate thermal wavelengths. The evolution governed by the Liouville operator preserves many properties of full quantum dynamics, including the conservation of total population, energy, and purity, and has shown quantitative agreement with exact quantum results for the expectation values of many observables where direct comparisons are feasible. However, since the QCLE density matrix operator is obtained from the partial Wigner transform of the full quantum density matrix, its matrix elements can have negative values, implying that the diagonal matrix elements behave as pseudo-densities rather than densities of classical phase space. Here, we compare phase-space distributions generated by exact quantum dynamics with those produced by QCLE evolution from pure quantum initial states. We show that resonance effects in the off-diagonal matrix elements differ qualitatively, particularly for low-energy states. Furthermore, numerical and analytical results for low-dimensional models reveal that the QCLE can violate the positivity of marginal phase-space densities, a property that should hold at all times for any physical system. A perturbative analysis of a model system confirms that such violations arise generically. We also show that the violations of positivity of the marginal densities vanish as the initial energy of the system increases relative to the energy gap between subsystem states. These findings suggest that a negativity index, quantifying deviations from positivity, may provide a useful metric for assessing the validity of mixed quantum\textendash{}classical descriptions.

quant-pharXiv:2512.11174Dec 15, 20250 citations

On Generated Dynamics for Open Quantum Systems: Spectral Analysis of Effective Liouville

Martin Janßen

We point out that the quantum dynamical map of an open quantum system can be generated by an effective Liouville operator. The effective Liouville shows the dynamical breaking of time reversibility. This breaking of reversibility is expressed by the effective Liouville's discrete spectrum having negative imaginary parts. This generated dynamics of open quantum systems is capable of memory effects described by a frequency dependence of the spectrum. When memory effects can be neglected or smoothed out, the effective Liouville generates the well known semi-group dynamics of open systems in the Markov approximation. The spectral analysis of the effective Liouville -with or without memory- allows to represent the quantum dynamical map and expectation values of physical quantities by metric expressions using complex eigenvalues and bi-orthonormal eigenmodes. The long time dynamics is dominated by the zero-mode, which forms the stationary state of the open system. The relaxation time scale $τ$ is set by the inverse of the smallest negative imaginary part of the non-vanishing effective eigenvalues of the effective Liouville. In generic cases of system-environment-coupling this time scale $τ$ is the relevant scale for both, the relaxation of diagonal elements (populations) and the decay of off-diagonal elements (coherences) of an initial density matrix in the stationary state's diagonal representation. It also sets the time scale after which the entropy tends to its stationary value. The effective Liouville and its spectral content thus give a framework to study the dynamic breaking of time reversibility, memory effects, decoherence, relaxation to a stationary state and the role of entropy in open quantum systems.We point out that the quantum dynamical map of an open quantum system can be generated by an effective Liouville operator. The effective Liouville shows the dynamical breaking of time reversibility. This breaking of reversibility is expressed by the effective Liouville's discrete spectrum having negative imaginary parts. This generated dynamics of open quantum systems is capable of memory effects described by a frequency dependence of the spectrum. When memory effects can be neglected or smoothed out, the effective Liouville generates the well known semi-group dynamics of open systems in the Markov approximation. The spectral analysis of the effective Liouville -with or without memory- allows to represent the quantum dynamical map and expectation values of physical quantities by metric expressions using complex eigenvalues and bi-orthonormal eigenmodes. The long time dynamics is dominated by the zero-mode, which forms the stationary state of the open system. The relaxation time scale $τ$ is set by the inverse of the smallest negative imaginary part of the non-vanishing effective eigenvalues of the effective Liouville. In generic cases of system-environment-coupling this time scale $τ$ is the relevant scale for both, the relaxation of diagonal elements (populations) and the decay of off-diagonal elements (coherences) of an initial density matrix in the stationary state's diagonal representation. It also sets the time scale after which the entropy tends to its stationary value. The effective Liouville and its spectral content thus give a framework to study the dynamic breaking of time reversibility, memory effects, decoherence, relaxation to a stationary state and the role of entropy in open quantum systems.

arXiv:1707.09660Jul 30, 20171 citation

Decoherence and quantum-classical master equation dynamics

Robbie Grunwald, Raymond Kapral

The conditions under which quantum-classical Liouville dynamics may be reduced to a master equation are investigated. Systems that can be partitioned into a quantum-classical subsystem interacting with a classical bath are considered. Starting with an exact non-Markovian equation for the diagonal elements of the density matrix, an evolution equation for the subsystem density matrix is derived. One contribution to this equation contains the bath average of a memory kernel that accounts for all coherences in the system. It is shown to be a rapidly decaying function, motivating a Markovian approximation on this term in the evolution equation. The resulting subsystem density matrix equation is still non-Markovian due to the fact that bath degrees of freedom have been projected out of the dynamics. Provided the computation of nonequilibrium average values or correlation functions is considered, the non-Markovian character of this equation can be removed by lifting the equation into the full phase space of the system. This leads to a trajectory description of the dynamics where each fictitious trajectory accounts for decoherence due to the bath degrees of freedom. The results are illustrated by computations of the rate constant of a model nonadiabatic chemical reaction.The conditions under which quantum-classical Liouville dynamics may be reduced to a master equation are investigated. Systems that can be partitioned into a quantum-classical subsystem interacting with a classical bath are considered. Starting with an exact non-Markovian equation for the diagonal elements of the density matrix, an evolution equation for the subsystem density matrix is derived. One contribution to this equation contains the bath average of a memory kernel that accounts for all coherences in the system. It is shown to be a rapidly decaying function, motivating a Markovian approximation on this term in the evolution equation. The resulting subsystem density matrix equation is still non-Markovian due to the fact that bath degrees of freedom have been projected out of the dynamics. Provided the computation of nonequilibrium average values or correlation functions is considered, the non-Markovian character of this equation can be removed by lifting the equation into the full phase space of the system. This leads to a trajectory description of the dynamics where each fictitious trajectory accounts for decoherence due to the bath degrees of freedom. The results are illustrated by computations of the rate constant of a model nonadiabatic chemical reaction.

arXiv:cond-mat/0612203The Journal of Chemical PhysicsMar 20, 200723 citations

Heat transfer statistics in mixed quantum-classical systems

Junjie Liu, Chang‐Yu Hsieh, Dvira Segal +1 more

The modelling of quantum heat transfer processes at the nanoscale is crucial for the development of energy harvesting and molecular electronic devices. Herein, we adopt a mixed quantum-classical description of a device, in which the open subsystem of interest is treated quantum mechanically and the surrounding heat baths are treated in a classical-like fashion. By introducing such a mixed quantum-classical description of the composite system, one is able to study the heat transfer between the subsystem and bath from a closed system point of view, thereby avoiding simplifying assumptions related to the bath time scale and subsystem-bath coupling strength. In particular, we adopt the full counting statistics approach to derive a general expression for the moment generating function of heat in systems whose dynamics are described by the quantum-classical Liouville equation (QCLE). From this expression, one can deduce expressions for the dynamics of the average heat and heat current, which may be evaluated using numerical simulations. Due to the approximate nature of the QCLE, we also find that the steady state fluctuation symmetry holds up to order ℏ for systems whose subsystem-bath couplings and baths go beyond bilinear and harmonic, respectively. To demonstrate the approach, we consider the nonequilibrium spin boson model and simulate its time-dependent average heat and heat current under various conditions.The modelling of quantum heat transfer processes at the nanoscale is crucial for the development of energy harvesting and molecular electronic devices. Herein, we adopt a mixed quantum-classical description of a device, in which the open subsystem of interest is treated quantum mechanically and the surrounding heat baths are treated in a classical-like fashion. By introducing such a mixed quantum-classical description of the composite system, one is able to study the heat transfer between the subsystem and bath from a closed system point of view, thereby avoiding simplifying assumptions related to the bath time scale and subsystem-bath coupling strength. In particular, we adopt the full counting statistics approach to derive a general expression for the moment generating function of heat in systems whose dynamics are described by the quantum-classical Liouville equation (QCLE). From this expression, one can deduce expressions for the dynamics of the average heat and heat current, which may be evaluated using numerical simulations. Due to the approximate nature of the QCLE, we also find that the steady state fluctuation symmetry holds up to order ℏ for systems whose subsystem-bath couplings and baths go beyond bilinear and harmonic, respectively. To demonstrate the approach, we consider the nonequilibrium spin boson model and simulate its time-dependent average heat and heat current under various conditions.

arXiv:1812.01064The Journal of Chemical PhysicsDec 11, 201815 citations

A Study of Quantum-classical Dynamics in the Mapping Basis

Ali M. Nassimi

Solving quantum dynamics is an exponentially difficult problem. Thus, an exact numerical solution is inaccessible for any condensed matter system.\nA promising approach is to divide the system into a quantum subsystem containing degrees of freedom which are of greater interest or those which have more profound quantum character (e.g., have smaller mass) and a classical bath containing the rest of the system.\nImposing such a partition and treating the bath classically results in quantum-classical dynamics. The quantum-classical Liouville equation is a general equation in the Hilbert space of quantum degrees of freedom while it resides in the phase space of the classical degrees of freedom.\nAny numerical solution to this equation requires representation of the quantum subsystem in some basis. Solutions to this equation have been already proposed in the subsystem, adiabatic and force bases, each with its own cons and pros.\n\nIn this work, the quantum-classical equations of motion are cast in the subsystem basis and subsequently mapped to a number of fictitious harmonic oscillators.\nThe result is quantum-classical dynamics in the mapping basis which treats both quantum and classical degrees of freedom on the same footing, i.e., in phase space. Neglecting a portion of the back reaction of the quantum-subsystem to classical bath results in an expression for the time evolution of an operator (density matrix) equal to its Poisson bracket with the Hamiltonian.\nThis equation can be solved in terms of characteristics to provide a computationally tractable method for calculating quantum-classical dynamical properties. The expressions for expectation values and correlation functions in this formalism are derived.\nCalculations on spin-boson system, barrier crossing models---the so called Tully models---and the Fenna-Mathews-Olson pigments show very good agreement between the results of this method and numerical solutions to the Schrödinger equation.Solving quantum dynamics is an exponentially difficult problem. Thus, an exact numerical solution is inaccessible for any condensed matter system.\nA promising approach is to divide the system into a quantum subsystem containing degrees of freedom which are of greater interest or those which have more profound quantum character (e.g., have smaller mass) and a classical bath containing the rest of the system.\nImposing such a partition and treating the bath classically results in quantum-classical dynamics. The quantum-classical Liouville equation is a general equation in the Hilbert space of quantum degrees of freedom while it resides in the phase space of the classical degrees of freedom.\nAny numerical solution to this equation requires representation of the quantum subsystem in some basis. Solutions to this equation have been already proposed in the subsystem, adiabatic and force bases, each with its own cons and pros.\n\nIn this work, the quantum-classical equations of motion are cast in the subsystem basis and subsequently mapped to a number of fictitious harmonic oscillators.\nThe result is quantum-classical dynamics in the mapping basis which treats both quantum and classical degrees of freedom on the same footing, i.e., in phase space. Neglecting a portion of the back reaction of the quantum-subsystem to classical bath results in an expression for the time evolution of an operator (density matrix) equal to its Poisson bracket with the Hamiltonian.\nThis equation can be solved in terms of characteristics to provide a computationally tractable method for calculating quantum-classical dynamical properties. The expressions for expectation values and correlation functions in this formalism are derived.\nCalculations on spin-boson system, barrier crossing models---the so called Tully models---and the Fenna-Mathews-Olson pigments show very good agreement between the results of this method and numerical solutions to the Schrödinger equation.

Aug 31, 20110 citations

The Quantum Mechanical Density Operator And Its Time Evolution: Quantum Dynamics Using The Quantum Liouville Equation

Abraham Nitzan

The starting point of the classical description of motion is the Newton equations that yield a phase space trajectory (rN (t), pN (t)) for a given initial condition (rN (0), pN (0)). Alternatively one may describe classical motion in the framework of the Liouville equation (Section (1.2.2)) that describes the time evolution of the phase space probability density f (rN , pN ; t). For a closed system fully described in terms of a well specified initial condition, the two descriptions are completely equivalent. Probabilistic treatment becomes essential in reduced descriptions that focus on parts of an overall system, as was demonstrated in Sections 5.1–5.3 for equilibrium systems, and in Chapters 7 and 8 that focus on the time evolution of classical systems that interact with their thermal environments. This chapter deals with the analogous quantum mechanical problem. Within the limitations imposed by its nature as expressed, for example, by Heisenbergtype uncertainty principles, the Schrödinger equation is deterministic. Obviously it describes a deterministic evolution of the quantum mechanical wavefunction. The analog of the phase space probability density f (rN , pN ; t) is now the quantum mechanical density operator (often referred to as the “density matrix”), whose time evolution is determined by the quantum Liouville equation. Again, when the system is fully described in terms of a well specified initial wavefunction, the two descriptions are equivalent. The density operator formalism can, however, be carried over to situations where the initial state of the system is not well characterized and/or a reduced description of part of the overall system is desired. Such situations are considered later in this chapter.The starting point of the classical description of motion is the Newton equations that yield a phase space trajectory (rN (t), pN (t)) for a given initial condition (rN (0), pN (0)). Alternatively one may describe classical motion in the framework of the Liouville equation (Section (1.2.2)) that describes the time evolution of the phase space probability density f (rN , pN ; t). For a closed system fully described in terms of a well specified initial condition, the two descriptions are completely equivalent. Probabilistic treatment becomes essential in reduced descriptions that focus on parts of an overall system, as was demonstrated in Sections 5.1–5.3 for equilibrium systems, and in Chapters 7 and 8 that focus on the time evolution of classical systems that interact with their thermal environments. This chapter deals with the analogous quantum mechanical problem. Within the limitations imposed by its nature as expressed, for example, by Heisenbergtype uncertainty principles, the Schrödinger equation is deterministic. Obviously it describes a deterministic evolution of the quantum mechanical wavefunction. The analog of the phase space probability density f (rN , pN ; t) is now the quantum mechanical density operator (often referred to as the “density matrix”), whose time evolution is determined by the quantum Liouville equation. Again, when the system is fully described in terms of a well specified initial wavefunction, the two descriptions are equivalent. The density operator formalism can, however, be carried over to situations where the initial state of the system is not well characterized and/or a reduced description of part of the overall system is desired. Such situations are considered later in this chapter.

Oxford University Press eBooksApr 6, 20061 citation

τ

is set by the inverse of the smallest negative imaginary part of the non-vanishing effective eigenvalues of the effective Liouville. In generic cases of system-environment-coupling this time scale

τ

is the relevant scale for both, the relaxation of diagonal elements (populations) and the decay of off-diagonal elements (coherences) of an initial density matrix in the stationary state's diagonal representation. It also sets the time scale after which the entropy tends to its stationary value. The effective Liouville and its spectral content thus give a framework to study the dynamic breaking of time reversibility, memory effects, decoherence, relaxation to a stationary state and the role of entropy in open quantum systems.