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

The Journal of Physical Chemistry B

Trotter-Based Simulation of Quantum-Classical Dynamics

Dónal MacKernan, Giovanni Ciccotti, Raymond Kapral

Published Dec 22, 2007• 105 citations
Publisher

Abstract

Quantum rate processes in condensed phase systems are often computed by combining quantum and classical descriptions of the dynamics. An algorithm for simulating the quantum-classical Liouville equation, which describes the dynamics of a quantum subsystem coupled to a classical bath, is presented in this paper. The algorithm is based on a Trotter decomposition of the quantum-classical propagator, in conjunction with Monte Carlo sampling of quantum transitions, to yield a surface-hopping representation of the dynamics. An expression for the nonadiabatic propagator that is responsible for quantum transitions and associated bath momentum changes is derived in a form that is convenient for Monte Carlo sampling and exactly conserves the total energy of the system in individual trajectories. The expectation values of operators or quantum correlation functions can be evaluated by initial sampling of quantum states and use of quantum-classical Liouville dynamics for the time evolution. The algorithm is tested by calculations on the spin-boson model, for which exact quantum results are available, and is shown to reproduce the exact results for stronger nonadiabatic coupling and much longer times using fewer trajectories than other schemes for simulating quantum-classical Liouville dynamics.

DOI
10.1021/jp0761416

Related papers

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

Surface-hopping dynamics of a spin-boson system

Dónal MacKernan, Giovanni Ciccotti, Raymond Kapral

The spin-boson model is solved within the framework of quantum-classical dynamics using our recently-developed surface-hopping scheme. The quantum-classical equation of motion is expressed in an adiabatic basis and its solution is constructed from an ensemble of trajectories which undergo nonadiabatic transitions and evolve coherently on the adiabatic surfaces. Details of the algorithm for the simulation of the dynamics are presented and the method of simple Monte Carlo sampling used to evaluate the expectation values of observables is discussed. The simulation method is applied to a spin-boson system with a harmonic bath composed of ten oscillators with an Ohmic spectral density. For the spin-boson model the present implementation of quantum-classical dynamics is exact and the results of our surface-hopping simulations are in accord with previous numerically exact results for this model.

The Journal of Chemical Physics•Feb 8, 2002•101 citations

Computer simulation of nonadiabatic dynamics by means of the quantum-classical Liouville equation.

Daniel A. Uken

Simulation of quantum dynamics for many-body systems is an open area of research. For interacting many-body quantum systems, the computer memory necessary to perform calculations has an astronomical value, so that approximated models are needed to reduce the required computational resources. A useful approximation that can often be made is that of quantum-classical dynamics, where the majority of the degrees are treated classically, while a few of them must be treated quantum mechanically. When energy is exchanged very quickly between the quantum subsystem and classical environment, the dynamics is nonadiabatic. Most theories for nonadiabatic dynamics are unsatisfactory, as they fail to properly describe the quantum backreaction of the subsystem on the environment. However, an approach based on the quantum-classical Liouville equation solves this problem. Even so, nonadiabatic dynamics is difficult to implement on a computer, and longer simulation times are often inaccessible due to statistical error. There is thus a need for improved algorithms for nonadiabatic dynamics. In this thesis, two algorithms that utilise the quantum-classical Liouville equation will be qualitatively and quantitatively compared. In addition, stochastic sampling schemes for nonadiabatic transitions will be studied, and a new sampling scheme is introduced [D. A. Uken et al., Phys. Rev. E. 88, 033301 (2013)] which proves to have a dramatic advantage over existing techniques, allowing far longer simulation times to be calculated reliably.

•Jan 1, 2013•0 citations

Mixed quantum-classical surface hopping dynamics

Steven O. Nielsen, Raymond Kapral, Giovanni Ciccotti

An algorithm is presented for the exact solution of the evolution of the density matrix of a mixed quantum-classical system in terms of an ensemble of surface hopping trajectories. The system comprises a quantum subsystem coupled to a classical bath whose evolution is governed by a mixed quantum-classical Liouville equation. The integral solution of the evolution equation is formulated in terms of a concatenation of classical evolution segments for the bath phase space coordinates separated by operators that change the quantum state and bath momenta. A hybrid Molecular Dynamics–Monte Carlo scheme which follows a branching tree of trajectories arising from the action of momentum derivatives is constructed to solve the integral equation. We also consider a simpler scheme where changes in the bath momenta are approximated by momentum jumps. These schemes are illustrated by considering the computation of the evolution of the density matrix for a two-level system coupled to a low dimensional classical bath.

The Journal of Chemical Physics•Apr 15, 2000•156 citations

Quantum-Classical Master Equation Dynamics: An Analysis of Decoherence and Surface-hopping Techniques

Robbie Grunwald

In this thesis quantum-classical dynamics is applied to the study of quantum condensed phase processes. This approach is based on the quantum-classical Liouville equation where the dynamics of a small subset of the degrees of freedom are treated quantum mechanically while the remaining degrees of freedom are treated by classical mechanics to a good approximation. We use this approach as it is computationally tractable, and the resulting equation of motion accurately accounts for the quantum and classical dynamics, as well as the coupling between these two components of the system. \n\nBy recasting the quantum-classical Liouville equation into the form of a generalized master equation we investigate connections to surface-hopping. The link between these approaches is decoherence arising from interaction of the subsystem with the environment. We derive an evolution equation for the subsystem which contains terms accounting for the effects of the environment. One of these terms involves a memory kernel that accounts for the coherent dynamics. If this term decays rapidly, a Markovian approximation can be made. By lifting the resulting subsystem master equation into the full phase space, we obtain a Markovian master equation that prescribes surface-hopping-like dynamics. Our analysis outlines the conditions under which such a description is valid.\n\nNext, we consider the calculation of the rate constant for a quantum mechanical barrier crossing process. Starting from the reactive-flux autocorrelation function, we derive a quantum-classical expression for the rate kernel. This expression involves quantum-classical evolution of a species operator averaged over the initial quantum equilibrium structure of the system making it possible to compute the rate constant via computer simulation. Using a simple model for a proton transfer reaction we compare the results of the rate calculation obtained by quantum-classical Liouville dynamics with that of master equation dynamics. The master equation provides a good approximation to the full quantum-classical Liouville calculation for our model and a more stable algorithm results due to the elimination of oscillating phase factors in the simulation. \n\nFinally, we make use of the theoretical framework established in this thesis to analyze some aspects of decoherence used in popular surface-hopping techniques.

•Jan 19, 2009•0 citations

Sequential short-time propagation of quantum classical dynamics

Dónal MacKernan, Raymond Kapral, Giovanni Ciccotti

An algorithm for the simulation of quantum-classical dynamics is presented. Quantum-classical evolution is effected by a propagator exp(i (L) over capt) involving the quantum classical Liouville operator (L) over cap that describes the evolution of a quantum subsystem coupled to a classical bath. Such a mixed description provides a means to study the dynamics of complex many-body systems where certain degrees of freedom are treated quantum mechanically. The algorithm is constructed by decomposing the time interval Deltat into small segments of length At and successively applying the propagator in the short time segments to obtain the evolution for long times. The algorithm is shown to be a discretization of the iterated Dyson form of the propagator whose direct solution is vexatious. The sequential short-time propagation algorithm is applied to the spin-boson model for a range of values of the Kondo parameter and shown to be effective.

Journal of Physics Condensed Matter•Sep 26, 2002•87 citations
Quantum rate processes in condensed phase systems are often computed by combining quantum and classical descriptions of the dynamics. An algorithm for simulating the quantum-classical Liouville equation, which describes the dynamics of a quantum subsystem coupled to a classical bath, is presented in this paper. The algorithm is based on a Trotter decomposition of the quantum-classical propagator, in conjunction with Monte Carlo sampling of quantum transitions, to yield a surface-hopping representation of the dynamics. An expression for the nonadiabatic propagator that is responsible for quantum transitions and associated bath momentum changes is derived in a form that is convenient for Monte Carlo sampling and exactly conserves the total energy of the system in individual trajectories. The expectation values of operators or quantum correlation functions can be evaluated by initial sampling of quantum states and use of quantum-classical Liouville dynamics for the time evolution. The algorithm is tested by calculations on the spin-boson model, for which exact quantum results are available, and is shown to reproduce the exact results for stronger nonadiabatic coupling and much longer times using fewer trajectories than other schemes for simulating quantum-classical Liouville dynamics.
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 spin-boson model is solved within the framework of quantum-classical dynamics using our recently-developed surface-hopping scheme. The quantum-classical equation of motion is expressed in an adiabatic basis and its solution is constructed from an ensemble of trajectories which undergo nonadiabatic transitions and evolve coherently on the adiabatic surfaces. Details of the algorithm for the simulation of the dynamics are presented and the method of simple Monte Carlo sampling used to evaluate the expectation values of observables is discussed. The simulation method is applied to a spin-boson system with a harmonic bath composed of ten oscillators with an Ohmic spectral density. For the spin-boson model the present implementation of quantum-classical dynamics is exact and the results of our surface-hopping simulations are in accord with previous numerically exact results for this model.
Simulation of quantum dynamics for many-body systems is an open area of research. For interacting many-body quantum systems, the computer memory necessary to perform calculations has an astronomical value, so that approximated models are needed to reduce the required computational resources. A useful approximation that can often be made is that of quantum-classical dynamics, where the majority of the degrees are treated classically, while a few of them must be treated quantum mechanically. When energy is exchanged very quickly between the quantum subsystem and classical environment, the dynamics is nonadiabatic. Most theories for nonadiabatic dynamics are unsatisfactory, as they fail to properly describe the quantum backreaction of the subsystem on the environment. However, an approach based on the quantum-classical Liouville equation solves this problem. Even so, nonadiabatic dynamics is difficult to implement on a computer, and longer simulation times are often inaccessible due to statistical error. There is thus a need for improved algorithms for nonadiabatic dynamics. In this thesis, two algorithms that utilise the quantum-classical Liouville equation will be qualitatively and quantitatively compared. In addition, stochastic sampling schemes for nonadiabatic transitions will be studied, and a new sampling scheme is introduced [D. A. Uken et al., Phys. Rev. E. 88, 033301 (2013)] which proves to have a dramatic advantage over existing techniques, allowing far longer simulation times to be calculated reliably.
An algorithm is presented for the exact solution of the evolution of the density matrix of a mixed quantum-classical system in terms of an ensemble of surface hopping trajectories. The system comprises a quantum subsystem coupled to a classical bath whose evolution is governed by a mixed quantum-classical Liouville equation. The integral solution of the evolution equation is formulated in terms of a concatenation of classical evolution segments for the bath phase space coordinates separated by operators that change the quantum state and bath momenta. A hybrid Molecular Dynamics–Monte Carlo scheme which follows a branching tree of trajectories arising from the action of momentum derivatives is constructed to solve the integral equation. We also consider a simpler scheme where changes in the bath momenta are approximated by momentum jumps. These schemes are illustrated by considering the computation of the evolution of the density matrix for a two-level system coupled to a low dimensional classical bath.
In this thesis quantum-classical dynamics is applied to the study of quantum condensed phase processes. This approach is based on the quantum-classical Liouville equation where the dynamics of a small subset of the degrees of freedom are treated quantum mechanically while the remaining degrees of freedom are treated by classical mechanics to a good approximation. We use this approach as it is computationally tractable, and the resulting equation of motion accurately accounts for the quantum and classical dynamics, as well as the coupling between these two components of the system. \n\nBy recasting the quantum-classical Liouville equation into the form of a generalized master equation we investigate connections to surface-hopping. The link between these approaches is decoherence arising from interaction of the subsystem with the environment. We derive an evolution equation for the subsystem which contains terms accounting for the effects of the environment. One of these terms involves a memory kernel that accounts for the coherent dynamics. If this term decays rapidly, a Markovian approximation can be made. By lifting the resulting subsystem master equation into the full phase space, we obtain a Markovian master equation that prescribes surface-hopping-like dynamics. Our analysis outlines the conditions under which such a description is valid.\n\nNext, we consider the calculation of the rate constant for a quantum mechanical barrier crossing process. Starting from the reactive-flux autocorrelation function, we derive a quantum-classical expression for the rate kernel. This expression involves quantum-classical evolution of a species operator averaged over the initial quantum equilibrium structure of the system making it possible to compute the rate constant via computer simulation. Using a simple model for a proton transfer reaction we compare the results of the rate calculation obtained by quantum-classical Liouville dynamics with that of master equation dynamics. The master equation provides a good approximation to the full quantum-classical Liouville calculation for our model and a more stable algorithm results due to the elimination of oscillating phase factors in the simulation. \n\nFinally, we make use of the theoretical framework established in this thesis to analyze some aspects of decoherence used in popular surface-hopping techniques.
An algorithm for the simulation of quantum-classical dynamics is presented. Quantum-classical evolution is effected by a propagator exp(i (L) over capt) involving the quantum classical Liouville operator (L) over cap that describes the evolution of a quantum subsystem coupled to a classical bath. Such a mixed description provides a means to study the dynamics of complex many-body systems where certain degrees of freedom are treated quantum mechanically. The algorithm is constructed by decomposing the time interval Deltat into small segments of length At and successively applying the propagator in the short time segments to obtain the evolution for long times. The algorithm is shown to be a discretization of the iterated Dyson form of the propagator whose direct solution is vexatious. The sequential short-time propagation algorithm is applied to the spin-boson model for a range of values of the Kondo parameter and shown to be effective.