Numerical sampling of nonadiabatic dynamics of quantum-classical systems.

Abstract

The simulation of the dynamics of quantum systems is very di cult, 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.