This thesis summarizes a course of investigation of various aspects of non-equilibrium dynamics in isolated quantum systems which can be controlled to the extent that one can speak of not just realizing but rather simulating a desired physical effect. The first subject considered concerns a general question of relaxation in a large class of physical models. It is rigorously proven that equilibration can occur for arbitrary local observables despite the entire system being perfectly isolated. Various mechanisms responsible for convergence to local equilibrium are highlighted. These involve in particular the memory loss of non-Gaussian correlations following an interaction quench, a notion of Gaussian ergodicity and a proof of the emergence of translation invariance of correlations due to the presence of this symmetry in the Hamiltonian governing the evolution. These results provide the long time and large system size asymptotics facilitating a thermodynamic limit, but at the same time are relevant for state-of-the-art quantum simulation experiments with large numbers of ultra-cold atoms: A related effect has been observed in a one-dimensional phononic quantum field simulator and additionally a method is provided to study relaxation dynamics of this type in optical lattice quantum simulators. Within the second theme explored in this thesis a novel quantum read-out method is proposed and applied in continuous field quantum simulators which allowed for the first time to measure experimentally various thermodynamical properties of one-dimensional quasi-condensates. In particular, tomographic results concerning thermal properties, non-commuting observables, momentum and time-resolved occupation numbers of phonons are presented. Finally, ideas for practical benchmarking of the dynamics of certain closed quantum systems are put forward, based on the concept of a fidelity witness. It is demonstrated that fidelity, despite being a sensitive measure for large systems, can be efficiently estimated for non-equilibrium dynamics in coherent quantum simulators implementing paradigmatic models of condensed matter physics. The method developed has already found an independent application in studies of variational quantum circuits aiming at achieving so-called quantum chemistry accuracy using the Sycamore quantum processor. The fact that all three themes of research laid out in this thesis have found an experimental realization hints at a prognosis for future developments in physics that it will become standard that quantum simulators will realize experimentally novel theoretical ideas on demand and the time between theoretical insights and experimental observations will be dramatically shortened.
This thesis summarizes a course of investigation of various aspects of non-equilibrium dynamics in isolated quantum systems which can be controlled to the extent that one can speak of not just realizing but rather simulating a desired physical effect. The first subject considered concerns a general question of relaxation in a large class of physical models. It is rigorously proven that equilibration can occur for arbitrary local observables despite the entire system being perfectly isolated. Various mechanisms responsible for convergence to local equilibrium are highlighted. These involve in particular the memory loss of non-Gaussian correlations following an interaction quench, a notion of Gaussian ergodicity and a proof of the emergence of translation invariance of correlations due to the presence of this symmetry in the Hamiltonian governing the evolution. These results provide the long time and large system size asymptotics facilitating a thermodynamic limit, but at the same time are relevant for state-of-the-art quantum simulation experiments with large numbers of ultra-cold atoms: A related effect has been observed in a one-dimensional phononic quantum field simulator and additionally a method is provided to study relaxation dynamics of this type in optical lattice quantum simulators. Within the second theme explored in this thesis a novel quantum read-out method is proposed and applied in continuous field quantum simulators which allowed for the first time to measure experimentally various thermodynamical properties of one-dimensional quasi-condensates. In particular, tomographic results concerning thermal properties, non-commuting observables, momentum and time-resolved occupation numbers of phonons are presented. Finally, ideas for practical benchmarking of the dynamics of certain closed quantum systems are put forward, based on the concept of a fidelity witness. It is demonstrated that fidelity, despite being a sensitive measure for large systems, can be efficiently estimated for non-equilibrium dynamics in coherent quantum simulators implementing paradigmatic models of condensed matter physics. The method developed has already found an independent application in studies of variational quantum circuits aiming at achieving so-called quantum chemistry accuracy using the Sycamore quantum processor. The fact that all three themes of research laid out in this thesis have found an experimental realization hints at a prognosis for future developments in physics that it will become standard that quantum simulators will realize experimentally novel theoretical ideas on demand and the time between theoretical insights and experimental observations will be dramatically shortened.
In this Thesis I present the results of four research projects I carried out during my Ph.D. training at SISSA. The broad topic of my research was the mechanism of relaxation in isolated, integrable quantum systems. Roughly speaking, integrable quantum systems are one-dimensional lattice models or field theories that feature an extensively large number of conservation laws beyond energy and momentum conservation, which greatly influence their dynamical behaviour. Integrable systems have enjoyed a considerable interest in the last decade, mainly due to the development of ultracold atomic techniques. These techniques made it possible to implement experimentally a wealth of one-dimensional quantum models whose significance had been purely theoretical before. The theoretical interest in integrable systems stems in part from the fact that in many aspects they are exactly solvable. Exact analytical formulas are available for their spectrum and their equilibrium thermodynamic properties, creating an incentive to study analytically their non-equilibrium behaviour as well. Exact results for the dynamics of integrable systems may serve as a benchmark for confirming theories in statistical mechanics or testing approximative methods. To this day, computing analytically the entire time evolution of observables in a genuinely interacting integrable system is prohibitively difficult in general. In this Thesis, I concentrate on the description of the system at long times after it is brought out of equilibrium. Homogeneous, translationally invariant integrable systems are expected to relax to a stationary state, in which local observables are described by a generalization of the Gibbs ensemble that takes into consideration all the conservation laws beyond the energy. Two of the projects presented here concern this relaxed state. In Chapter 2, a method is presented for computing the Rényi entropies of the Generalized Gibbs Ensemble describing the relaxed state in the integrable Heisen- berg XXZ spin chain. Rényi entropies provide information about the relaxed state and its entanglement properties, going beyond the ordinary von Neumann thermo- dynamic entropy. The method presented there is very general; in principle it works in any integrable model and any translationally invariant thermodynamic macrostate. An alternative to the GGE in describing relaxed states after a quantum quench is the Quench Action approach. This approach is based on the explicit analytical knowledge of the overlaps between the initial state and the Hamiltonian eigenstates. In Chapter 3, this approach is generalized to a spin-1 integrable spin chain, the SU(3) invariant Lai–Sutherland model. The evolution of entanglement is then studied using the quasiparticle picture, which features ballistically propagating classical particle pairs carrying entanglement. If the requirement of translational invariance of the initial state is relaxed, various transport phenomena arise in integrable systems. The evolution of slowly varying, smoothly distributed states is described by a generalization of hydrodynamics, which takes into consideration the infinitely many local conservation laws. At the scale of non-diffusive, ballistic transport, the main feature of this generalized hydrodynamics (GHD) is the existence of continuity equations for the occupation number of each mode. Further two projects are presented in this Thesis which are concerned with GHD.I n Chapter 4, a numerical method is presented for computing the evolution of entanglement entropies within GHD. This method combines the quasiparticle picture of entanglement with a molecular dynamics simulation of GHD, the so-called fleagas. Numerical evidence is shown that this method reproduces previously known analytical results in GHD. Its versatility is demonstrated by computing another entanglement-related quantity, the mutual information, for which no analytical results are available. Finally, in Chapter 5, the GHD framework is generalized to a system of spin-1/2 fermions with repulsive contact interaction, the Yang–Gaudin model. It is shown that at low temperatures, the dynamics of the system is described by the superposition of two uncoupled, inhomogeneous conformal field theories. This fact can be interpreted as a separation of the spin and the charge degrees of freedom of the system.
This thesis investigates quantum-many body systems out of equilibrium in a variety of physically relevant and intriguing settings. It provides an overview of the recent developments in this field and based on the author’s recent review illuminates many interesting connections and summarises open questions. In this general framework, the various results of the author are embedded. While they often rely on advanced mathematics, great care is taken to present them in an intuitive way and most technical material is discussed separately in appendices. In the context of equilibration, work created as part of this thesis is able to capture the Gaussification of correlated initial states for a large class of free models. This allows to significantly extend equilibration results of free models, which are of crucial importance, as they provide reasonable relaxation time scales and immediately connect to our physical intuition in terms of ballistic spreading. Building on two surprising connections between static and dynamic features of Hamiltonians, which in general are hard to obtain, the results of this thesis further greatly contribute to the recent debate on the true nature of quantum many- body localisation. For those interacting models, we derive the localised structure of eigenstates from a dynamical suppression of information propagation on the low-energy sector. Further, we show how a non-degenerate spectrum, indicating the presence of interactions, and the existence of an approximately local constant of motion are sufficient to show information propagation, if arbitrary energies are allowed. In the context of quantum phase transitions, we continue such connections between static and dynamic features. We specifically investigate the Mott-superfluid transition of the Bose-Hubbard model and ask to what extent it has a universal dynamical signature. In a joint experimental, numerical and analytical effort, complex behaviour of these dynamics is uncovered, thus challenging the common believe that the Kibble-Zurek mechanism is sufficient to fully capture such transitions. We embed such experimental investigation in the recent debate on quantum simulators, which give the exciting perspective of solving notoriously and even provably hard problems efficiently in the laboratory. For such devices, we argue that the final read-out of the results is an important out- standing problem to be solved. We approach this issue in the case of a continuous quantum field, in which the notion of reconstructing the quantum state is even conceptually unclear. Based on tensor network methods, we demonstrate that in an experiment of ultra-cold atoms in a continuous setup, states can nevertheless efficiently be obtained. Thus, this thesis constitutes not only an important review of the field of quantum many-body systems out of equilibrium, but, using advanced mathematical and numerical tools, has significantly contributed to our understanding of various important out- standing questions in interacting many-body systems.
The study of non-equilibrium behaviors of quantum matter is one of the research directions at the frontiers of condensed matter physics and, in general, of theoretical and experimental physics. Understanding the mechanisms which prevent many quantum particles from reaching thermal equilibrium can ultimately shed light on the foundations of statistical physics and provide insight into the anomalous long-time collective behaviors observed in experiments. In this thesis, we present new theoretical developments in this direction. In the first part of the thesis we discuss the non-equilibrium dynamics of the two-dimensional quantum Ising model, the prototypical model of magnetism. In particular, we investigate the dynamics of interfaces separating different spin domains. In some regimes, we can provide an exact description of their unitary time evolution, but we also discuss the robustness of our predictions away from the exactly-solvable limit. Moreover, we address the same problem in the presence of disorder and provide a hydrodynamic interpretation of our findings. In the second part of the thesis, we investigate the localization properties of disordered quantum systems. In particular, we develop a renormalization-group picture of the Anderson localization transition in generic geometries. This provides a new framework to understand the phase diagram of such systems and their critical behavior beyond the usual finite-size scaling analysis. We also provide further insights on this topic both by studying resonances and their connection with observables and by employing unsupervised non-parametric learning techniques to probe disordered systems in an experimentally accessible form. The results presented in this thesis constitute a step forward in understanding the non-equilibrium phases and regimes of quantum matter and they pave the way to interesting future developments.
In the last decade the non-equilibrium dynamics of isolated quantum systems has raised a huge interest for both its conceptual and practical implications. Simulation of quantum dynamics represents by now a very active research effort, with highly controllable quantum devices composed of hundreds of qubits already being realized in numerous setups. These offer the unprecedented possibility to access real-time quantum dynamics and strongly correlated quantum many-body states, opening numerous avenues not only for future quantum technologies, but also for long-standing theoretical problems. One of the most sought-after applications is the computation of real-time evolution in theories of fundamental interactions. Since Wilson's proposal in 1974, the equilibrium phase diagrams of lattice gauge theories have been successfully investigated with Monte Carlo simulations in a wide range of parameters. Nonequilibrium properties, on the other hand, cannot be accessed with conventional techniques. Thanks to the recent advances in numerics (e.g. tensor networks) and experiments (quantum simulation), new promising tools are now available for tackling these questions. Remarkable works have already demonstrated that applying these techniques to simple theories, such as one-dimensional quantum electrodynamics, can help to detect the signatures of confinement and pair creation. Despite these impressive achievements, investigating the non-equilibrium dynamics of gauge theories with quantum simulators turned out to be extremely challenging. This is mostly due to the complex constrained dynamics of these theories, which make their simulation with realistic atomic platforms particularly difficult. At the same time, the dynamical constraints can play a key role in determining novel interesting phenomena and peculiar non-equilibrium properties. These phenomena -- which include, for example, quantum scars, fractons, Hilbert space fragmentation -- can lead to anomalous dynamics, eventually evading the paradigm of thermalization, which predicts the relaxation of physical observable to their thermodynamic values. The aim of this thesis is twofold: first, we want to address through a combination of numerical and analytical tools some of the open questions regarding the non-equilibrium dynamics of quantum many-body systems with constrained Hilbert spaces, including lattice gauge theories. In addition, the work reported here is part of a large effort of the community towards the quantum simulation of complex phenomena of condensed matter and high-energy physics.
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.
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.