No Arabic abstract
We focus on the many-body eigenstates across a localization-delocalization phase transition. To characterize the robustness of the eigenstates, we introduce the eigenstate overlaps $mathcal{O}$ with respect to the different boundary conditions. In the ergodic phase, the average of eigenstate overlaps $bar{mathcal{O}}$ is exponential decay with the increase of the system size indicating the fragility of its eigenstates, and this can be considered as an eigenstate-version butterfly effect of the chaotic systems. For localized systems, $bar{mathcal{O}}$ is almost size-independent showing the strong robustness of the eigenstates and the broken of eigenstate thermalization hypothesis. In addition, we find that the response of eigenstates to the change of boundary conditions in many-body localized systems is identified with the single-particle wave functions in Anderson localized systems. This indicates that the eigenstates of the many-body localized systems, as the many-body wave functions, may be independent of each other. We demonstrate that this is consistent with the existence of a large number of quasilocal integrals of motion in the many-body localized phase. Our results provide a new method to study localized and delocalized systems from the perspective of eigenstates.
We propose a new approach to probing ergodicity and its breakdown in quantum many-body systems based on their response to a local perturbation. We study the distribution of matrix elements of a local operator between the systems eigenstates, finding a qualitatively different behaviour in the many-body localized (MBL) and ergodic phases. To characterize how strongly a local perturbation modifies the eigenstates, we introduce the parameter ${cal G}(L)=langle ln (V_{nm}/delta) rangle$, which represents a disorder-averaged ratio of a typical matrix element of a local operator $V$ to the energy level spacing, $delta$; this parameter is reminiscent of the Thouless conductance in the single-particle localization. We show that the parameter ${cal G}(L)$ decreases with system size $L$ in the MBL phase, and grows in the ergodic phase. We surmise that the delocalization transition occurs when ${cal G}(L)$ is independent of system size, ${cal G}(L)={cal G}_csim 1$. We illustrate our approach by studying the many-body localization transition and resolving the many-body mobility edge in a disordered 1D XXZ spin-1/2 chain using exact diagonalization and time-evolving block decimation methods. Our criterion for the MBL transition gives insights into microscopic details of transition. Its direct physical consequences, in particular logarithmically slow transport at the transition, and extensive entanglement entropy of the eigenstates, are consistent with recent renormalization group predictions.
The level statistics in the transition between delocalized and localized {phases of} many body interacting systems is {considered}. We recall the joint probability distribution for eigenvalues resulting from the statistical mechanics for energy level dynamics as introduced by Pechukas and Yukawa. The resulting single parameter analytic distribution is probed numerically {via Monte Carlo method}. The resulting higher order spacing ratios are compared with data coming from different {quantum many body systems}. It is found that this Pechukas-Yukawa distribution compares favorably with {$beta$--Gaussian ensemble -- a single parameter model of level statistics proposed recently in the context of disordered many-body systems.} {Moreover, the Pechukas-Yukawa distribution is also} only slightly inferior to the two-parameter $beta$-h ansatz shown {earlier} to reproduce {level statistics of} physical systems remarkably well.
We prove that all eigenstates of many-body localized symmetry protected topological systems with time reversal symmetry have four-fold degenerate entanglement spectra in the thermodynamic limit. To that end, we employ unitary quantum circuits where the number of sites the gates act on grows linearly with the system size. We find that the corresponding matrix product operator representation has similar local symmetries as matrix product ground states of symmetry protected topological phases. Those local symmetries give rise to a $mathbb{Z}_2$ topological index, which is robust against arbitrary perturbations so long as they do not break time reversal symmetry or drive the system out of the fully many-body localized phase.
The many-body localization transition (MBLT) between ergodic and many-body localized phase in disordered interacting systems is a subject of much recent interest. Statistics of eigenenergies is known to be a powerful probe of crossovers between ergodic and integrable systems in simpler examples of quantum chaos. We consider the evolution of the spectral statistics across the MBLT, starting with mapping to a Brownian motion process that analytically relates the spectral properties to the statistics of matrix elements. We demonstrate that the flow from Wigner-Dyson to Poisson statistics is a two-stage process. First, fractal enhancement of matrix elements upon approaching the MBLT from the metallic side produces an effective power-law interaction between energy levels, and leads to a plasma model for level statistics. At the second stage, the gas of eigenvalues has local interaction and level statistics belongs to a semi-Poisson universality class. We verify our findings numerically on the XXZ spin chain. We provide a microscopic understanding of the level statistics across the MBLT and discuss implications for the transition that are strong constraints on possible theories.
Numerical studies of amorphous silicon in harmonic approximation show that the highest 3.5% of vibrational normal modes are localized. As vibrational frequency increases through the boundary separating localized from delocalized modes, near omega_c=70meV, (the mobility edge) there is a localization-delocalization (LD) transition, similar to a second-order thermodynamic phase transition. By a numerical study on a system with 4096 atoms, we are able to see exponential decay lengths of exact vibrational eigenstates, and test whether or not these diverge at omega_c. Results are consistent with a localization length xi which diverges above omega_c as (omega-omega_c)^{-p} where the exponent is p = 1.3 +/- 0.5. Below the mobility edge we find no evidence for a diverging correlation length. Such an asymmetry would contradict scaling ideas, and we suppose it is a finite-size artifact. If the scaling regime is narrower than our 1 meV resolution, then it cannot be seen directly on our finite system.