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.
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.
In this PhD thesis we investigate some properties of one-dimensional quantum systems, focusing on two important aspects of integrable models: Their entanglement properties at equilibrium and their dynamical correlators after a quantum quench. The first part of the thesis will therefore be devoted to the study of the entanglement entropy in one-dimensional integrable systems, with a special focus on the XYZ spin-1/2 chain, which, in addition to being integrable, is also an interacting theory with non perturbative solutions. We derive its bipartite Renyi entropies in the thermodynamic limit and its behaviour in different phases and for different values of the mass-gap is analysed, both analytically and numerically. In particular it is worth mentioning that the numerical analysis of the entropies of the XYZ model presented in this dissertation has never been published in literature. In the second part of the thesis we study the dynamics of correlators after a quantum quench, preparing the system in a squeezed coherent initial state. The emphasis will be on the Transverse Field Ising Chain and the O(3) non-linear sigma model, the latter studied by means of the semi-classical theory, the former by a form-factor approach. Moreover in the last chapter we outline a general result about the dynamics of correlation functions of local observables after a quantum quench. In particular we show that if there are not long-range interactions in the final Hamiltonian, then the dynamics of the model at long times (non equal- time correlations) is described by the same statistical ensemble that describes its statical properties (equal-time correlations). For the Transverse Field Ising Chain this result means that its dynamics is described by a Generalized Gibbs Ensemble.
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.
The primary area, on which we focus in the thesis, are the out-of-equilibrium properties of quantum and classical systems. The playground we use to study those properties, are interacting integrable and solitonic systems, whose structure enables us to perform analytical calculations.
In the context of the out-of-equilibrium physics we focus on two questions. The first problem is related to the set of stationary states in quantum systems. The question is especially important in the context of the homogeneous quantum quenches, since it makes possible the calculations of expectation values of local observables. The building blocks of stationary states are quasilocal conserved quantities which constitute the information that is preserved under the time evolution. In the first part of the thesis we deal with the construction of quasilocal conserved quantities in one dimensional anisotropic Heisenberg model and bosonic models. Quasilocal charges are put in the context of the standard theory of integrability, which in turn allows for the calculation of physical quantities.
The second question we deal with are the transport properties of one dimensional systems. The connection between the local integrals of motion and ideal transport was established in 90's. Here we extend the connection to the diffusive transport. This enables us to obtain explicit lower bound on diffusion constant in Heisenberg model, which proves that the Heisenberg model is not an insulator at finite temperatures.
In the last chapter we deal with the transport properties of classical celular automata. The first example we consider is the gas of charged hard core interacting particles, and the second one the gas of interacting solitons. Explicit solution of the dynamics of local observables allows for the explicit calculation of transport coeffcients, and the solution of the inhomogeneous quench problem. Despite their simplicity, the models exhibit wide range of transport phenomena ranging from insulating, through diffusive to ballistic transport.
In this thesis, we investigated the evolution of low-dimensional quantum many-body systems through the lenses of their dynamics. The main tools we used are unitary dynamics and kinetic theory. The interest in this class of systems derives from the fact that they are dominated by the effects of interactions between fundamental constituents. In one dimension, the situation is even more dramatic: even the smallest interaction leads to a description that cannot be cast into a small perturbation of the initial model. Even though, a particular set of systems still allows for exact solutions; hence, they are deemed integrable and they allow for an infinite number of conserved currents. A simple way, that can also be performed experimentally, to drive many-body systems outside of the equilibrium is called quantum quench. It consists in switching instantaneously the parameters of the theory after the system is at the equilibrium and measuring an observable after the driving has occurred. Our contributions to the field amount to providing a full characterization of the evolution of the observables of an integrable system after a global quantum quench. We were also able to extend these results, valid in principle for relativistic theories, to the Lieb-Liniger model, prominent for the description of experimental setups, by the means of non-relativistic limit. For two-dimensional models, interactions are not changing the system from its foundations, although there are unexpected effects. A paradigmatic example is graphene: in a two-dimensional system where carriers are charged and massless, namely their dispersion relation does not have any gap and particles and holes can be freely excited from the thermal state, it is predicted that they are interacting via Coulomb force hundreds of times more than normal electrons in void. Hence, the system behaves in a completely different way from electrons in quantum electrodynamics: for instance, there are additional modes which are called plasmons that are also massless but neutral. They represent collective excitations of the model, as oscillations in the local density of charge. This system is not integrable but still strongly interacting; solving it as we do for one-dimensional systems is not possible. The role of conservation laws is a common theme of this thesis. In the example above exposed, total charge, momentum and energy are conserved in the pristine model, but in presence of impurities and defects they dissipate. The evolution of such quantities can though be derived exactly by the means of hydrodynamics. We were able to derive the equations that rule their evolution in time for any electronic system with two bands touching at the Dirac point. They are generalizations of the celebrated Navier-Stokes equations, and they differ from them only outside the compressible regime. We are able to show that the transport properties in the hydrodynamic regime strongly differ from the usual ballistic picture popularized by Drude and his model; this is captured for instance by violations of the Wiedemann-Franz ratio. Moreover, we prove that the role of interaction is prominent in the viscosity of these fluids.
Recently, the possibility of investigating single molecule, or single spin observables, as well as the necessity of a better understanding of the mechanisms underlying quantum dynamics in order to obtain nanoscale devices and nanostructered materials suitable for quantum computing tasks, have revived the interest in foundational aspects of quantum statistical mechanics. This thesis aims to give a contribution to this field by re-considering the statistical characterization of a quantum system at the light of some paradigmatic changes in our understanding of quantum theory which have taken place in the last two decades. In particular the impressive development of quantum information theory has changed the perceptions of quantum entanglement: for a long time it has been considered a somewhat paradoxical property of the matter at the atomic scale, but now it is regarded as an essential and ubiquitous phenomenon whose consequences are affecting the very macroscopic world that we experience. Still the decoherence program has brought out the importance of considering a quantum system together with its environment in order to clarify some key aspects of quantum dynamics. Thus, we start from the idea that quantum correlations are ubiquitous and somewhat uncontrollable in systems with many degrees of freedom which are typically considered in statistical mechanics. As a consequence we assume the standpoint that quantum statistical mechanics has not to be based on the underlying idea of a collection of many, independent quantum systems but rather it has to emerge at the level of a global wavefunction (pure state) which describes the system as well as its environment as a whole. In order to investigate the consequences of these assumptions we study the equilibrium distribution of an isolated quantum system. This is defined, in analogy with the ergodic foundations of classical statistical mechanics, on the basis of the time evolution of the quantum state. Then, we study the emergence of thermodynamic properties in a quantum system by studying the probability distribution of some function of interests, as the entropy and the equilibrium state of a subsystem, on Ensembles of Pure States. Such a probability distribution is derived from the geometry of the Hilbert space, and the theoretical tools suitable for its characterization are developed. On the one hand we perform a numerical sampling of the ensemble distributions by employing Monte Carlo techniques, on the other hand simpler analytical approximation of the geometrical distributions are derived by means of a maximum entropy principle. Model systems composed of an ensemble of spins are chosen to illustrate the salient features which emerge from the developed theoretical framework: the main point is that the Ensemble Distributions of “thermodynamic observables” (entropy or equilibrium state of a subsystem) are sharply peaked around a typical value. From the analysis it emerges that each of the overwhelming majority of the wavefunctions which has appreciable weight in the considered ensemble, is characterized by the same value of the “macroscopic” functions. This is a striking evidence of the “typicality” of these properties. In the essence, our impossibility to know the state of the system in detail does not matter, just for the remarkable fact that almost all quantum states behave essentially in the same way. By virtue of this typicality the study of the behaviour of the typical values of the thermodynamic function become meaningful. Notably, under certain conditions, one recovers the results of standard statistical mechanics, that is, the equilibrium average of the state of a subsystem can be cast in the Boltzmann canonical form at the temperature given by the usual thermodynamical relation . In the second part of the thesis we consider the dynamical aspects of the equilibrium state of a subsystem interacting with its environment. The fluctuations around the equilibrium average critically depends on the entanglement between the system and the environment and on the form of the interaction Hamiltonian. The connection between the dynamics of the fluctuations of an observable at the equilibrium and the relaxation toward the equilibrium from a “non typical” initial value is also investigated with the aid of simple model systems. The study presented in this thesis was partly motivated by a critical analysis of the statistical methods available for the theoretical modelling of magnetic resonance experiments. One of these, the Stochastic Liouville Equation, has been employed in a work completed during the first year of my Ph.D. program in order to interpret some feature of a two dimensional electron spin resonance experiment, [Fresch B., Frezzato D., Moro G. J., Kothe G., Freed J. H.; J. Phys. Chem. B., 110, 24238, (2006)].
Among the few lessons that Physics teaches us on a daily basis, one in particular is hard to miss: thermalization processes are ubiquitous. With the same degree of certainty, one can predict an apple to fall onto the ground if its stalk is cut, or a cup of tea to be found cold if left on the table for too long. In this respect, generic classical and quantum mechanical systems seem to behave in the same way. This is also one of the main difficulties in the experimental observation of the bizarre effects predicted from the quantum theory, as several of these are known to disappear at finite temperatures. It has been known for a long time that thermalization is associated with “chaotic” behavior at the microscopic level, but recently systematic efforts, both theoretical and experimental, have allowed us to understand its mechanisms at an unprecedented level of accuracy. Indeed, it has been realized that thermalization generally occurs also in isolated systems, where the absence of interactions with the environment allows for first principle theoretical investigations1. In this case thermalization takes place locally, as the whole system acts as a thermal bath for its own subsystems. Theoretical research has been motivated and paralleled by exciting experimental progress in cold-atom physics, which has provided robust, nearly ideal realizations of several theoretical models2. As an established piece of knowledge, recent research has now confirmed that two outstanding exceptions actually exist in the quantum realm, where thermalization is prevented to occur. These are many-body localized systems3, where disorder plays a crucial role, and integrable ones, which are protected by the existence of higher conservation laws. By their own nature, these systems display exceptional non-equilibrium features, which are not washed out by the onset of thermalization as relaxation occurs. At the core of the present thesis lie the remarkable properties of integrable systems out of equilibrium, with a special attention to physical effects which could be observed in cold-atom realizations. The work collected here is part of a large theoretical effort4 which in the past decade has focused on the study of relatively simple protocols to bring a quantum system out of equilibrium, such as the so-called quantum quench. In these “theoretical laboratories” it has been possible to provide quantitative predictions which helped us to develop an intuition on general questions regarding non-equilibrium and thermalization processes. Furthermore, in some cases these studies have highlighted interesting physical effects which could be directly probed experimentally within cold-atom settings. The contribution of the present thesis is two-fold. On the one hand, we have developed new technical tools for the study of integrable systems out of equilibrium, and for the computation of measurable physical quantities such as correlation functions. The techniques employed are mainly analytical and rooted within the mathematical structures of integrability. On the other hand, we have singled out physically relevant nonequilibrium situations where exotic, non-thermal states of matter emerge after relaxation occurs, and provided quantitative analytical predictions in these cases. A comprehensive discussion of the motivations and results of our work will be presented in Chapter 1, where we provide a complete overview of this thesis and discuss the organization of its content.