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.
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.
Simulating quantum spin systems at finite temperatures is an open challenge in many-body physics. This work studies the temperature-dependent spin dynamics of a pivotal compound,
FeI2, to determine if universal quantum effects can be accounted for by a phenomenological renormalization of the dynamical spin structure factor $S(\mathbf{q},\ensuremath{\omega})$ measured by inelastic neutron scattering. Renormalization schemes based on the quantum-to-classical correspondence principle are commonly applied at low temperatures to the harmonic oscillators describing normal modes. However, it is not clear how to extend this renormalization to arbitrarily high temperatures. Here we introduce a temperature-dependent normalization of the classical moments, the magnitude of which is determined by imposing the quantum sum rule, e.g., $\ensuremath{\int}d\ensuremath{\omega}d\mathbf{q}S(\mathbf{q},\ensuremath{\omega})={N}_{S}S(S+1)$ for
NS dipolar magnetic moments. We show that this simple renormalization scheme significantly improves the agreement between the calculated and measured $S(\mathbf{q},\ensuremath{\omega})$ for
FeI2 at all temperatures. Due to the coupled dynamics of dipolar and quadrupolar moments in that material, this renormalization procedure is extended to classical theories based on SU(3) coherent states, and by extension, to any
SU(N) coherent state representation of local multipolar moments.
Quantum mechanics is not logically closed with respect to the classical world.Its formalism unfolds as the quantization of a sub-set of classical Hamiltonians.The interpretation of quantum theory in terms of the measurement process inevitably requires to deal with systems composed by a mixture of both classical and quantum degrees of freedom.Moreover, when energy can flow between the quantum and classical degrees of freedom (i.e., in the case of nonadiabatic dynamics), there are more theoretical difficulties in order to obtain a fully consistent quantum-classical formalism.In order to perform calculations, one can renounce to the usual Lie algebraic structure of well-established physical theories, adopt non-Hamiltonian brackets, and obtain a formalism for the dynamics and statistics of quantum-classical systems that has an affordable computational complexity.Recent progress in the algorithms for the sampling of nonadiabatic dynamics of quantum-classical systems at long time is reviewed here.