We introduce and study a toy model for anomalous transport and Griffiths effects in one-dimensional (1D) quantum disordered isolated systems near the many-body localization (MBL) transition. The model is constituted by a collection of 1D tight-binding chains with on-site random energies, locally coupled to a weak perturbation, corresponding to a random matrix which belongs to the Gaussian orthogonal ensemble (GOE) and mimics the effect of thermal inclusions due to delocalizing interactions by providing a local broadening of the Poisson spectrum. While in the absence of such a coupling the model is localized as expected for the one-dimensional Anderson model, increasing the coupling with the GOE perturbation we find a transition to a conducting phase driven by the proliferation of quantum avalanches, which does not fit the standard paradigm of the Anderson transition. In particular an intermediate Griffiths region emerges, where exponentially distributed insulating segments coexist with a few rare resonances. Typical correlations decay exponentially fast, while average correlations decay as stretched exponentials and diverge with the length of the chain, indicating that the conducting inclusions have a fractal structure and that the localization length is broadly distributed at the critical point. This behavior is consistent with a Kosterlitz-Thouless-like criticality of the MBL transition. Transport and relaxation are dominated by rare resonances and rare strong insulating regions, and show anomalous behaviors strikingly similar to those observed in recent simulations and experiments in the bad metal delocalized phase preceding MBL. In particular, we find subdiffusive transport and power-law decay of the return probability at large times, with exponents that gradually change as one moves across the intermediate region. Concomitantly, the ac conductivity vanishes near zero frequency with an anomalous power law.
The study of spectrum statistics, such as the consecutive-gap ratio distribution, has revealed many interesting properties of many-body complex systems. Here we propose a two-parameter surmise expression for such distribution to describe the crossover between the Gaussian orthogonal ensemble (GOE) and Poisson statistics. This crossover is observed in the isotropic Heisenberg spin-
1/2 chain with disordered local field, exhibiting the Many-Body Localization (MBL) transition. Inspired by the analysis of stability in dynamical systems, this crossover is presented as a flow pattern in the parameter space, with the Poisson statistics being the fixed point of the system, which represents the MBL phase. We also analyze an isotropic Heisenberg spin-
1/2 chain with disordered local exchange coupling and a zero magnetic field. In this case, the system never achieves the MBL phase because of the spin rotation symmetry. This case is more sensitive to finite-size effects than the previous one, and thus the flow pattern resembles a two-dimensional random walk close to its fixed point. We propose a system of linearized stochastic differential equations to estimate this fixed point. We study the continuous-state Markov process that governs the probability of finding the system close to this fixed point as the disorder strength increases. In addition, we discuss the conditions under which the stationary probability distribution is given by a bivariate normal distribution.
Recent studies on disorder-induced many-body localization (MBL) in non-Hermitian quantum systems have attracted great interest. However, the non-Hermitian disorder-free MBL still needs to be clarified. We consider a one-dimensional interacting Stark model with nonreciprocal hoppings having time-reversal symmetry, the properties of which are boundary dependent. Under periodic boundary conditions (PBCs), such a model exhibits three types of phase transitions: the real-complex transition of eigenenergies, the topological phase transition, and the non-Hermitian Stark MBL transition. The real-complex and topological phase transitions occur at the same point in the thermodynamic limit but do not coincide with the non-Hermitian Stark MBL transition, which is quite different from the non-Hermitian disordered cases. By the level statistics, the system transitions from the Ginibre ensemble (GE) to the Gaussian orthogonal ensemble (GOE) to the Possion ensemble with the increase of the linear tilt potential's strength. The real-complex transition of the eigenvalues is accompanied by the GE-to-GOE transition in the ergodic regime. Moreover, the second transition of the level statistics corresponds to the occurrence of non-Hermitian Stark MBL. We demonstrate that the non-Hermitian Stark MBL is robust and shares many similarities with disorder-induced MBL, which several existing characteristic quantities of the spectral statistics and eigenstate properties can confirm. The dynamical evolutions of the entanglement entropy and the density imbalance can distinguish the real-complex and Stark MBL transitions. Finally, we find that our system under open boundary conditions lacks a real-complex transition, and the transition of non-Hermitian Stark MBL is the same as that under PBCs.