Negative Marginal Densities in Mixed Quantum-Classical Liouville Dynamics

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.

• The paper demonstrates that QCLE dynamics can evolve physically valid states into nonphysical negative marginal densities, raising concerns about classical probability violations.
• It employs Tully’s dual avoided crossing model to quantify negativity using a 'negativity index' and highlights energy-dependent deviations.
• The analysis contrasts local derivative approximations in QCLE with full quantum evolution, outlining regimes where quantum-classical methods may fail.

Negative Marginal Densities in Mixed Quantum-Classical Liouville Dynamics

Overview

The study addresses a fundamental limitation of the mixed quantum-classical Liouville equation (QCLE) in the context of nonadiabatic quantum dynamics. QCLE provides a core framework for mixed quantum-classical systems, widely applied in simulating nonadiabatic processes where quantum subsystems interact with a classical environment, such as electron/proton transfer and energy transport. The paper demonstrates that the QCLE can violate positive semidefiniteness of marginal phase-space densities, both numerically and analytically, in low-dimensional models—contradicting the physical requirement of positivity for classical distributions and raising questions about what constitutes a physically consistent quantum-classical description.

QCLE Formulation and Conservation Laws

QCLE emerges from a systematic expansion of the quantum Liouville-von Neumann equation via a partial Wigner transform, separating quantum and classical degrees of freedom based on thermal de Broglie wavelengths. Unlike heuristic models, the QCLE preserves probability, energy, and purity due to its anti-Hermitian Liouvillian structure. In the classical limit, QCLE reduces appropriately, yet it captures quantum coherence and nonadiabatic transitions—making it a principal model for benchmarking more approximate mixed-quantum-classical dynamical methods. Formally, QCLE is

$$\frac{\partial \hat{\rho}_W}{\partial t} = -\frac{i}{\hbar}[\hat{H}_W, \hat{\rho}_W] + \frac{1}{2}\left( \{ \hat{H}_W, \hat{\rho}_W \} - \{ \hat{\rho}_W, \hat{H}_W \} \right)$$

where the first term represents quantum commutator evolution and the second embodies leading-order quantum-classical coupling corrections.

Marginal Densities and their Positivity

Marginal densities for classical coordinates, defined as integrals of the pseudo-density over complementary degrees of freedom, should reflect observable position or momentum probability distributions. In quantum mechanics, marginals derived from the density operator are always positive. However, because the QCLE employs the partial Wigner transform, its diagonal components can attain negative values, and thus the classical “pseudo-density” $\rho(R,P,t)$ is not pointwise positive.

While the negativity of the full phase-space quasi-probability distribution can be interpreted as a quantum effect (analogous to the full Wigner function in quantum mechanics), physical marginal densities such as $$n(R,t) = \int \mathrm{Tr}\, \hat{\rho}_W(R,P,t)\, dP$$ and $$\eta(P,t) = \int \mathrm{Tr}\, \hat{\rho}_W(R,P,t)\, dR$$ should remain non-negative to be physically interpretable as probabilities. The paper rigorously demonstrates that this expectation fails in QCLE dynamics.

Numerical Demonstration: Tully's Dual Avoided Crossing Model

The central numerical result is based on Tully’s dual avoided crossing (DAC) model—a canonical two-level, one-dimensional system that allows direct comparison between exact quantum and QCLE trajectories. Initial pure quantum states are propagated according to both schemes. While QCLE reproduces final populations and mean observables with high accuracy—even at low energy—maps of phase-space pseudo-densities and especially their marginals show pronounced negative regions at low and intermediate energies. These negative regions persist even when the QCLE accurately predicts population transfer and expectation values.

A quantitative metric, the negativity index $\mathcal{N}[f]$, is used to assess the extent of positivity violation in marginals. This index is nonzero for QCLE-generated marginals below a critical initial bath energy; it decays rapidly as the energy increases relative to the system’s energy gap.

Analytical Results: Constant Nonadiabatic Coupling

An analytically tractable force-free model with constant diabatic and adiabatic couplings is constructed to elucidate the origin of marginal negativity. By comparing QCLE to the full quantum Liouville equation, it is evident that QCLE replaces nonlocal (momentum-shifted) terms with local derivatives. This local approximation is accurate only in the limit where the nonadiabatic coupling is weak relative to the wavepacket dispersion, or the system is far from resonance. When these conditions are not met, incorrect negative values generically develop in the marginal distributions.

Perturbative solutions for the marginal densities confirm that nonphysical negative regions occur precisely where the QCLE’s locality approximation deviates significantly from the exact quantum propagation, especially as higher-order terms in the expansion become non-negligible.

Implications for Mixed Quantum-Classical Simulation

The identification of negative marginal densities under QCLE evolution, despite preservation of average observables and overall probability, demonstrates an intrinsic failure distinct from previously known pathologies in open quantum system dynamics (such as those encountered with Redfield or Lindblad-type master equations). Importantly, the negativity does not originate from Markovianity or irreversible approximations, but from the local truncation of phase-space dynamics.

For high-dimensional, condensed phase systems with large classical baths, the fractional volume with negative pseudo-density is expected to become exponentially small, meaning these violations are likely negligible in practice. However, for low-dimensional or strongly coupled systems, violation of marginal positivity is a diagnostic of the breakdown of quantum-classical separation and indicates regimes where QCLE-based surface-hopping, mapping, or related approximate methods may yield unphysical results. The proposed negativity index provides a tool to assess and benchmark the validity of QCLE-based simulations in such cases.

Conclusion

The study robustly establishes that the QCLE can evolve physically valid initial quantum states into nonphysical pseudo-densities with negative marginal distributions for classical degrees of freedom. This occurs generically for low-dimensional, resonant, or strongly nonadiabatic models and persists until the initial energy dominates all relevant gap scales. Theoretical analysis links this violation to the locality approximation inherent to the QCLE. For high-dimensional systems or high energy initial states, negativity rapidly decays and the QCLE is consistent with quantum-classical expectations. The work motivates use of marginal-based negativity diagnostics in the assessment of hybrid quantum-classical methods and establishes a clear boundary for the regime of validity of the QCLE.

Reference:

"Negative Marginal Densities in Mixed Quantum-Classical Liouville Dynamics" (2512.11174)