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:2610.16864

Exact Gaussian--Hermite closure of the mapping quantum--classical Liouville equation

Sahand Nikzat Ali Nassimi, Rony Modak Parai

physics.chem-phphysics.comp-phquant-phSubmitted Oct 12, 2026• Announced Oct 13, 2026• 0 citations• Revised Oct 15, 2026 (v2)
v1 · Oct 12, 2026v2 · Oct 15, 2026
arXivHTMLPDF

Abstract

In the Meyer--Miller--Stock--Thoss mapping representation, the Liouville operator of the quantum--classical Liouville equation contains an excess coupling term in addition to the Poisson bracket. Neglecting this term yields the Poisson bracket mapping equation (PBME). We expand the mapping dependence of the density in Hermite functions about the Gaussian factor common to all physical mapping densities. In this basis, the excess coupling term acts purely as a raising operator and exactly cancels the raising contribution generated by the back--reaction force. Consequently, the exact dynamics does not generate higher Hermite orders and closes on a finite, sparse system of equations for the coefficient fields. Within the physical sector, these equations reduce to the subsystem quantum--classical Liouville equation. In contrast, the PBME retains the raising contribution of the force, modifying the back reaction within the physical sector while continuously generating higher Hermite orders. For Tully's dual avoided crossing model, direct propagation of the coefficient equations confirms this finite closure and agrees with independent solutions of the subsystem equation. Comparison with converged full quantum dynamics further shows that, across all regimes examined, the PBME is less accurate than the quantum--classical description for populations, coherences, and momentum distributions.

Related papers

Analysis of the quantum-classical Liouville equation in the mapping basis

Ali M. Nassimi, Sara Bonella, Raymond Kapral

The quantum-classical Liouville equation provides a description of the dynamics of a quantum subsystem coupled to a classical environment. Representing this equation in the mapping basis leads to a continuous description of discrete quantum states of the subsystem and may provide an alternate route to the construction of simulation schemes. In the mapping basis the quantum-classical Liouville equation consists of a Poisson bracket contribution and a more complex term. By transforming the evolution equation, term-by-term, back to the subsystem basis, the complex term (excess coupling term) is identified as being due to a fraction of the back reaction of the quantum subsystem on its environment. A simple approximation to quantum-classical Liouville dynamics in the mapping basis is obtained by retaining only the Poisson bracket contribution. This approximate mapping form of the quantum-classical Liouville equation can be simulated easily by Newtonian trajectories. We provide an analysis of the effects of neglecting the presence of the excess coupling term on the expectation values of various types of observables. Calculations are carried out on nonadiabatic population and quantum coherence dynamics for curve crossing models. For these observables, the effects of the excess coupling term enter indirectly in the computation and good estimates are obtained with the simplified propagation.

arXiv:1006.2075•The Journal of Chemical Physics•Oct 6, 2010•71 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.

•Aug 31, 2011•0 citations

Mapping quantum-classical Liouville equation: Projectors and trajectories

Aaron Kelly, Ramses van Zon, Jeremy Schofield +1 more

The evolution of a mixed quantum-classical system is expressed in the mapping formalism where discrete quantum states are mapped onto oscillator states, resulting in a phase space description of the quantum degrees of freedom. By defining projection operators onto the mapping states corresponding to the physical quantum states, it is shown that the mapping quantum-classical Liouville operator commutes with the projection operator so that the dynamics is confined to the physical space. It is also shown that a trajectory-based solution of this equation can be constructed that requires the simulation of an ensemble of entangled trajectories. An approximation to this evolution equation which retains only the Poisson bracket contribution to the evolution operator does admit a solution in an ensemble of independent trajectories but it is shown that this operator does not commute with the projection operators and the dynamics may take the system outside the physical space. The dynamical instabilities, utility, and domain of validity of this approximate dynamics are discussed. The effects are illustrated by simulations on several quantum systems.

arXiv:1201.1042•The Journal of Chemical Physics•Feb 22, 2012•121 citations

Hermitization and the Poisson bracket-commutator correspondence as a consequence of averaging

Adriana I. Pesci, Raymond E. Goldstein, Hermann Uys

Here we present a study of the solutions and mathematical structure of the momentum-averaged Liouville (or Collisionless Boltzmann) equation in Fourier space. We show that the averaging procedure leads to a formalism identical to that of the density matrix of quantum mechanics. This mathematical mapping leaves the averages of all quantities unaltered and provides a unique way to construct the Hermitian version of a given operator. This seems to be the only method that resolves the ambiguity of Hermitization of operators that contain products of non-commuting variables. We also present a systematic perturbation scheme to evaluate correctly the classical solutions from the quantum ones and a formal proof of the approximate correspondence between the Poisson brackets and commutators.

Journal of Physics A Mathematical and General•Jan 11, 2006•1 citation

Quantum-classical Liouville dynamics in the mapping basis

Hyojoon Kim, Ali Nassimi, Raymond Kapral

The quantum-classical Liouville equation describes the dynamics of a quantum subsystem coupled to a classical environment. It has been simulated using various methods, notably, surface-hopping schemes. A representation of this equation in the mapping Hamiltonian basis for the quantum subsystem is derived. The resulting equation of motion, in conjunction with expressions for quantum expectation values in the mapping basis, provides another route to the computation of the nonadiabatic dynamics of observables that does not involve surface-hopping dynamics. The quantum-classical Liouville equation is exact for the spin-boson system. This well-known model is simulated using an approximation to the evolution equation in the mapping basis, and close agreement with exact quantum results is found.

arXiv:0903.0847•The Journal of Chemical Physics•Aug 25, 2008•156 citations

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\textendash{}classical descriptions.

quant-pharXiv:2512.11174•Dec 15, 2025•0 citations
In the Meyer--Miller--Stock--Thoss mapping representation, the Liouville operator of the quantum--classical Liouville equation contains an excess coupling term in addition to the Poisson bracket. Neglecting this term yields the Poisson bracket mapping equation (PBME). We expand the mapping dependence of the density in Hermite functions about the Gaussian factor common to all physical mapping densities. In this basis, the excess coupling term acts purely as a raising operator and exactly cancels the raising contribution generated by the back--reaction force. Consequently, the exact dynamics does not generate higher Hermite orders and closes on a finite, sparse system of equations for the coefficient fields. Within the physical sector, these equations reduce to the subsystem quantum--classical Liouville equation. In contrast, the PBME retains the raising contribution of the force, modifying the back reaction within the physical sector while continuously generating higher Hermite orders. For Tully's dual avoided crossing model, direct propagation of the coefficient equations confirms this finite closure and agrees with independent solutions of the subsystem equation. Comparison with converged full quantum dynamics further shows that, across all regimes examined, the PBME is less accurate than the quantum--classical description for populations, coherences, and momentum distributions.
The quantum-classical Liouville equation provides a description of the dynamics of a quantum subsystem coupled to a classical environment. Representing this equation in the mapping basis leads to a continuous description of discrete quantum states of the subsystem and may provide an alternate route to the construction of simulation schemes. In the mapping basis the quantum-classical Liouville equation consists of a Poisson bracket contribution and a more complex term. By transforming the evolution equation, term-by-term, back to the subsystem basis, the complex term (excess coupling term) is identified as being due to a fraction of the back reaction of the quantum subsystem on its environment. A simple approximation to quantum-classical Liouville dynamics in the mapping basis is obtained by retaining only the Poisson bracket contribution. This approximate mapping form of the quantum-classical Liouville equation can be simulated easily by Newtonian trajectories. We provide an analysis of the effects of neglecting the presence of the excess coupling term on the expectation values of various types of observables. Calculations are carried out on nonadiabatic population and quantum coherence dynamics for curve crossing models. For these observables, the effects of the excess coupling term enter indirectly in the computation and good estimates are obtained with the simplified propagation.
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.
The evolution of a mixed quantum-classical system is expressed in the mapping formalism where discrete quantum states are mapped onto oscillator states, resulting in a phase space description of the quantum degrees of freedom. By defining projection operators onto the mapping states corresponding to the physical quantum states, it is shown that the mapping quantum-classical Liouville operator commutes with the projection operator so that the dynamics is confined to the physical space. It is also shown that a trajectory-based solution of this equation can be constructed that requires the simulation of an ensemble of entangled trajectories. An approximation to this evolution equation which retains only the Poisson bracket contribution to the evolution operator does admit a solution in an ensemble of independent trajectories but it is shown that this operator does not commute with the projection operators and the dynamics may take the system outside the physical space. The dynamical instabilities, utility, and domain of validity of this approximate dynamics are discussed. The effects are illustrated by simulations on several quantum systems.
Here we present a study of the solutions and mathematical structure of the momentum-averaged Liouville (or Collisionless Boltzmann) equation in Fourier space. We show that the averaging procedure leads to a formalism identical to that of the density matrix of quantum mechanics. This mathematical mapping leaves the averages of all quantities unaltered and provides a unique way to construct the Hermitian version of a given operator. This seems to be the only method that resolves the ambiguity of Hermitization of operators that contain products of non-commuting variables. We also present a systematic perturbation scheme to evaluate correctly the classical solutions from the quantum ones and a formal proof of the approximate correspondence between the Poisson brackets and commutators.
The quantum-classical Liouville equation describes the dynamics of a quantum subsystem coupled to a classical environment. It has been simulated using various methods, notably, surface-hopping schemes. A representation of this equation in the mapping Hamiltonian basis for the quantum subsystem is derived. The resulting equation of motion, in conjunction with expressions for quantum expectation values in the mapping basis, provides another route to the computation of the nonadiabatic dynamics of observables that does not involve surface-hopping dynamics. The quantum-classical Liouville equation is exact for the spin-boson system. This well-known model is simulated using an approximation to the evolution equation in the mapping basis, and close agreement with exact quantum results is found.
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.