No Arabic abstract
We numerically study the level statistics of the Gaussian $beta$ ensemble. These statistics generalize Wigner-Dyson level statistics from the discrete set of Dyson indices $beta = 1,2,4$ to the continuous range $0 < beta < infty$. The Gaussian $beta$ ensemble covers Poissonian level statistics for $beta to 0$, and provides a smooth interpolation between Poissonian and Wigner-Dyson level statistics. We establish the physical relevance of the level statistics of the Gaussian $beta$ ensemble by showing near-perfect agreement with the level statistics of a paradigmatic model in studies on many-body localization over the entire crossover range from the thermal to the many-body localized phase. In addition, we show similar agreement for a related Hamiltonian with broken time-reversal symmetry.
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.
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.
The intriguing phenomenon of many-body localization (MBL) has attracted significant interest recently, but a complete characterization is still lacking. In this work, we introduce the total correlations, a concept from quantum information theory capturing multi-partite correlations, to the study of this phenomenon. We demonstrate that the total correlations of the diagonal ensemble provides a meaningful diagnostic tool to pin-down, probe, and better understand the MBL transition and ergodicity breaking in quantum systems. In particular, we show that the total correlations has sub-linear dependence on the system size in delocalized, ergodic phases, whereas we find that it scales extensively in the localized phase developing a pronounced peak at the transition. We exemplify the power of our approach by means of an exact diagonalization study of a Heisenberg spin chain in a disordered field.
In [Van Beeumen, et. al, HPC Asia 2020, https://www.doi.org/10.1145/3368474.3368497] a scalable and matrix-free eigensolver was proposed for studying the many-body localization (MBL) transition of two-level quantum spin chain models with nearest-neighbor $XX+YY$ interactions plus $Z$ terms. This type of problem is computationally challenging because the vector space dimension grows exponentially with the physical system size, and averaging over different configurations of the random disorder is needed to obtain relevant statistical behavior. For each eigenvalue problem, eigenvalues from different regions of the spectrum and their corresponding eigenvectors need to be computed. Traditionally, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large number of spins $L$, e.g., one needs thousands of compute nodes on modern high performance computing infrastructures to go beyond $L = 24$. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication imbalance. We present a few strategies to reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver. To optimize the communication performance, we leverage the consistent space runtime, CSPACER, and show its efficiency in accelerating the MBL irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.
Recent developments in matrix-product-state (MPS) investigations of many-body localization (MBL) are reviewed, with a discussion of benefits and limitations of the method. This approach allows one to explore the physics around the MBL transition in systems much larger than those accessible to exact diagonalization. System sizes and length scales that can be controllably accessed by the MPS approach are comparable to those studied in state-of-the-art experiments. Results for 1D, quasi-1D, and 2D random systems, as well as 1D quasi-periodic systems are presented. On time scales explored (up to $t approx 300$ in units set by the hopping amplitude), a slow, subdiffusive transport in a rather broad disorder range on the ergodic side of the MBL transition is found. For 1D random spin chains, which serve as a standard model of the MBL transition, the MPS study demonstrates a substantial drift of the critical point $W_c(L)$ with the system size $L$: while for $L approx 20$ we find $W_c approx 4$, as also given by exact diagonalization, the MPS results for $L = 50$--100 provide evidence that the critical disorder saturates, in the large-$L$ limit, at $W_c approx 5.5$. For quasi-periodic systems, these finite-size effects are much weaker, which suggests that they can be largely attributed to rare events. For quasi-1D ($dtimes L$, with $d ll L$) and 2D ($Ltimes L$) random systems, the MPS data demonstrate an unbounded growth of $W_c$ in the limit of large $d$ and $L$, in agreement with analytical predictions based on the rare-event avalanche theory.