Methods for simulating the dynamics of composite systems, where part of the system is treated quantum mechanically and its environment is treated classically, are discussed. Such quantum-classical sys- tems arise in many physical contexts where certain degrees of freedom have an essential quantum character while the other degrees of freedom to which they are coupled may be treated classically to a good approxi- mation. The dynamics of these composite systems are governed by a quantum-classical Liouville equation for either the density matrix or the dynamical variables which are operators in the Hilbert space of the quantum subsystem and functions of the classical phase space variables of the classical environment. Solutions of the evolution equations may be formulated in terms of surface-hopping dynamics involving ensembles of tra- jectory segments interspersed with quantum transitions. The surface-hopping schemes incorporate quantum coherence and account for energy exchanges between the quantum and classical degrees of freedom. Various simulation algorithms are discussed and illustrated with calculations on simple spin-boson models but the methods described here are applicable to realistic many- body environments.
Methods for simulating the dynamics of composite systems, where part of the system is treated quantum mechanically and its environment is treated classically, are discussed. Such quantum-classical sys- tems arise in many physical contexts where certain degrees of freedom have an essential quantum character while the other degrees of freedom to which they are coupled may be treated classically to a good approxi- mation. The dynamics of these composite systems are governed by a quantum-classical Liouville equation for either the density matrix or the dynamical variables which are operators in the Hilbert space of the quantum subsystem and functions of the classical phase space variables of the classical environment. Solutions of the evolution equations may be formulated in terms of surface-hopping dynamics involving ensembles of tra- jectory segments interspersed with quantum transitions. The surface-hopping schemes incorporate quantum coherence and account for energy exchanges between the quantum and classical degrees of freedom. Various simulation algorithms are discussed and illustrated with calculations on simple spin-boson models but the methods described here are applicable to realistic many- body environments.
Solving quantum dynamics is an exponentially difficult problem. Thus, an exact numerical solution is inaccessible for any condensed matter system.\nA promising approach is to divide the system into a quantum subsystem containing degrees of freedom which are of greater interest or those which have more profound quantum character (e.g., have smaller mass) and a classical bath containing the rest of the system.\nImposing such a partition and treating the bath classically results in quantum-classical dynamics. The quantum-classical Liouville equation is a general equation in the Hilbert space of quantum degrees of freedom while it resides in the phase space of the classical degrees of freedom.\nAny numerical solution to this equation requires representation of the quantum subsystem in some basis. Solutions to this equation have been already proposed in the subsystem, adiabatic and force bases, each with its own cons and pros.\n\nIn this work, the quantum-classical equations of motion are cast in the subsystem basis and subsequently mapped to a number of fictitious harmonic oscillators.\nThe result is quantum-classical dynamics in the mapping basis which treats both quantum and classical degrees of freedom on the same footing, i.e., in phase space. Neglecting a portion of the back reaction of the quantum-subsystem to classical bath results in an expression for the time evolution of an operator (density matrix) equal to its Poisson bracket with the Hamiltonian.\nThis equation can be solved in terms of characteristics to provide a computationally tractable method for calculating quantum-classical dynamical properties. The expressions for expectation values and correlation functions in this formalism are derived.\nCalculations on spin-boson system, barrier crossing models---the so called Tully models---and the Fenna-Mathews-Olson pigments show very good agreement between the results of this method and numerical solutions to the Schrödinger equation.
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 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 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.
This thesis develops mixed quantum-classical methods to model dynamical nonequilibrium quantum phenomena in condensed phase systems, and also explores properties of 2-d systems near a quantum critical point using field theoretic techniques and purely classical simulation. Mixed quantum-classical methods involve dividing the system into a quantum sub-system coupled to a classical bath, which is useful for examining the dynamics of systems in the condensed phase since the full quantum treatment of a large number of degrees of freedom is often computationally infeasible. Although the bath is treated classically, its coupling to the quantum system precludes a Newtonian description. We implement an approach that uses the quantum-classical Liouville equation to propagate the quantum subsystem density matrix, which results in an ensemble of surface-hopping trajectories. We test this method on model systems and compare with a linearized path-integral approach. We also develop a novel approach to propagating the reduced density matrix that incorporates ideas from both methods while
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.