We study the statistical properties of the off-diagonal matrix elements of observables in the energy eigenstates of integrable quantum systems. They have been found to be dense in the spin-1/2 XXZ chain, while they are sparse in noninteracting systems. We focus on the quasimomentum occupation of hard-core bosons in one dimension and show that the distributions of the off-diagonal matrix elements are well described by generalized Gamma distributions, in both the presence and absence of translational invariance but not in the presence of localization. We also show that the results obtained for the off-diagonal matrix elements of observables in the spin-1/2 XXZ model are well described by a generalized Gamma distribution.
We study the statistical properties of the off-diagonal matrix elements of observables in the energy eigenstates of integrable quantum systems. They have been found to be dense in the spin-1/2 XXZ chain, while they are sparse in noninteracting systems. We focus on the quasimomentum occupation of hard-core bosons in one dimension and show that the distributions of the off-diagonal matrix elements are well described by generalized Gamma distributions, in both the presence and absence of translational invariance but not in the presence of localization. We also show that the results obtained for the off-diagonal matrix elements of observables in the spin-1/2 XXZ model are well described by a generalized Gamma distribution.
We study the off-diagonal matrix elements of observables that break the translational symmetry of a spin-chain Hamiltonian, and as such connect energy eigenstates from different total quasimomentum sectors. We consider quantum-chaotic and interacting integrable points of the Hamiltonian, and focus on average energies at the center of the spectrum. In the quantum-chaotic model, we find that there is eigenstate thermalization; specifically, the matrix elements are Gaussian distributed with a variance that is a smooth function of ω=E_{α}-E_{β} (E_{α} are the eigenenergies) and scales as 1/D (D is the Hilbert space dimension). In the interacting integrable model, we find that the matrix elements exhibit a skewed log-normal-like distribution and have a variance that is also a smooth function of ω that scales as 1/D. We study in detail the low-frequency behavior of the variance of the matrix elements to unveil the regimes in which it exhibits diffusive or ballistic scaling. We show that in the quantum-chaotic model the behavior of the variance is qualitatively similar for matrix elements that connect eigenstates from the same versus different quasimomentum sectors. We also show that this is not the case in the interacting integrable model for observables whose translationally invariant counterpart does not break integrability if added as a perturbation to the Hamiltonian.
We study the bipartite von Neumann entanglement entropy and matrix elements of local operators in the eigenstates of an interacting integrable Hamiltonian (the paradigmatic spin-1/2 XXZ chain), and we contrast their behavior with that of quantum chaotic systems. We find that the leading term of the average (over all eigenstates in the zero magnetization sector) eigenstate entanglement entropy has a volume-law coefficient that is smaller than the universal (maximal entanglement) one in quantum chaotic systems. This establishes the entanglement entropy as a powerful measure to distinguish integrable models from generic ones. Remarkably, our numerical results suggest that the volume-law coefficient of the average entanglement entropy of eigenstates of the spin-1/2 XXZ Hamiltonian is very close to, or the same as, the one for translationally invariant quadratic fermionic models. We also study matrix elements of local operators in the eigenstates of the spin-1/2 XXZ Hamiltonian at the center of the spectrum. For the diagonal matrix elements, we show evidence that the support does not vanish with increasing system size, while the average eigenstate-to-eigenstate fluctuations vanish in a power-law fashion. For the off-diagonal matrix elements, we show that they follow a distribution that is close to (but not quite) log-normal, and that their variance is a well-defined function of ω=E_{α}-E_{β} ({E_{α}} are the eigenenergies) proportional to 1/D, where D is the Hilbert space dimension.
The thermalization phenomenon and many-body quantum statistical properties are studied on the example of several observables in isolated spin-chain systems, both integrable and generic nonintegrable. While diagonal matrix elements for nonintegrable models comply with the eigenstate thermalization hypothesis, the integrable systems show evident deviations and similarity to properties of noninteracting many-fermion models. The finite-size scaling reveals that the crossover between the two regimes is given by a scale closely related to the scattering length. Low-frequency off-diagonal matrix elements related to dc transport quantities also follow in a generic system a behavior analogous to the eigenstate thermalization hypothesis, however, unrelated to one of the diagonal matrix elements.
We investigate off-diagonal matrix elements of local operators in integrable spin chains, focusing on the isotropic spin- 1 / 2 Heisenberg chain ( X X X chain). We employ state-of-the-art algebraic Bethe ansatz results, which allow us to efficiently compute matrix elements of operators with support up to two sites between generic energy eigenstates. We consider both matrix elements between eigenstates that are in the same thermodynamic macrostate, as well as eigenstates that belong to different macrostates. In the former case, focusing on thermal states we numerically show that matrix elements are compatible with the exponential decay as exp ( − L | M i j O | ) . The probability distribution functions of M i j O depend on the observable and on the macrostate, and are well described by Gumbel distributions. On the other hand, matrix elements between eigenstates in different macrostates decay faster as exp ( − | M ′ i j O | L 2 ) , with M ′ i j O , again, compatible with a Gumbel distribution.
In the time evolution of isolated quantum systems out of equilibrium, local observables generally relax to a long-time asymptotic value, governed by the expectation values (diagonal matrix elements) of the corresponding operator in the eigenstates of the system. The temporal fluctuations around this value, response to further perturbations, and the relaxation toward this asymptotic value are all determined by the off-diagonal matrix elements. Motivated by this nonequilibrium role, we present generic statistical properties of off-diagonal matrix elements of local observables in two families of interacting many-body systems with local interactions. Since integrability (or lack thereof) is an important ingredient in the relaxation process, we analyze models that can be continuously tuned to integrability. We show that, for generic nonintegrable systems, the distribution of off-diagonal matrix elements is a Gaussian centered at zero. As one approaches integrability, the peak around zero becomes sharper, so the distribution is approximately a combination of two Gaussians. We characterize the proximity to integrability through the deviation of this distribution from a Gaussian shape. We also determine the scaling dependence on system size of the average magnitude of off-diagonal matrix elements.
We study the matrix elements of local operators in the eigenstates of the integrable XXZ chain and of the quantum-chaotic model obtained by locally perturbing the XXZ chain with a magnetic impurity. We show that, at frequencies that are polynomially small in the system size, the behavior of the variances of the off-diagonal matrix elements can be starkly different depending on the operator. In the integrable model we find that, as the frequency $\ensuremath{\omega}\ensuremath{\rightarrow}0$, the variances are either nonvanishing (generic behavior) or vanishing (for a special class of operators). In the quantum-chaotic model, on the other hand, we find the variances to be nonvanishing as $\ensuremath{\omega}\ensuremath{\rightarrow}0$ and to indicate diffusive dynamics. We highlight which properties of the matrix elements of local operators are different between the integrable and quantum-chaotic models independently of the specific operator selected.