We suggest a method to compute approximations to temporal correlation functions of few-body observables in chaotic many-body systems in the thermodynamic limit based on the respective Lanczos coefficients. Given the knowledge of these Lanczos coefficients, the method is very cheap. Usually, accuracy increases with more Lanczos coefficients taken into account; however, we numerically find and analytically argue that the convergence is rather quick if the Lanczos coefficients exhibit a smoothly increasing structure. For pertinent examples, we compare with data from dynamical typicality computations for large but finite systems and find good agreement if few Lanczos coefficients are taken into account. From the method, it is evident that in these cases the correlation functions are well described by a low number of damped oscillations.
We suggest a method to compute approximations to temporal correlation functions of few-body observables in chaotic many-body systems in the thermodynamic limit based on the respective Lanczos coefficients. Given the knowledge of these Lanczos coefficients, the method is very cheap. Usually, accuracy increases with more Lanczos coefficients taken into account; however, we numerically find and analytically argue that the convergence is rather quick if the Lanczos coefficients exhibit a smoothly increasing structure. For pertinent examples, we compare with data from dynamical typicality computations for large but finite systems and find good agreement if few Lanczos coefficients are taken into account. From the method, it is evident that in these cases the correlation functions are well described by a low number of damped oscillations.
We study the equilibration times T_{eq} of local observables in quantum chaotic systems by considering their autocorrelation functions. Based on the recursion method, we suggest a scheme to estimate T_{eq} from the corresponding Lanczos coefficients that is expected to hold in the thermodynamic limit. We numerically find that, if the observable eventually shows smoothly growing Lanczos coefficients, a finite number of the former is sufficient for a reasonable estimate of the equilibration time. This implies that equilibration occurs on a realistic time scale much shorter than the life of the Universe. The numerical findings are further supported by analytical arguments.
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 establish rigorously the scaling properties of the Lanczos process applied to an arbitrary extensive many-body system which is carried to convergence n→∞ and the thermodynamic limit N→∞ taken. In this limit the solution for the limiting Lanczos coefficients are found exactly and generally through two equivalent sets of equations, given initial knowledge of the exact cumulant generating function. The measure and the orthogonal polynomial system associated with the Lanczos process in this regime are also given explicitly. Some important representations of these Lanczos functions are given, including Taylor series expansions, and theorems controling their general properties are proven.
We review various theoretical methods that have been used in recent years to calculate dynamical correlation functions of many-body systems. Time-dependent correlation functions and their associated frequency spectral densities are the quantities of interest, for they play a central role in both the theoretical and experimental understanding of dynamic properties. In particular, dynamic correlation functions appear in the fluctuation-dissipation theorem, where the response of a many-body system to an external perturbation is given in terms of the relaxation function of the unperturbed system, provided the disturbance is small. The calculation of the relaxation function is rather difficult in most cases of interest, except for a few examples where exact analytic expressions are allowed. For most of systems of interest approximation schemes must be used. The method of recurrence relation has, at its foundation, the solution of Heisenberg equation of motion of an operator in a many-body interacting system. Insights have been gained from theorems that were discovered with that method. For instance, the absence of pure exponential behavior for the relaxation functions of any Hamiltonian system. The method of recurrence relations was used in quantum systems such as dense electron gas, transverse Ising model, Heisenberg model, XY model, Heisenberg model with Dzyaloshinskii-Moriya interactions, as well as classical harmonic oscillator chains. Effects of disorder were considered in some of those systems. In the cases where analytical solutions were not feasible, approximation schemes were used, but are highly model-dependent. Another important approach is the numericallly exact diagonalizaton method. It is used in finite-sized systems, which sometimes provides very reliable information of the dynamics at the infinite-size limit. In this work, we discuss the most relevant applications of the method of recurrence relations and numerical calculations based on exact diagonalizations. The method of recurrence relations relies on the solution to the coefficients of a continued fraction for the Laplace transformed relaxation function. The calculation of those coefficients becomes very involved and, only a few cases offer exact solution. We shall concentrate our efforts on the cases where extrapolation schemes must be used to obtain solutions for long times (or low frequency) regimes. We also cover numerical work based on the exact diagonalization of finite sized systems. The numerical work provides some thermodynamically exact results and identifies some difficulties intrinsic to the method of recurrence relations.
We present a method to compute dynamical correlation functions in quantum many-body systems at zero temperature. The ground state of a finite system is evaluated exactly by a modified Lanczos method and from it the real-frequency correlation functions are obtained by a projective technique for the memory-function formalism. We apply the method to the density-density correlation function of a one-dimensional spinless fermion system. Exact results are recovered in the noninteracting limit; in the strong-coupling limit our results are compared with Monte Carlo simulation and analytical approximations.
We review various theoretical methods that have been used in recent years to calculate dynamical correlation functions of many-body systems. Time-dependent correlation functions and their associated frequency spectral densities are the quantities of interest, for they play a central role in both the theoretical and experimental understanding of dynamic properties. The calculation of the relaxation function is rather difficult in most cases of interest, except for a few examples where exact analytic expressions are allowed. For most of systems of interest approximation schemes must be used. The method of recurrence relation has, at its foundation, the solution of Heisenberg equation of motion of an operator in a many-body interacting system. Insights have been gained from theorems that were discovered with that method. For instance, the absence of pure exponential behavior for the relaxation functions of any Hamiltonian system. The method of recurrence relations was used in quantum systems such as dense electron gas, transverse Ising model, Heisenberg model, XY model, Heisenberg model with Dzyaloshinskii-Moriya interactions, as well as classical harmonic oscillator chains. Effects of disorder were considered in some of those systems. In the cases where analytical solutions were not feasible, approximation schemes were used, but are highly model-dependent. Another important approach is the numerically exact diagonalization method. It is used in finite-sized systems, which sometimes provides very reliable information of the dynamics at the infinite-size limit. In this work, we discuss the most relevant applications of the method of recurrence relations and numerical calculations based on exact diagonalizations.