Possible strategies to describe analytically the dynamics of many-body quantum systems out of equilibrium include the use of solvable models and of full random matrices. None of the two approaches represent actual realistic systems, but they serve as references for the studies of these ones. We take the second path and obtain analytical expressions for the survival probability, density imbalance, and out-of-time-ordered correlator. Using these findings, we then propose an approximate expression that matches very well numerical results for the evolution of realistic finite quantum systems that are strongly chaotic and quenched far from equilibrium. In the case of the survival probability, the expression proposed covers all different time scales, from the moment the system is taken out of equilibrium to the moment it reaches a new equilibrium. The realistic systems considered are described by one-dimensional spin-1/2 models.
Possible strategies to describe analytically the dynamics of many-body quantum systems out of equilibrium include the use of solvable models and of full random matrices. None of the two approaches represent actual realistic systems, but they serve as references for the studies of these ones. We take the second path and obtain analytical expressions for the survival probability, density imbalance, and out-of-time-ordered correlator. Using these findings, we then propose an approximate expression that matches very well numerical results for the evolution of realistic finite quantum systems that are strongly chaotic and quenched far from equilibrium. In the case of the survival probability, the expression proposed covers all different time scales, from the moment the system is taken out of equilibrium to the moment it reaches a new equilibrium. The realistic systems considered are described by one-dimensional spin-1/2 models.
We study numerically and analytically the quench dynamics of isolated many-body quantum systems. Using full random matrices from the Gaussian orthogonal ensemble, we obtain analytical expressions for the evolution of the survival probability, density imbalance, and out-of-time-ordered correlator. They are compared with numerical results for a one-dimensional-disordered model with two-body interactions and shown to bound the decay rate of this realistic system. Power-law decays are seen at intermediate times, and dips below the infinite time averages (correlation holes) occur at long times for all three quantities when the system exhibits level repulsion. The fact that these features are shared by both the random matrix and the realistic disordered model indicates that they are generic to nonintegrable interacting quantum systems out of equilibrium. Assisted by the random matrix analytical results, we propose expressions that describe extremely well the dynamics of the realistic chaotic system at different time scales.
Despite its importance to experiments, numerical simulations, and the development of theoretical models, self-averaging in many-body quantum systems out of equilibrium remains underinvestigated. Usually, in the chaotic regime, self-averaging is taken for granted. The numerical and analytical results presented here force us to rethink these expectations. They demonstrate that self-averaging properties depend on the quantity and also on the time scale considered. We show analytically that the survival probability in chaotic systems is not self-averaging at any time scale, even when evolved under full random matrices. We also analyze the participation ratio, R\'enyi entropies, the spin autocorrelation function from experiments with cold atoms, and the connected spin-spin correlation function from experiments with ion traps. We find that self-averaging holds at short times for the quantities that are local in space, while at long times, self-averaging applies for quantities that are local in time. Various behaviors are revealed at intermediate time scales.
We study analytically and numerically the non-equilibrium dynamics of an isolated interacting many-body quantum system following a random quench. We model the system Hamiltonian by Embedded Gaussian Orthogonal Ensemble (EGOE) of random matrices with one plus few-body interactions for fermions. EGOE are paradigmatic models to study the crossover from integrability to chaos in interacting many-body quantum systems. We obtain a generic formulation, based on spectral variances, for describing relaxation dynamics of survival probabilities as a function of rank of interactions. Our analytical results are in good agreement with numerics.
Recent experimental achievements in controlling ultracold gases in optical lattices open a new perspective on quantum many-body physics. In these experimental setups, it is possible to study coherent time evolution of isolated quantum systems. These dynamics reveal new physics beyond the low-energy properties that are usually relevant in solid-state many-body systems. In this paper, we study the time evolution of antiferromagnetic order in the Heisenberg chain after a sudden change of the anisotropy parameter, using various numerical and analytical methods. As a generic result, we find that the order parameter, which can show oscillatory or non-oscillatory dynamics, decays exponentially except for the effectively non-interacting case of the XX limit. For weakly ordered initial states, we also find evidence for an algebraic correction to the exponential law. The study is based on numerical simulations using a numerical matrix product method for infinite system sizes (iMPS), for which we provide a detailed description and an error analysis. Additionally, we investigate in detail the exactly solvable XX limit. These results are compared to approximative analytical approaches including an effective description by the XZ model as well as by mean-field, Luttinger-liquid and sine-Gordon theories. The comparison reveals which aspects of non-equilibrium dynamics can, as in equilibrium, be described by low-energy theories and which are the novel phenomena specific to quantum quench dynamics. The relevance of the energetically high part of the spectrum is illustrated by means of a full numerical diagonalization of the Hamiltonian.
In this work, we study the non-equilibrium dynamics of low-dimensional strongly correlated quantum systems, which in our case are generated by sudden and time-dependent quantum quenches. Such non-equilibrium dynamics are a particularly difficult problem, because the time evolution of the initial many-body quantum state, which is governed by the post-quench Hamiltonian, is highly non-trivial. On top of that, applying an operator to this many-body quantum state in order to calculate its expectation values, typically creates an intricate set of excitations. Throughout the thesis, we furthermore contrast the results that we have obtained using exact numerical methods for two microscopic models, namely the Heisenberg and Fermi-Hubbard chains, with results derived from the Luttinger liquid theory, which, in equilibrium, describes the universal behavior of the two models.In the second part of the thesis, we reduce the dimensionality even further and employ complementary numerical and analytical tools to examine the unusual non-equilibrium behavior of a ring-shaped quantum impurity hosting interacting spin-less fermions.
We study numerically and analytically isolated interacting quantum systems that are taken out of equilibrium instantaneously (quenched). The probability of finding the initial state in time, the so-called fidelity, decays fastest for systems described by full random matrices, where simultaneous many-body interactions are implied. In the realm of realistic systems with two-body interactions, the dynamics is slower and depends on the interplay between the initial state and the Hamiltonian characterizing the system. The fastest fidelity decay in this case is Gaussian and can persist until saturation. A simple general picture, in which the fidelity plays a central role, is also achieved for the short-time dynamics of few-body observables. It holds for initial states that are eigenstates of the observables. We also discuss the need to reassess analytical expressions that were previously proposed to describe the evolution of the Shannon entropy. Our analyses are mainly developed for initial states that can be prepared in experiments with cold atoms in optical lattices.