Complex environments, such as molecular matrices and biological material, play a fundamental role in many important dynamic processes in condensed phases. Because it is extremely difficult to conduct full quantum dynamics simulations on such environments due to their many degrees of freedom, here, we treat in detail the environment only around the main system of interest (the subenvironment), while the other degrees of freedom needed to maintain the equilibrium temperature are described by a simple harmonic bath, which we call a quantum thermostat. The noise generated by the subenvironment is spatially non-local and non-Gaussian and cannot be characterized by the fluctuation-dissipation theorem. We describe this model by simulating the dynamics of a two-level system (TLS) that interacts with a subenvironment consisting of a one-dimensional XXZ spin chain. The hierarchical Schrödinger equations of motion are employed to describe the quantum thermostat, allowing for time-irreversible simulations of the dynamics at arbitrary temperature. To see the effects of a quantum phase transition of the subenvironment, we investigate the decoherence and relaxation processes of the TLS at zero and finite temperatures for various values of the spin anisotropy. We observed the decoherence of the TLS at finite temperature even when the anisotropy of the XXZ model is enormous. We also found that the population-relaxation dynamics of the TLS changed in a complex manner with the change in the anisotropy and the ferromagnetic or antiferromagnetic orders of spins.
Complex environments, such as molecular matrices and biological material, play a fundamental role in many important dynamic processes in condensed phases. Because it is extremely difficult to conduct full quantum dynamics simulations on such environments due to their many degrees of freedom, here, we treat in detail the environment only around the main system of interest (the subenvironment), while the other degrees of freedom needed to maintain the equilibrium temperature are described by a simple harmonic bath, which we call a quantum thermostat. The noise generated by the subenvironment is spatially non-local and non-Gaussian and cannot be characterized by the fluctuation-dissipation theorem. We describe this model by simulating the dynamics of a two-level system (TLS) that interacts with a subenvironment consisting of a one-dimensional XXZ spin chain. The hierarchical Schrödinger equations of motion are employed to describe the quantum thermostat, allowing for time-irreversible simulations of the dynamics at arbitrary temperature. To see the effects of a quantum phase transition of the subenvironment, we investigate the decoherence and relaxation processes of the TLS at zero and finite temperatures for various values of the spin anisotropy. We observed the decoherence of the TLS at finite temperature even when the anisotropy of the XXZ model is enormous. We also found that the population-relaxation dynamics of the TLS changed in a complex manner with the change in the anisotropy and the ferromagnetic or antiferromagnetic orders of spins.
When modelling real physical systems, one should always consider the influence of the environment.The study of open quantum systems aims precisely at the assessment of the role of the environment in the evolution of these systems.The computational study of these systems is in general a formidable task.However, in many instances one can approximate the environment by means of classical dynamics.For such cases, when the environment is in a thermal state, one can generalize non-Hamiltonian equations of motion to quantum-classical dynamics in order to control the thermodynamic temperature.Such a generalization is achieved within the partial Wigner representation of quantum mechanics, and leads to a quantum-classical evolution of the system.We have studied how to simulate a thermal bath, by means of the least possible number of additional extended variables [1].This has been achieved upon reformulating the Nosé-Hoover Power (NHP) thermostat in quantum-classical theory.When applied to the dissipative evolution of a quantum spin the NHP thermostat, obtained numerical results in agreement with those obtained using Nosè-Hoover chains.However, since a fewer number of variables are needed to achieve the correct sampling of the canonical distribution at equilibrium, the NHP thermostat promises to be better suited for the simulation of low dimensional open quantum system on discrete grids.Following this the quantum dynamics of a XXZ two-spin chain interacting with multiple bosonic baths have been studied, within the mixed Wigner-Heisenberg representation of quantum mechanics.In particular, we have simulated the dynamics of the reduced density matrix ρ s (t), along with the time evolution of the quantum entanglement of the two spins via the concurrence C (ρ).
An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.
An open quantum system is a system coupled to an environment that can describe time-irreversible dynamics through which the system evolves toward the thermal equilibrium state.We present a quantum mechanically rigorous theory in order to help the analysis of spectra obtained from the advanced nuclear magnetic resonance (NMR) and muon spin rotation, relaxation or resonance (μSR) techniques.Our approach is based on the numerically "Exact" hierarchical equations of motion (HEOM) approach, which allows us to study the reduced system dynamics for non-perturbative and non-Markovian system-bath interactions at finite temperature even under strong time-dependent perturbations.We demonstrate the present theory to analyze μSR and low-field NMR spectra, as an extension of the Kubo-Toyabe theory focusing on the effects of temperature and anisotropy of a local magnetic field on spectra, to help further the development of these experimental means.
We explore some aspects of nonequilibrium statistical mechanics of classical and quantum systems. Two chapters are devoted to fluctuation theorems which were originally derived for classical systems. The main challenge in formulating them in quantum mechanics is the fact that fundamental quantities of interest, like work, are defined via the classical concept of a phase space trajectory. We utilize the decoherent histories conceptual framework, in which classical trajectories emerge in quantum mechanics as a result of coarse graining, and provide a first-principles analysis of the nonequilibrium work relation of Jarzynski and Crooks's fluctuation theorem for a quantum system interacting with a general environment based on the quantum Brownian motion (QBM) model. We indicate a parameter range at low temperatures where the theorems might fail in their original form. Fluctuation theorems of Jarzynski and Crooks for systems obeying classical Hamiltonian dynamics are derived under the assumption that the initial conditions are sampled from a canonical ensemble, even though the equilibrium state of an isolated system is typically associated with the microcanonical ensemble. We address this issue through an exact analysis of the classical Brownian motion model. We argue that a stronger form of ensemble equivalence than usually discussed in equilibrium statistical mechanics is required for these theorems to hold in the infinite environment limit irrespective of the ensemble used, and proceed to prove it for this model. An exact expression for the probability distribution of work is obtained for finite environments. Intuitively one expects a system to relax to an equilibrium state when brought into contact with a thermal environment. Yet it is important to have rigorous results that provide conditions for equilibration and characterize the equilibrium state. We consider the dynamics of open quantum systems using the Langevin and master equations and rigorously show that under fairly general conditions quantum systems interacting with a heat bath relax to the equilibrium state defined as the reduced thermal state of the system plus environment, even in the strong coupling regime. Our proof is valid to second-order in interaction strength for general systems and exact for the linear QBM model, for which we also show the equivalence of multi-time correlations. In the final chapter we give a sampling of our investigations into macroscopic quantum phenomena. We work out in detail a specific example of how and under what conditions the center of mass (CoM) coordinate of a macroscopic object emerges as the relevant degree of freedom. Interaction patterns are studied in terms of the couplings they induce between the CoM and relative coordinates of two macroscopic objects. We discuss the implications of these interaction patterns on macroscopic entanglement.
A novel method is introduced in order to treat the dissipative dynamics of quantum systems interacting with a bath of classical degrees of freedom. The method is based upon an extension of the Nosè–Hoover chain (constant temperature) dynamics to quantum-classical systems. Both adiabatic and nonadiabatic numerical calculations on the relaxation dynamics of the spin-boson model show that the quantum-classical Nosè–Hoover chain dynamics represents the thermal noise of the bath in an accurate and simple way. Numerical comparisons, both with the constant-energy calculation and with the quantum-classical Brownian motion treatment of the bath, show that the quantum-classical Nosè–Hoover chain dynamics can be used to introduce dissipation in the evolution of a quantum subsystem even with just one degree of freedom for the bath. The algorithm can be computationally advantageous in modelling, within computer simulation, the dynamics of a quantum subsystem interacting with complex molecular environments.
The work presented in this thesis is divided into two main parts: strong zero modes in quantum spin chains with open boundaries and large deviations in classical spin systems. The results concerning the strong zero modes are presented in chapters 2 and 3, while chapter 4 contains the results about large deviations in systems of many, but independent, classical degrees of freedom. These two distinct classes of problems in spin chains are both about aspects of slow relaxation to equilibrium in many-body systems. Nowadays slow dynamics and non-ergodicity in many-body quantum systems are receiving much attention due to their notable technological potential for storing and processing, for very long times, the quantum information encoded in boundary degrees of freedom. In the first part of the work presented in this thesis we mainly focus on quantum spin chains with open boundaries, which can display coherence times for edge spins, that diverge with the system size, as a consequence of the emergence, in the thermodynamic limit, of conserved local quantities, called strong zero modes. First, we discuss the fate of these coherence times when dynamics becomes dissipative. Understanding the behavior of these long correlation times in open quantum systems, undergoing Markovian (memoryless) dynamics, would allow for the possibility of exploiting the protected information encoded in the boundary spins for applications in quantum devices in realistic non-equilibrium settings. Second, we study one possible simple generalization of the transverse field Ising model to a class of ladder spin systems with plaquette interactions in the presence of a transverse field. Also for these systems one can explicitly construct the strong zero mode operators and understand the behavior of the consequential long coherence times for the edge degrees of freedom under isolated (unitary) dynamics. This is possible because the models we consider, while appearing complex superficially, can be brought to simple forms by exploiting an extensive number of conserved quantities they possess. The results presented in chapters 2 and 3 are closely related, even if they focus on different models and investigate different aspects of their dynamics. Nevertheless, the underlying concept they are based on is the same one, that is, an eigenstate phase transition for the strong zero modes that has interesting consequences for the full energy spectrum. In the second part of the thesis we focus on a different type of phase transitions. Here we address the issue of dynamical large deviation phase transitions from a disordered phase to an ordered one, in classical systems which consist in a collection of independent Ising spins evolving stochastically. Surprisingly, by reweighing trajectories accordingly to a certain type of time-integrated observable, even a system of independently evolving spins can display some critical behavior, at the fluctuation level, in trajectory space in the large size and long time limit. By making use of the tilted generator techniques we study the statistics of two particular time-integrated observables. Their associated tilted generators have an exact Kramers-Wannier duality and they look like, up to a sign and an additive constant, some specific well known quantum Hamiltonians. In this sense, the dynamical large deviation phase transitions we observe in these classical stochastic systems translate in ground state quantum phase transitions taking place at the self dual point of the tilted generators we study.