This review is devoted to the problem of thermalization in a small isolated conglomerate of interacting constituents. A variety of physically important systems of intensive current interest belong to this category: complex atoms, molecules (including biological molecules), nuclei, small devices of condensed matter and quantum optics on nano-and micro-scale, cold atoms in optical lattices, ion traps. Physical implementations of quantum computers, where there are many interacting qubits, also fall into this group. Statistical regularities come into play through inter-particle interactions, which have two fundamental components: mean field, that along with external conditions, forms the regular component of the dynamics, and residual interactions responsible for the complex structure of the actual stationary states. At sufficiently high level density, the stationary states become exceedingly complicated superpositions of simple quasiparticle excitations. At this stage, regularities typical of quantum chaos emerge and bring in signatures of thermalization. We describe all the stages and the results of the processes leading to thermalization, using analytical and massive numerical examples for realistic atomic, nuclear, and spin systems, as well as for models with random parameters. The structure of stationary states, strength functions of simple configurations, and concepts of entropy and temperature in application to isolated mesoscopic systems are discussed in detail. We conclude with a schematic discussion of the time evolution of such systems to equilibrium.

By means of full exact diagonalization, we study level statistics and the structure of the eigenvectors of one-dimensional gapless bosonic and fermionic systems across the transition from integrability to quantum chaos. These systems are integrable in the presence of only nearest-neighbor terms, whereas the addition of next-nearest-neighbor hopping and interaction may lead to the onset of chaos. We show that the strength of the next-nearest-neighbor terms required to observe clear signatures of nonintegrability is inversely proportional to the system size. Interestingly, the transition to chaos is also seen to depend on particle statistics, with bosons responding first to the integrability breaking terms. In addition, we discuss the use of delocalization measures as main indicators for the crossover from integrability to chaos and the consequent viability of quantum thermalization in isolated systems.

We study one-dimensional lattices of interacting spins-1/2 and show that the effects of quenching the amplitude of a local magnetic field applied to a single site of the lattice can be comparable to the effects of a global perturbation applied instantaneously to the entire system. Both quenches take the system to the chaotic domain, the energy distribution of the initial states approaches a Breit-Wigner shape, the fidelity (Loschmidt echo) decays exponentially, and thermalization becomes viable.

We study the onset of chaos and statistical relaxation in two isolated dynamical quantum systems of interacting spins 1/2, one of which is integrable and the other chaotic. Our approach to identifying the emergence of chaos is based on the level of delocalization of the eigenstates with respect to the energy shell, the latter being determined by the interaction strength between particles or quasiparticles. We also discuss how the onset of chaos may be anticipated by a careful analysis of the Hamiltonian matrices, even before diagonalization. We find that despite differences between the two models, their relaxation processes following a quench are very similar and can be described analytically with a theory previously developed for systems with two-body random interactions. Our results imply that global features of statistical relaxation depend on the degree of spread of the eigenstates within the energy shell and may happen to both integrable and nonintegrable systems.

We study the transition to chaos and the emergence of statistical relaxation in isolated dynamical quantum systems of interacting particles. Our approach is based on the concept of delocalization of the eigenstates in the energy shell, controlled by the Gaussian form of the strength function. We show that, although the fluctuations of the energy levels in integrable and nonintegrable systems are different, the global properties of the eigenstates are quite similar, provided the interaction between particles exceeds some critical value. In this case, the statistical relaxation of the systems is comparable, irrespective of whether or not they are integrable. The numerical data for the quench dynamics manifest excellent agreement with analytical predictions of the theory developed for systems of two-body interactions with a completely random character.

We study how the proximity to an integrable point or to localization as one approaches the atomic limit, as well as the mixing of symmetries in the chaotic domain, may affect the onset of thermalization in finite one-dimensional systems. We consider systems of hard-core bosons at half-filling with nearest-neighbor hopping and interaction, and next-nearest-neighbor interaction. The latter breaks integrability and induces a ground-state superfluid to insulator transition. By full exact diagonalization, we study chaos indicators and few-body observables. We show that when different symmetry sectors are mixed, chaos indicators associated with the eigenvectors, contrary to those related to the eigenvalues, capture the onset of chaos. The results for the complexity of the eigenvectors and for the expectation values of few-body observables confirm the validity of the eigenstate thermalization hypothesis in the chaotic regime, and therefore the occurrence of thermalization. We also study the properties of the off-diagonal matrix elements of few-body observables in relation to the transition from integrability to chaos and from chaos to localization.

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. INTRODUCTIONDespite the ubiquity of many-body quantum systems out of equilibrium, they are much less understood than quantum systems in equilibrium. To advance our understanding and to construct a general picture, it is necessary to identify the elements that lead to similar dynamics. Determining how fast these systems evolve in time [1-5] is also essential for the development of algorithms for quantum optimal control [6]. In these two contexts, the unitary evolution of isolated many-body quantum systems is of particular interest due, in part, to the connection with current experiments in optical lattices [7][8][9][10][11][12][13][14][15][16][17]. The latter are quasi-isolated systems, where coherent evolutions can be studied for very long times.The evolution of an isolated system can be initiated by changing instantaneously the parameters of a certain initial Hamiltonian which is brought into a new final Hamiltonian. This abrupt perturbation is referred to as a quench. The system starts off in an eigenstate of the initial Hamiltonian. The fidelity (return probability) [18,19], which is defined as the overlap between the initial state and its evolved counterpart, is a way to characterize the system evolution. This quantity is related to the Loschmidt echo. It is also analogous to the characteristic function of the probability distribution of work [20][21][22][23] and is therefore likely to find applications in quantum thermodynamics, particularly in studies related with the quantification of the work done to take quantum systems out of equilibrium. The fidelity decays exponentially when the final Hamiltonian is chaotic [24][25][26][27][28][29][30][31][32]. In fact, this behavior is expected to hold even in integrable Hamiltonians provided the initial state be sufficiently delocalized in the energy eigenbasis [33][34][35].Here, we extend the results obtained in Ref. [36] and show that the fidelity can have a faster than exponential behavior. The fidelity ...

The isolated one-dimensional Heisenberg model with static random magnetic fields has become paradigmatic for the analysis of many-body localization. Here, we study the dynamics of this system initially prepared in a highly-excited nonstationary state. Our focus is on the probability for finding the initial state later in time, the so-called survival probability. Two distinct behaviors are identified before equilibration. At short times, the decay is very fast and equivalent to that of clean systems. It subsequently slows down and develops a power-law behavior with an exponent that coincides with the multifractal dimension of the eigenstates.