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

Modern Physics Letters B

MOMENTUM SHIFT IN NONADIABATIC DYNAMICS

Daniel A. Uken, Alessandro Sergi

Published Mar 31, 2011• 2 citations
Publisher

Abstract

We review the momentum-jump approximation in the quantum-classical Liouville formulation of nonadiabatic dynamics. We consider different momentum shift rules and study their numerical importance in computer simulations of the spin-boson model. Our findings might help the design of generalized sampling schemes of nonadiabatic dynamics at long time.

DOI
10.1142/s021798491102622x

Related papers

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

Numerical sampling of nonadiabatic dynamics of quantum-classical systems.

Daniel A. Uken

The simulation of the dynamics of quantum systems is very difficult, due to the fact that, in general, it cannot be calculated exactly for interacting many-body systems. Brute force simulations of quantum dynamics are simply not feasible, and approximations need to be made. In many instances a quantum system can be approximated as a quantum-classical system, where only a subsystem of interest is treated quantum mechanically, and the rest is considered as a classical bath. When energy is free to be exchanged between the subsystem and its environment, the dynamics that occur is said to be nonadiabatic. This type of dynamics is challenging to calculate on a computer, as it can lead to large statistical errors at long times. Hence, there is a need for improved algorithms for nonadiabatic dynamics. In this thesis, a recently introduced nonadiabatic sampling scheme [A. Sergi and F. Petruccione, Phys. Rev. E 81, 032101 (2010)] is used to calculate the long-time dynamics of a model system comprising a quantum spin coupled to a bath of harmonic oscillators. Also, various technical aspects of the algorithm are investigated.

•Jan 1, 2010•0 citations

Filtering schemes in the quantum-classical Liouville approach to nonadiabatic dynamics

Daniel A. Uken, Alessandro Sergi, Francesco Petruccione

We study a number of filtering schemes for the reduction of the statistical error in nonadiabatic calculations by means of the quantum-classical Liouville equation. In particular, we focus on a scheme based on setting a threshold value on the sampling weights, so that when the threshold is overcome the value of the weight is reset, and on another approach which prunes the ensemble of the allowed nonadiabatic transitions according to a generalized sampling probability. Both methods have advantages and drawbacks, however, their combination drastically improves the performance of an algorithm known as the sequential short-time step propagation [MacKernan et al., J. Phys: Condens. Matter 14, 9069 (2002)], which is derived from a simple first order expansion of the quantum-classical propagator. Such an algorithm together with the combined filtering procedures produce results that compare very well with those obtained by means of numerically accurate path integral quantum calculations for the spin-boson model, even for intermediate and strong coupling regimes.

arXiv:1306.5205•Physical Review E•Sep 4, 2013•20 citations

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

Stochastic simulation of long-time nonadiabatic dynamics

Daniel A. Uken, Alessandro Sergi, Francesco Petruccione

Using a generalized transition probability, we show how nonadiabatic calculations, within the Wigner–Heisenberg representation of quantum mechanics, can be reliably extended to far longer times than those allowed by a primitive sampling scheme. Tackling the spin-boson model as a paradigmatic example, we provide substantial numerical evidence that the statistical error can be reduced for a wide range of parameter values.

Physica Scripta•Feb 1, 2011•11 citations

Trotter-Based Simulation of Quantum-Classical Dynamics

Dónal MacKernan, Giovanni Ciccotti, Raymond Kapral

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 Journal of Physical Chemistry B•Dec 22, 2007•105 citations
We review the momentum-jump approximation in the quantum-classical Liouville formulation of nonadiabatic dynamics. We consider different momentum shift rules and study their numerical importance in computer simulations of the spin-boson model. Our findings might help the design of generalized sampling schemes of nonadiabatic dynamics at long time.
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.
The simulation of the dynamics of quantum systems is very difficult, due to the fact that, in general, it cannot be calculated exactly for interacting many-body systems. Brute force simulations of quantum dynamics are simply not feasible, and approximations need to be made. In many instances a quantum system can be approximated as a quantum-classical system, where only a subsystem of interest is treated quantum mechanically, and the rest is considered as a classical bath. When energy is free to be exchanged between the subsystem and its environment, the dynamics that occur is said to be nonadiabatic. This type of dynamics is challenging to calculate on a computer, as it can lead to large statistical errors at long times. Hence, there is a need for improved algorithms for nonadiabatic dynamics. In this thesis, a recently introduced nonadiabatic sampling scheme [A. Sergi and F. Petruccione, Phys. Rev. E 81, 032101 (2010)] is used to calculate the long-time dynamics of a model system comprising a quantum spin coupled to a bath of harmonic oscillators. Also, various technical aspects of the algorithm are investigated.
We study a number of filtering schemes for the reduction of the statistical error in nonadiabatic calculations by means of the quantum-classical Liouville equation. In particular, we focus on a scheme based on setting a threshold value on the sampling weights, so that when the threshold is overcome the value of the weight is reset, and on another approach which prunes the ensemble of the allowed nonadiabatic transitions according to a generalized sampling probability. Both methods have advantages and drawbacks, however, their combination drastically improves the performance of an algorithm known as the sequential short-time step propagation [MacKernan et al., J. Phys: Condens. Matter 14, 9069 (2002)], which is derived from a simple first order expansion of the quantum-classical propagator. Such an algorithm together with the combined filtering procedures produce results that compare very well with those obtained by means of numerically accurate path integral quantum calculations for the spin-boson model, even for intermediate and strong coupling regimes.
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.
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.