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:1201.1042 · The Journal of Chemical Physics

Mapping quantum-classical Liouville equation: Projectors and trajectories

Aaron Kelly, Ramses van Zon, Jeremy Schofield, Raymond Kapral

Published Feb 22, 2012• 121 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1063/1.3685420

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

Nonadiabatic dynamics in open quantum-classical systems: Forward-backward trajectory solution

Chang‐Yu Hsieh, Raymond Kapral

A new approximate solution to the quantum-classical Liouville equation is derived starting from the formal solution of this equation in forward-backward form. The time evolution of a mixed quantum-classical system described by this equation is obtained in a coherent state basis using the mapping representation, which expresses N quantum degrees of freedom in a 2N-dimensional phase space. The solution yields a simple dynamics in which a set of N coherent state coordinates evolves in forward and backward trajectories, while the bath coordinates evolve under the influence of the mean potential that depends on these forward and backward trajectories. It is shown that the solution satisfies the differential form of the quantum-classical Liouville equation exactly. Relations to other mixed quantum-classical and semi-classical schemes are discussed.

arXiv:1205.1283•The Journal of Chemical Physics•Jul 18, 2012•112 citations

Nonadiabatic dynamics in open quantum-classical systems: Forward-backward trajectory solution

Chang‐Yu Hsieh, Raymond Kapral

A new approximate solution to the quantum-classical Liouville equation is derived starting from the formal solution of this equation in forward-backward form. The time evolution of a mixed quantum-classical system described by this equation is obtained in a coherent state basis using the mapping representation, which expresses N quantum degrees of freedom in a 2N-dimensional phase space. The solution yields a simple non-Hamiltonian dynamics in which a set of N coherent state coordinates evolve in forward and backward trajectories while the bath coordinates evolve under the influence of the mean potential that depends on these forward and backward tra-jectories. It is shown that the solution satisfies the differential form of the quantum-classical Liouville equation exactly. Relations to other mixed quantum-classical and semi-classical schemes are discussed. 1 ar

•Aug 17, 2016•70 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

Dynamics of quantum–classical systems in nonequilibrium environments

Jeremy Schofield, Raymond Kapral

The dynamics of a quantum system coupled to a classical environment and subject to constraints that drive it out of equilibrium are described. The evolution of the system is governed by the quantum-classical Liouville equation. Rather than evaluating the evolution of the mixed quantum-classical density operator, we derive exact equations of motion for the nonequilibrium average values of a set of operators or variables, along with correlation function expressions for the dissipative coefficients that enter these equations. These equations are obtained by requiring that the exact nonequilibrium averages are equal to local nonequilibrium averages that depend on auxiliary fields whose values satisfy evolution equations obtained using projection operator methods. The results are illustrated by deriving reaction-diffusion equations coupled to fluid hydrodynamic equations for a solution of quantum particles that can exist in two metastable states. Nonequilibrium steady states are discussed along with the reaction rate and diffusion correlation functions that characterize such states.

arXiv:2411.18713•The Journal of Chemical Physics•Feb 24, 2025•0 citations

A healthier semi-classical dynamics

Isaac Layton, Jonathan Oppenheim, Zachary Weller-Davies

We study the back-reaction of quantum systems onto classical ones. Taking the starting point that semi-classical physics should be described at all times by a point in classical phase space and a quantum state in Hilbert space, we consider an unravelling approach, describing the system in terms of a classical-quantum trajectory. We derive the general form of the dynamics under the assumptions that the classical trajectories are continuous and the evolution is autonomous, and the requirement that the dynamics is linear and completely positive in the combined classical-quantum state. This requirement is necessary in order to consistently describe probabilities, and forces the dynamics to be stochastic when the back-reaction is non-zero. The resulting equations of motion are natural generalisations of the standard semi-classical equations of motion, but since the resulting dynamics is linear in the combined classical-quantum state, it does not lead to the pathologies which usually follow from evolution laws based on expectation values. In particular, the evolution laws we present account for correlations between the classical and quantum system, which resolves issues associated with other semi-classical approaches. In addition, despite a breakdown of predictability in the classical degrees of freedom, the quantum state evolves deterministically conditioned on the classical trajectory, provided a trade-off between decoherence and diffusion is saturated. As a result, the quantum state remains pure when conditioned on the classical trajectory. To illustrate these points, we numerically simulate a number of semi-classical toy models, including one of vacuum fluctuations as a source driving the expansion of the universe. Finally, we discuss the application of these results to semi-classical gravity, and the black-hole information problem.

arXiv:2208.11722•Quantum•Dec 16, 2024•20 citations
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.
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.
A new approximate solution to the quantum-classical Liouville equation is derived starting from the formal solution of this equation in forward-backward form. The time evolution of a mixed quantum-classical system described by this equation is obtained in a coherent state basis using the mapping representation, which expresses N quantum degrees of freedom in a 2N-dimensional phase space. The solution yields a simple dynamics in which a set of N coherent state coordinates evolves in forward and backward trajectories, while the bath coordinates evolve under the influence of the mean potential that depends on these forward and backward trajectories. It is shown that the solution satisfies the differential form of the quantum-classical Liouville equation exactly. Relations to other mixed quantum-classical and semi-classical schemes are discussed.
A new approximate solution to the quantum-classical Liouville equation is derived starting from the formal solution of this equation in forward-backward form. The time evolution of a mixed quantum-classical system described by this equation is obtained in a coherent state basis using the mapping representation, which expresses N quantum degrees of freedom in a 2N-dimensional phase space. The solution yields a simple non-Hamiltonian dynamics in which a set of N coherent state coordinates evolve in forward and backward trajectories while the bath coordinates evolve under the influence of the mean potential that depends on these forward and backward tra-jectories. It is shown that the solution satisfies the differential form of the quantum-classical Liouville equation exactly. Relations to other mixed quantum-classical and semi-classical schemes are discussed. 1 ar
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 dynamics of a quantum system coupled to a classical environment and subject to constraints that drive it out of equilibrium are described. The evolution of the system is governed by the quantum-classical Liouville equation. Rather than evaluating the evolution of the mixed quantum-classical density operator, we derive exact equations of motion for the nonequilibrium average values of a set of operators or variables, along with correlation function expressions for the dissipative coefficients that enter these equations. These equations are obtained by requiring that the exact nonequilibrium averages are equal to local nonequilibrium averages that depend on auxiliary fields whose values satisfy evolution equations obtained using projection operator methods. The results are illustrated by deriving reaction-diffusion equations coupled to fluid hydrodynamic equations for a solution of quantum particles that can exist in two metastable states. Nonequilibrium steady states are discussed along with the reaction rate and diffusion correlation functions that characterize such states.
We study the back-reaction of quantum systems onto classical ones. Taking the starting point that semi-classical physics should be described at all times by a point in classical phase space and a quantum state in Hilbert space, we consider an unravelling approach, describing the system in terms of a classical-quantum trajectory. We derive the general form of the dynamics under the assumptions that the classical trajectories are continuous and the evolution is autonomous, and the requirement that the dynamics is linear and completely positive in the combined classical-quantum state. This requirement is necessary in order to consistently describe probabilities, and forces the dynamics to be stochastic when the back-reaction is non-zero. The resulting equations of motion are natural generalisations of the standard semi-classical equations of motion, but since the resulting dynamics is linear in the combined classical-quantum state, it does not lead to the pathologies which usually follow from evolution laws based on expectation values. In particular, the evolution laws we present account for correlations between the classical and quantum system, which resolves issues associated with other semi-classical approaches. In addition, despite a breakdown of predictability in the classical degrees of freedom, the quantum state evolves deterministically conditioned on the classical trajectory, provided a trade-off between decoherence and diffusion is saturated. As a result, the quantum state remains pure when conditioned on the classical trajectory. To illustrate these points, we numerically simulate a number of semi-classical toy models, including one of vacuum fluctuations as a source driving the expansion of the universe. Finally, we discuss the application of these results to semi-classical gravity, and the black-hole information problem.