From the motion of hurricanes to the cooling of a cup of tea, dissipative dynamics are ubiquitous in nature. Unfortunately, when studying lattice-based strongly correlated quantum systems, a common focus of experiments and highly relevant for material science, the curse of dimensionality prevents us from directly solving the exact equations describing the dissipative motion. We instead often use approximations and numerical analysis, as few analytical techniques currently exist which can provide useful insights. This thesis studies the constraints that non-abelian algebraic structures place on the dynamics of such non-equilibrium many-body quantum systems. In doing so, we develop novel algebraic methods that provide previously absent analytic insight into a range of problems. We apply these new techniques to study dynamical phases of quantum matter that are of both theoretical and experimental interest. We study how simple algebraic structures imply persistent non-stationarity in dissipative many-body quantum systems. These results enable us to develop a theory of quantum synchronisation based on symmetries and to explore a novel class of time crystals where order is induced by noise. These dynamical phases of matter have been gaining recent interest for their insights into thermalisation and their potential application in quantum technologies. We therefore use our results to explore experimentally realisable models. We also use the general principles of prethermalisation to propose a new theory to potentially explain recent experiments studying light-induced superconductivity in terms of eta-pairing and approximate symmetries. We discuss the merits and weaknesses of our proposal in relation to the experimental evidence. We also indicate techniques that, within our theory, stabilise superconducting states by harnessing the quantum Zeno effect For certain systems, known as integrable, there are enough symmetries to allow for exact solutions using the Bethe ansatz. We demonstrate that Bethe ansatz techniques can be extended and applied to non-equilibrium problems. As an example, we use this method to analytically study the dynamical phase transitions of an XXZ spin chain with boundary loss.
From the motion of hurricanes to the cooling of a cup of tea, dissipative dynamics are ubiquitous in nature. Unfortunately, when studying lattice-based strongly correlated quantum systems, a common focus of experiments and highly relevant for material science, the curse of dimensionality prevents us from directly solving the exact equations describing the dissipative motion. We instead often use approximations and numerical analysis, as few analytical techniques currently exist which can provide useful insights. This thesis studies the constraints that non-abelian algebraic structures place on the dynamics of such non-equilibrium many-body quantum systems. In doing so, we develop novel algebraic methods that provide previously absent analytic insight into a range of problems. We apply these new techniques to study dynamical phases of quantum matter that are of both theoretical and experimental interest. We study how simple algebraic structures imply persistent non-stationarity in dissipative many-body quantum systems. These results enable us to develop a theory of quantum synchronisation based on symmetries and to explore a novel class of time crystals where order is induced by noise. These dynamical phases of matter have been gaining recent interest for their insights into thermalisation and their potential application in quantum technologies. We therefore use our results to explore experimentally realisable models. We also use the general principles of prethermalisation to propose a new theory to potentially explain recent experiments studying light-induced superconductivity in terms of eta-pairing and approximate symmetries. We discuss the merits and weaknesses of our proposal in relation to the experimental evidence. We also indicate techniques that, within our theory, stabilise superconducting states by harnessing the quantum Zeno effect For certain systems, known as integrable, there are enough symmetries to allow for exact solutions using the Bethe ansatz. We demonstrate that Bethe ansatz techniques can be extended and applied to non-equilibrium problems. As an example, we use this method to analytically study the dynamical phase transitions of an XXZ spin chain with boundary loss.
Remarkable experimental progress and theoretical advancements in the study of strongly corre- lated systems have opened up new doors towards a better understanding of interesting states of matter. Such states could possibly be used to develop novel technological applications based on the principles of quantum mechanics. In recent years, the investigation of both complex solid state and ultracold atom systems has lead to a number of fascinating discoveries with regards to their non-equilibrium properties. One route to probe non-equilibrium physics is to consider the evolution of open systems defined by the coupling to an external environment introducing dissipation. Engineering the environment according to predefined requirements can then be used to dynamically prepare states which are difficult to access otherwise. This thesis explores the dynamics of many-body open quantum systems aiming at under- standing the effects resulting from the interplay of kinetic motion, interactions and dissipation in different physical settings. First, we build upon recent developments concerning the coupling of ultracold atoms to optical resonators which are subject to photon losses. We show how the dissipation can trigger an attractor dynamics towards a topologically non-trivial steady state that is robust against perturbations. The numerically obtained full open system evolution en- ables us to not only characterize in detail the steady state but also to analyze the features of the approaching dynamics. Second, we study the propagation of two-time correlation functions in interacting quantum spin-1/2 chains in contact with an environment causing dephasing. Pro- vided with the numerically exact full time-evolution, we identify an algebraic scaling regime where time-translation invariance breaks down and uncover the emergence of aging dynam- ics. Finally, we study a system of interacting fermionic atoms in a one-dimensional lattice with particle loss at the central site. There we report on the dissipation-induced deceleration of the non-equilibrium dynamics known as the quantum Zeno effect. As in the majority of cases quantum many-body systems are not solvable analytically, another central aspect of this thesis deals with adapting numerical methods based on matrix product state algorithms to investigate open quantum systems dynamics.
Abstract We develop a Bethe ansatz based approach to study dissipative systems experiencing loss. The method allows us to exactly calculate the spectra of interacting, many-body Liouvillians. We discuss how the dissipative Bethe ansatz opens the possibility of analytically calculating the dynamics of a wide range of experimentally relevant models including cold atoms subjected to one and two body losses, coupled cavity arrays with bosons escaping the cavity, and cavity quantum electrodynamics. As an example of our approach we study the relaxation properties in a boundary driven XXZ spin chain. We exactly calculate the Liouvillian gap and find different relaxation rates with a novel type of dynamical dissipative phase transition. This physically translates into the formation of a stable domain wall in the easy-axis regime despite the presence of loss. Such analytic results have previously been inaccessible for systems of this type.
In this thesis we study problems relating the the structure and simulation of entangled many body quantum systems, their utility in adiabatic quantum computation, and the influence of the environment in thermalizing the system and degrading the usefulness of quantum dynamics in this setting. We then study a particular strongly coupled many body quantum system in order to better understand when quantum systems do not thermalize in this manner. In chapter 2 of this thesis we study the properties of quantum dynamics restricted to an efficiently representable sub-manifold of quantum states both the finite and infinite chain of spin- 1=2 subsystems. We investigate the trade-off between gains in efficiency due to this restriction against losses in fidelity. We find the integration to be very stable and shows significant gains in efficiency compared to the naively related matrix product states. However, much of this advantage is offset by a significant reduction in fidelity. We investigate the effect of explicit symmetry breaking in the ansatz and formulate the principles for determining when correlator product states may be a useful tool. We find that scaling with overlap/bond order may be more stable with correlator product states allowing a more efficient extraction of critical exponents and present an example in which the use of correlator product states is orders of magnitude quicker than matrix product states. In chapters 3, 4 and 5 we extend this picture to allow for the study of the dissipative and decohering dynamics of a quantum system interacting with a bath, and pay particular reference to its effect on adiabatic quantum computation. In chapter 3 we consider a system of mutually interacting superconducting flux qubits coupled to a thermal bath that generalises the dissipative model of Landau-Lifschitz-Gilbert to the case of anisotropic bath couplings. We show that the dissipation acts to bias the quantum trajectories towards a reduced phase space. We study the model in the context of the D-Wave computing device and recover dynamics closely related to several models proposed on phenomenological grounds. In chapter 4 we extend this analysis to study explicitly the influence of dissipative dynamics on the lifetime of entanglement. In chapter 5 we apply this understanding to develop a methodology for benchmarking the quantum correlations harnessed by an adiabatic computation and apply this process to the D-Wave Vesuvius machine. Further developing this interest in the effect of thermalisation of quantum dynamics in chapter 6 we consider systems which fail to thermalise even in the presence of strong coupling to their surroundings. This many body localised behaviour has been recently established to be a robust phase of matter in the presence of strong disorder in one dimension. Here we show the the low lying energy states of a many body system contain immobile excitations, this immobility results in an transition in the character of low lying eigenstates at arbitrarily weak disorder. This represents a novel appearance of localising behaviour in many body systems. Finally we consider possible avenues for future work stemming from this thesis.
The study of quantum dynamics is an increasingly important aspect of atomic, molecular and optical physics, as well as condensed-matter physics and other disciplines. It is essential in order to understand non-equilibrium phenomena which do not fit within the framework of equilibrium thermodynamics, such as many-body localization, dynamical quantum phase transitions, and non-equilibrium steady states in driven-dissipative systems. However, the exact calculation of dynamical observables is extremely difficult in general for several reasons, which makes it challenging for theoretical predictions and analyses to keep pace with rapidly advancing experiments. In this thesis I use and develop several computational methods to investigate the dynamics of quantum many-body systems, with the two goals of better understanding the evolution of local observables and correlation functions in spin and Fermion systems as well as creating and testing more effective numerical methods for solving quantum dynamics problems. I explore quantum dynamics by using Hamiltonian parameter quenches in several commonly-used models in AMO physics and elsewhere: the Fermi-Hubbard model, the XXZ model, and the transverse-Ising model. The first represents a system of short-range interacting Fermions, while the latter two represent spin systems with different spin-spin couplings and magnetic field terms. All of these models can be realized in current experiments. First, I derive an exact solution for an extreme interaction quench in the Fermi-Hubbard model, from a strongly-interacting initial state to a noninteracting final Hamiltonian. From this result I show that there are nontrivial transient connected correlations, despite the fact that the system has very high temperature and a noninteracting Hamiltonian. Those features would suppress such correlations in equilibrium, but the dynamical correlations are still significant. Then, I develop and examine two numerical methods for solving dynamics of the spin expectation values and correlation functions in the transverse-Ising and XXZ models. Both methods rely on approximating a thermodynamically large system with a much smaller and more computationally tractable cluster of lattice sites, or a series of clusters. The first method is a dynamical numerical linked cluster expansion (NLCE). NLCEs have already been used successfully on systems in thermodynamic equilibrium, and I extended the technique to directly compute the dynamics of an observable after a quantum quench. I review the general NLCE method, describe my dynamical variant of it, and compute the quench dynamics of spin and correlation observables in the transverse-Ising and XXZ models. The second method, which I term the adaptive-boundary cluster method, relies on recently refined bounds on information propagation in quantum systems to choose a cluster of lattice sites that minimizes the finite size error bound of a dynamical computation. I briefly review these bounds, and describe a new algorithm which uses them to construct a cluster whose boundary is chosen according to a related error estimate. Such an adaptive-boundary cluster has the property that there exists a computable upper bound on the finite size error of observables measured at the cluster center, and this error bound is nearly minimal when compared to that of other clusters of the same size. I then compute the quench dynamics of spin and correlation observables in the transverse-Ising model using the adaptive boundary clusters. I apply this method to both uniform and strongly disordered systems. To evaluate the effectiveness of both dynamical NLCE and adaptive-boundary clusters, I compare the results obtained using them to results from a standard and general purpose numerical method: exact solvers applied a rectangular cluster with periodic boundary conditions. Both of the new methods usually converge faster and are more accurate than the standard method for comparable cluster sizes. This thesis illustrates that choosing numerical methods which take advantage of the geometric and entanglement structure of the initial state and Hamiltonian can substantially enhance our ability to study quantum dynamics, by reducing the computational cost of numerically exact solutions. I conclude with a discussion of various avenues for future research to build on this work.
We use the Bethe Ansatz technique to study dissipative systems experiencing loss. The method allows us to exactly calculate the Liouvillian spectrum. This opens the possibility of analytically calculating the dynamics of a wide range of experimentally relevant models including cold atoms subjected to one and two body losses, coupled cavity arrays with bosons escaping the cavity, and cavity quantum electrodynamics. As an example of our approach we study the relaxation properties in a boundary driven XXZ spin chain. We exactly calculate the Liouvillian gap and find different relaxation rates with a novel type of dynamical dissipative phase transition. This physically translates into the formation of a stable domain wall in the easy-axis regime despite the presence of loss. Such analytic results have previously been inaccessible for systems of this type.
This thesis studies non-equilibrium open quantum systems where the dissipation is crucial to the achievement of novel physical regimes. We focus on atomic systems which allow for the coupling of a ground state to a Rydberg state, relying on the strong interactions between Rydberg atoms to produce the collective behaviour that we aim to investigate. For atoms in an optical lattice undergoing standard dissipation forms, e.g. loss and dephasing, we find these simple settings allow for the production of models contained in the non-equilibrium realm. \nWe start by looking at a system with engineered pair dissipation on a one-dimensional lattice. When the dissipation is strong relative to a tunnelling process it creates a quantum Zeno effect which projects the system onto a Zeno-subspace. This subspace is found to contain complexes which experience a binding due to the dissipation. The properties of these complexes are found to feature spin-orbit coupling and, in certain instances, a flat band. \nWe then study what kinetically constrained models (KCMs) can be reproduced in a lattice system. KCMs are models which typically feature trivial steady states, but a complex relaxation dynamics. These models appear in the fields of glasses and soft matter physics. We find a general framework for the consideration of a quantum Hamiltonian and a classical potential with strong dephasing noise. We then focus on a model mimicking volume excluded KCMs and find characteristic constrained behaviour, such as ergodicity breaking. \nWe apply this framework to the decay of a many-body localised state in an open system with interactions in which we find the decay to be classical in the two interaction limits. For weak interactions, it follows a stretched exponential form due to pair relaxation, while for strong interactions the decay follows a compressed exponential, now being modelled as an Avrami process due to the correlated relaxation. We also find that on-site loss only affects the strong interacting limit. \nWe then move on to the study of universal non-equilibrium behaviour in the directed percolation (DP) class. We consider on-site atomic loss and gain as a substitute for the standard decay channel. We show that this replaces the absorbing state with an enlarged absorbing space, leading to a loss of the DP transition at lower average densities. This class of DP-like systems has received little study, and we present a method of experimentally realising it in current set-ups. \nWe finish with a look at a quantum DP model, where we consider its quantum and classical limits. We find that the transition changes from first to second order as the system becomes more classical, featuring a bi-critical point. We then numerically demonstrate that the same transitions are visible in idealised and Rydberg models.