This cumulative doctoral thesis, based on the works [P1–P5], addresses various topical concepts as well as foundational questions in the nonequilibrium physics of quantum many-body systems. Publication [P1] concerns itself with the question of time’s arrow in quantum mechanics. Based on the recursion method, we attribute autocorrelation functions of few-body observables at infinite temperature (almost) monotonously decreasing functions, constructed from temporal derivatives of the correlation function that indicate a directed approach to equilibrium. We numerically find that in chaotic models, these constructions comprise a notion of simplicity of the underlying dynamics. Addressing the ubiquity of simple relaxation dynamics in otherwise complex chaotic systems, in Ref. [P3] we introduce a scheme based on Lanczos coefficients to approximate correlation functions in the thermodynamic limit. Examining several pertinent condensed-matter type models and benchmarking with state-of-the-art numerics, we find the method to yield good results already for a low number of included Lanczos coefficients, provided that the latter exhibit a smooth and increasing structure. We underpin these findings with analytical arguments on the method’s applicability. Ref. [P2] provides an estimate on the equilibration time of such infinite-temperature correlation functions in the thermodynamic limit. Similarly to Ref. [P3], we find the estimate to capture the equilibration times very reasonably, given that the corresponding Lanczos coefficients exhibit a smooth growth. Substantiated by analytical considerations, the convergence of the estimate implies that equilibration timescales of the relaxation processes are finite and conceivably short. In Ref. [P5] we scrutinize the interpretative power of the Lanczos coefficients concerning the system’s integrability or chaoticity, operator growth and the indicativeness of physical behavior. To this end, we analytically deter- mine the Lanczos coefficients for a Majorana observable in a strictly quadratic model and illustrate that within this model any structure of the Lanczos coefficients may be realised. Different to the other works, Ref. [P4] is not cast in the language of the recursion method. Focusing on the full eigenstate thermalization hypothesis, it explores the connection between thermalization dynamics of multi-point correlation functions and the scaling of thermal free cumulants in setups mimicking a system-bath situation.
This cumulative doctoral thesis, based on the works [P1–P5], addresses various topical concepts as well as foundational questions in the nonequilibrium physics of quantum many-body systems. Publication [P1] concerns itself with the question of time’s arrow in quantum mechanics. Based on the recursion method, we attribute autocorrelation functions of few-body observables at infinite temperature (almost) monotonously decreasing functions, constructed from temporal derivatives of the correlation function that indicate a directed approach to equilibrium. We numerically find that in chaotic models, these constructions comprise a notion of simplicity of the underlying dynamics. Addressing the ubiquity of simple relaxation dynamics in otherwise complex chaotic systems, in Ref. [P3] we introduce a scheme based on Lanczos coefficients to approximate correlation functions in the thermodynamic limit. Examining several pertinent condensed-matter type models and benchmarking with state-of-the-art numerics, we find the method to yield good results already for a low number of included Lanczos coefficients, provided that the latter exhibit a smooth and increasing structure. We underpin these findings with analytical arguments on the method’s applicability. Ref. [P2] provides an estimate on the equilibration time of such infinite-temperature correlation functions in the thermodynamic limit. Similarly to Ref. [P3], we find the estimate to capture the equilibration times very reasonably, given that the corresponding Lanczos coefficients exhibit a smooth growth. Substantiated by analytical considerations, the convergence of the estimate implies that equilibration timescales of the relaxation processes are finite and conceivably short. In Ref. [P5] we scrutinize the interpretative power of the Lanczos coefficients concerning the system’s integrability or chaoticity, operator growth and the indicativeness of physical behavior. To this end, we analytically deter- mine the Lanczos coefficients for a Majorana observable in a strictly quadratic model and illustrate that within this model any structure of the Lanczos coefficients may be realised. Different to the other works, Ref. [P4] is not cast in the language of the recursion method. Focusing on the full eigenstate thermalization hypothesis, it explores the connection between thermalization dynamics of multi-point correlation functions and the scaling of thermal free cumulants in setups mimicking a system-bath situation.
We report an implementation of the recursion method that addresses quantum many-body dynamics in the nonperturbative regime. The method essentially amounts to constructing a Lanczos basis in the space of operators and solving coupled Heisenberg equations in this basis. The reported implementation has two key ingredients: a computer-algebraic routine for the symbolic calculation of nested commutators and a procedure to extrapolate the sequence of Lanczos coefficients according to the universal operator growth hypothesis. We apply the method to calculate infinite-temperature correlation functions for spin-1/2 systems on one- and two-dimensional lattices. In two dimensions the accessible timescale is large enough to essentially embrace the relaxation to equilibrium. The method allows one to accurately calculate transport coefficients. As an illustration, we compute the diffusion constant for the transverse-field Ising model on a square lattice.
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.
Technological and scientific advances have given rise to an era in which coherent quantum-mechanical phenomena can be probed and experimentally-realised over unprecedented timescales in condensed matter physics. In turn, scientific interest in non-equilibrium dynamics and irreversibility signatures of thermodynamics has taken place in recent decades, particularly in relation to cold-atom platforms and thermoelectric devices. In this PhD thesis I summarise some of the most important results obtained over the duration of my PhD, on the topic of thermodynamics involving interacting many-body quantum systems. The topics of discussion encompass three main themes: spin/particle transport in non-integrable systems, explorations in eigenstate thermalisation and finite-temperature transport in autonomous thermal machines. By questioning the effect of local integrability-breaking perturbations, I describe the subtle effects that may rise from conservation laws and their connection to linear response transport using microscopic models. Then, I describe thermalisation in the context of the eigenstate thermalisation hypothesis and its consequence to multipartite entanglement and high-order correlation functions. Finally, motivated by the necessity of describing quantum thermal machines in the finite-temperature regime, I introduce a novel tensor-network based method to investigate thermodynamic properties of autonomous machines that employ interacting many-body systems as working media.
We study different nonequilibrium phenomena of isolated quantum systems ranging from few- to many-body interacting bosons. Firstly, we have suggested the dynamics of the center-of-mass motion to sensitively detect unconverged numerical many-body dynamics in potential with separable quantum motion of the center of mass. As an example, we consider the time evolution of attractive bosons in a homogenous background and use it to benchmark a specific numerical method based on variational multimode expansion of the many-body wave function - the Multicon gurational time-dependent Hartree for bosons (MCTDHB). We demonstrate that the simplified convergence criterion based on a threshold value for the least occupied mode function fails to assure qualitatively correct result while our suggested convergence test based on the center-of-mass motion correctly detects the deviation of numerical results from the exact results. Recent technological progress in manipulating low-entropy quantum states has motivated us to study the phenomenon of interaction blockade in bosonic systems. We propose an experimental protocol to observe the expected bosonic enhancement factor in this blockade regime. Specifically, we suggest the use of an asymmetric double-well potential constructed by superposition of multiple optical tweezer laser beams. Numerical simulations using the MCTDHB method predict that the relevant states and the expected enhancement factor can be observed. In the second half of the thesis, we have investigated the onset of quantum thermalization in a two-level generalization of the Bose-Hubbard dimer. To this end, the relaxation dynamics following a quench is studied using two numerical methods: We study different nonequilibrium phenomena of isolated quantum systems ranging from few- to many-body interacting bosons. Firstly, we have suggested the dynamics of the center-of-mass motion to sensitively detect unconverged numerical many-body dynamics in potential with separable quantum motion of the center of mass. As an example, we consider the time evolution of attractive bosons in a homogenous background and use it to benchmark a specific numerical method based on variational multimode expansion of the many-body wave function - the Multicon gurational time-dependent Hartree for bosons (MCTDHB). We demonstrate that the simplified convergence criterion based on a threshold value for the least occupied mode function fails to assure qualitatively correct result while our suggested convergence test based on the center-of-mass motion correctly detects the deviation of numerical results from the exact results. Recent technological progress in manipulating low-entropy quantum states has motivated us to study the phenomenon of interaction blockade in bosonic systems. We propose an experimental protocol to observe the expected bosonic enhancement factor in this blockade regime. Specifically, we suggest the use of an asymmetric double-well potential constructed by superposition of multiple optical tweezer laser beams. Numerical simulations using the MCTDHB method predict that the relevant states and the expected enhancement factor can be observed. In the second half of the thesis, we have investigated the onset of quantum thermalization in a two-level generalization of the Bose-Hubbard dimer. To this end, the relaxation dynamics following a quench is studied using two numerical methods: (1) full quantum dynamics and (2) semiclassical phase-space method. We rely on arguments based on the eigenstate thermalization hypothesis (ETH), quantum chaos as seen from the distribution of level spacings, and the concept of chaotic eigenstates in demonstrating equilibration dynamics of local observables in the system after an integrability-breaking quench. The same issue on quantum thermalization can be viewed from a different perspective using semiclassical phase-space methods. In particular, we employ the truncated Wigner approximation (TWA) to simulate the quantum dynamics. In this case, we show that the marginal distributions of the individual trajectories which sample the initial Wigner distribution are in good agreement with the corresponding microcanonical distribution.
This thesis addresses different aspects of the dynamics of isolated quantum many-body systems far from equilibrium, ranging from transient phenomena at short times to questions about relaxation and equilibration at late times. The first part deals with nonequilibrium critical phenomena. Dynamical quantum phase transitions are characterized in a two-dimensional spin system by studying the distribution of Fisher zeros in the complex time plane. Moreover, a new topologically driven transition in the nonequilibrium steady state of quenched Chern insulators is investigated, resulting in the identification of universal nonanalytic behavior of the Hall conductance as a function of the quench parameter. The second part of the thesis comprises a study of transient dynamics of the transverse-field Ising model in two and three dimensions based on a perturbative derivation of classical networks to encode the time-evolved wave function; this construction is also applicable to represent the state as an artificial neural network. In the last part of the thesis a possible definition of the irreversibility of quantum many-body dynamics is introduced, which is based on the echo dynamics of observables under imperfect effective time reversal. Echo dynamics is studied in generic and integrable quantum systems as well as in the semiclassical limit. The relation to out-of-time-order correlators is explored, which were recently introduced to probe the quantum butterfly effect.
The recursion method, which solves coupled Heisenberg equations in a Lanczos operator basis, has recently emerged as a powerful nonperturbative tool for computing dynamical correlation functions in strongly correlated two- and three-dimensional quantum many-body systems. Motivated by this success, we investigate whether the method can be extended to expectation values of observables following a quantum quench. We find that such an extension encounters an obstacle absent in the computation of dynamical correlation functions. The latter are fully determined by the Lanczos coefficients b_n b n , which in generic systems exhibit universal behavior, enabling reliable extrapolation from the first few dozens of explicitly computed coefficients. In contrast, quench dynamics additionally requires “quench coefficients” c_n c n , defined as overlaps of Lanczos basis operators with the initial state. We show that, unlike the Lanczos coefficients, the quench coefficients exhibit no universal structure and cannot be reliably extrapolated, thereby limiting the time up to which the method yields accurate results. The behavior of quench coefficients is highly state-dependent, ranging from decaying to irregular or even growing sequences; typically, the less regular the sequence c_n c n , the shorter the accessible timescale. Nevertheless, for favorable initial states, the method remains competitive with state-of-the-art approaches. Moreover, its symbolic implementation allows a single computation to be reused across different Hamiltonian parameters and initial states, making it particularly advantageous in studies requiring extensive scans over Hamiltonian parameters or initial states.