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

Stochastic Simulation of Nonadiabatic Dynamics at Long Time

Daniel A. Uken, Alessandro Sergi, Francesco Petruccione

Published Mar 16, 2010• 3 citations
arXivHTMLPDF

Abstract

Using a generalized energy-conserving transition probability, it is shown 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, substantial numerical evidence is provided that effective integration of the dynamics can be achieved for a wide range of temperatures and friction.

DOI
10.48550/arxiv.1003.3117

Related papers

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

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

MOMENTUM SHIFT IN NONADIABATIC DYNAMICS

Daniel A. Uken, Alessandro Sergi

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.

Modern Physics Letters B•Mar 31, 2011•2 citations

Sampling of quantum dynamics at long time

Alessandro Sergi, Francesco Petruccione

The principle of energy conservation leads to a generalized choice of transition probability in a piecewise adiabatic representation of quantum(-classical) dynamics. Significant improvement (almost an order of magnitude, depending on the parameters of the calculation) over previous schemes is achieved. Perspectives for theoretical calculations in coherent many-body systems are opened.

arXiv:0912.1113•Physical Review E•Mar 3, 2010•18 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
Using a generalized energy-conserving transition probability, it is shown 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, substantial numerical evidence is provided that effective integration of the dynamics can be achieved for a wide range of temperatures and friction.
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.
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.