Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2512.11174

Negative Marginal Densities in Mixed Quantum-Classical Liouville Dynamics

Kai Gu, Jeremy Schofield

quant-phphysics.chem-phSubmitted Dec 11, 2025• Announced Dec 15, 2025• 0 citations
arXivHTMLPDF

Abstract

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.

Comments
53 pages, 13 figures

Related papers

Negative marginal densities in mixed quantum–classical Liouville dynamics

Kai Gu, Jeremy Schofield

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-classical descriptions.

The Journal of Chemical Physics•Feb 19, 2026•0 citations

Mixed quantum-classical dynamics

Raymond Kapral, Giovanni Ciccotti

Mixed quantum-classical equations of motion are derived for a quantum subsystem of light (mass m) particles coupled to a classical bath of massive (mass M) particles. The equation of motion follows from a partial Wigner transform over the bath degrees of freedom of the Liouville equation for the full quantum system, followed by an expansion in the small parameter μ=(m/M)1/2 in analogy with the theory of Brownian motion. The resulting mixed quantum-classical Liouville equation accounts for the coupled evolution of the subsystem and bath. The quantum subsystem is represented in an adiabatic (or other) basis and the series solution of the Liouville equation leads to a representation of the dynamics in an ensemble of surface-hopping trajectories. A generalized Pauli master equation for the evolution of the diagonal elements of the density matrix is derived by projection operator methods and its structure is analyzed in terms of surface-hopping trajectories.

The Journal of Chemical Physics•May 8, 1999•677 citations

Partial Wigner transforms and the quantum-classical Liouville equation

Christof Schütte

In molecular dynamics applications there is a growing interest in mixed quantum-classical models. The quantum-classical Liouville equation (QCL) describes most atoms of the molecular system under consideration by means of a classical phase space density but an important, small portion of the system by means of quantum mechanics. The QCL is derived from the full quantum dynamical (QD) description by applying the Wigner transform to the \\classical part" of the system only. We discuss the conditions under which the QCL model approximates the full QD evolution of the system. First, analysis of the asymptotic properties of the Wigner transform shows that solving the QCL yields a rst order approximation of full quantum dynamics. Second, we discuss the adiabatic limit of the QCL. This discussion shows that the QCL solutions may be interpretated as classical phase space densities, at least near the adiabatic limit. Third, it is demonstrated that the QCL yields good approximations of non-adia...

•Apr 12, 1999•2 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.

arXiv:1707.09660•Jul 30, 2017•1 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.

arXiv:cond-mat/0612203•The Journal of Chemical Physics•Mar 20, 2007•23 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.

arXiv:1812.01064•The Journal of Chemical Physics•Dec 11, 2018•15 citations
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-classical descriptions.
Mixed quantum-classical equations of motion are derived for a quantum subsystem of light (mass m) particles coupled to a classical bath of massive (mass M) particles. The equation of motion follows from a partial Wigner transform over the bath degrees of freedom of the Liouville equation for the full quantum system, followed by an expansion in the small parameter μ=(m/M)1/2 in analogy with the theory of Brownian motion. The resulting mixed quantum-classical Liouville equation accounts for the coupled evolution of the subsystem and bath. The quantum subsystem is represented in an adiabatic (or other) basis and the series solution of the Liouville equation leads to a representation of the dynamics in an ensemble of surface-hopping trajectories. A generalized Pauli master equation for the evolution of the diagonal elements of the density matrix is derived by projection operator methods and its structure is analyzed in terms of surface-hopping trajectories.
In molecular dynamics applications there is a growing interest in mixed quantum-classical models. The quantum-classical Liouville equation (QCL) describes most atoms of the molecular system under consideration by means of a classical phase space density but an important, small portion of the system by means of quantum mechanics. The QCL is derived from the full quantum dynamical (QD) description by applying the Wigner transform to the \\classical part" of the system only. We discuss the conditions under which the QCL model approximates the full QD evolution of the system. First, analysis of the asymptotic properties of the Wigner transform shows that solving the QCL yields a rst order approximation of full quantum dynamics. Second, we discuss the adiabatic limit of the QCL. This discussion shows that the QCL solutions may be interpretated as classical phase space densities, at least near the adiabatic limit. Third, it is demonstrated that the QCL yields good approximations of non-adia...
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.
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 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.