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 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 explore the role of the initial state on the onset of thermalization in isolated quantum many-body systems after a quench. The initial state is an eigenstate of an initial Hamiltonian ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{I}$ and it evolves according to a different final Hamiltonian ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$. If the initial state has a chaotic structure with respect to ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$, i.e., if it fills the energy shell ergodically, thermalization is certain to occur. This happens when ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{I}$ is a full random matrix, because its states projected onto ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$ are fully delocalized. The results for the observables then agree with those obtained with thermal states at infinite temperature. However, finite real systems with few-body interactions, as the ones considered here, are deprived of fully extended eigenstates, even when described by a nonintegrable Hamiltonian. We examine how the initial state delocalizes as it gets closer to the middle of the spectrum of ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}}_{F}$, causing the observables to approach thermal averages, be the models integrable or chaotic. Our numerical studies are based on initial states with energies that cover the entire lower half of the spectrum of one-dimensional Heisenberg spin-1/2 systems.
Using the ergodicity principle for the expectation values of several types of observables, we investigate the thermalization process in isolated fermionic systems.These are described by the two-body random ensemble, which is a paradigmatic model to study quantum chaos and especially the dynamical transition from integrability to chaos.By means of exact diagonalizations we analyze the relevance of the eigenstate thermalization hypothesis as well as the influence of other factors, such as the energy and structure of the initial state, or the dimension of the Hilbert space.We also obtain analytical expressions linking the degree of thermalization for a given observable with the so-called number of principal components for transition strengths originating at a given energy, with the dimensions of the whole Hilbert space and microcanonical energy shell, and with the correlations generated by the observable.As the strength of the residual interaction is increased, an order-to-chaos transition takes place, and we show that the onset of Wigner spectral fluctuations, which is the standard signature of chaos, is not sufficient to guarantee thermalization in finite systems.When all the signatures of chaos are fulfilled, including the quasicomplete delocalization of eigenfunctions, the eigenstate thermalization hypothesis is the mechanism responsible for the thermalization of certain types of observables, such as (linear combinations of) occupancies and strength function operators.Our results also suggest that fully chaotic systems will thermalize relative to most observables in the thermodynamic limit.
Using the ergodicity principle for the expectation values of several types of observables, we investigate the thermalization process in isolated fermionic systems. These are described by the two-body random ensemble, which is a paradigmatic model to study quantum chaos and especially the dynamical transition from integrability to chaos. By means of exact diagonalizations we analyze the relevance of the eigenstate thermalization hypothesis as well as the influence of other factors, such as the energy and structure of the initial state, or the dimension of the Hilbert space. We also obtain analytical expressions linking the degree of thermalization for a given observable with the so-called number of principal components for transition strengths originating at a given energy, with the dimensions of the whole Hilbert space and microcanonical energy shell, and with the correlations generated by the observable. As the strength of the residual interaction is increased, an order-to-chaos transition takes place, and we show that the onset of Wigner spectral fluctuations, which is the standard signature of chaos, is not sufficient to guarantee thermalization in finite systems. When all the signatures of chaos are fulfilled, including the quasicomplete delocalization of eigenfunctions, the eigenstate thermalization hypothesis is the mechanism responsible for the thermalization of certain types of observables, such as (linear combinations of) occupancies and strength function operators. Our results also suggest that fully chaotic systems will thermalize relative to most observables in the thermodynamic limit.
Numerically, we study the time fluctuations of few-body observables after relaxation in isolated dynamical quantum systems of interacting particles. Our results suggest that they decay exponentially with system size in both regimes, integrable and chaotic. The integrable systems considered are solvable with the Bethe ansatz and have a highly nondegenerate spectrum. This is in contrast with integrable Hamiltonians mappable to noninteracting ones. We show that the coefficient of the exponential decay depends on the level of delocalization of the initial state with respect to the energy shell.
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.