We consider the random-field Heisenberg model, a paradigmatic model for many-body localization (MBL), and add a Markovian dephasing bath coupled to the Anderson orbitals of the model's noninteracting limit. We map this system to a classical facilitated hopping model that is computationally tractable for large system sizes and investigate its dynamics. The classical model exhibits a robust crossover between an ergodic (thermal) phase and a frozen (localized) phase. The frozen phase is destabilized by thermal subregions (bubbles), which thermalize surrounding sites by providing a fluctuating interaction energy and so enable off-resonance particle transport. Investigating steady-state transport, we observe that the interplay between thermal and frozen bubbles leads to a clear transition between diffusive and subdiffusive regimes. This phenomenology both describes the MBL system coupled to a bath and provides a classical analog for the MBL transition in the corresponding quantum model, in that the classical model displays long local memory times. It also highlights the importance of the details of the bath coupling in studies of MBL systems coupled to thermal environments.
We consider the random-field Heisenberg model, a paradigmatic model for many-body localization (MBL), and add a Markovian dephasing bath coupled to the Anderson orbitals of the model's noninteracting limit. We map this system to a classical facilitated hopping model that is computationally tractable for large system sizes and investigate its dynamics. The classical model exhibits a robust crossover between an ergodic (thermal) phase and a frozen (localized) phase. The frozen phase is destabilized by thermal subregions (bubbles), which thermalize surrounding sites by providing a fluctuating interaction energy and so enable off-resonance particle transport. Investigating steady-state transport, we observe that the interplay between thermal and frozen bubbles leads to a clear transition between diffusive and subdiffusive regimes. This phenomenology both describes the MBL system coupled to a bath and provides a classical analog for the MBL transition in the corresponding quantum model, in that the classical model displays long local memory times. It also highlights the importance of the details of the bath coupling in studies of MBL systems coupled to thermal environments.
We study transport in a one-dimensional boundary-driven Anderson insulator (the XX spin chain with onsite disorder) with randomly positioned onsite dephasing, observing a transition from diffusive to subdiffusive spin transport below a critical density of sites with dephasing. This model is intended to mimic the passage of an excitation through (many-body) insulating regions or ergodic bubbles, therefore providing a toy model for the diffusion-subdiffusion transition observed in the disordered Heisenberg model by \ifmmode \check{Z}\else \v{Z}\fi{}nidari\ifmmode \check{c}\else \v{c}\fi{} et al. [Phys. Rev. Lett. 117, 040601 (2016)]. We also present the exact solution of a semiclassical model of conductors and insulators introduced by Agarwal et al. [Phys. Rev. Lett. 114, 160401 (2015)], which exhibits both diffusive and subdiffusive phases, and qualitatively reproduces the results of the quantum system. The critical properties of both models, when passing from diffusion to subdiffusion, are interpreted in terms of ``Griffiths effects.'' We show that the finite-size scaling comes from the interplay of three characteristic lengths: one associated with disorder (the localization length), one with dephasing, and the third with the percolation problem defining large, rare, insulating regions. We conjecture that the latter, which grows logarithmically with system size, may potentially be responsible for the fact that heavy-tailed resistance distributions typical of Griffiths effects have not been observed in subdiffusive interacting systems.
It is known that strong disorder in closed quantum systems leads to many-body localization (MBL), and that this quantum phase can be destroyed by coupling to an infinitely large Markovian environment. However, the stability of the MBL phase is less clear when the system and environment are of finite and comparable size. Here, we study the stability and eventual localization properties of a disordered Heisenberg spin chain coupled to a finite environment, and extensively explore the effects of environment disorder, geometry, initial state, and system-bath coupling strength, using the steady-state value of magnetization as a probe. Focusing on nonequilibrium dynamics and steady-state properties, our results indicate that within system sizes amenable to exact diagonalization, a strongly localized system interacting in a junction configuration retains remnant information on its initial state at long times despite coupling to a finite ergodic environment. In contrast, in a ladder configuration, strong dependencies on the initial state and coupling strength are observed, which can lead to either the loss or retention of information. Finally, we highlight and discuss discrepancies that can arise when similar methodologies are employed to infer localization or thermalization, revealing the need for careful interpretation.
We study transport in the boundary-driven XX spin chain with onsite disorder and randomly positioned onsite dephasing, observing a transition from diffusive to subdiffusive spin transport below a critical density of sites with dephasing. We then present an exactly solvable semiclassical model of conductors and insulators, which exhibits both diffusive and subdiffusive phases, and qualitatively reproduces the results of the quantum system. The subdiffusion in these models is a consequence of rare insulating regions, and they therefore exhibit, in a more controlled setup, the physics of which have been conjectured to cause the subdiffusive transport observed in interacting many-body localizable systems. For the quantum model we show that finite-size effects come from the interplay of three characteristic lengths: one associated with disorder (the localization length), one with dephasing, and the third with the percolation problem defining large, rare, insulating regions. The latter grows logarithmically with system size, and we conjecture that this may be the reason why the heavy-tailed distributions typical of Griffiths effects have not been observed in subdiffusive interacting systems.
We examine the standard model of many-body localization (MBL), i.e., the disordered chain of interacting spinless fermions, by representing it as the network in the many-body (MB) basis of noninteracting localized Anderson states. By studying eigenstates of the full Hamiltonian, for strong disorders we find that the dynamics is confined up to very long times to disconnected MB clusters in the Fock space. By keeping only resonant contributions and simplifying the quantum problem to rate equations (REs) for MB states, in analogy with percolation problems, the MBL transition is located via the universal cluster distribution and the emergence of the macroscopic cluster. On the ergodic side, our approximate RE approach to the relaxation processes captures well the diffusion transport, as found for the full quantum model. In a broad transient regime, we find an anomalous, i.e., subdiffusivelike, transport, emerging from weak links between MB states.
We study high temperature spin transport in a disordered Heisenberg chain in the ergodic regime. By employing a density matrix renormalization group technique for the study of the stationary states of the boundary-driven Lindblad equation we are able to study extremely large systems (400 spins). We find both a diffusive and a subdiffusive phase depending on the strength of the disorder and on the anisotropy parameter of the Heisenberg chain. Studying finite-size effects, we show numerically and theoretically that a very large crossover length exists that controls the passage of a clean-system dominated dynamics to one observed in the thermodynamic limit. Such a large length scale, being larger than the sizes studied before, explains previous conflicting results. We also predict spatial profiles of magnetization in steady states of generic nondiffusive systems.
Many-body localization(MBL) provides a mechanism by which isolated interacting quantum systems with disorder can avoid thermalization unlike ergodic systems satisfying the eigenstate thermalization hypothesis(ETH). Many-body localized systems retain signatures of their initial conditions at long times, whereas systems obeying ETH lose such information as they approach thermal equilibrium. Studying Non-equilibrium dynamics across ETH-MBL crossover is an important problem in condensed matter physics. Adiabatic control of parameters in interacting disordered systems provides a powerful framework to investigate MBL phases and their dynamical robustness. Using exact diagonalization and time-dependent numerical methods we study the effects of adiabatically ramped interactions in a disordered spin-1/2 XXZ chain, a paradigmatic model for exploring the many-body localization transition. By monitoring diagonal entropy density and entanglement entropy density growth across various ramp speeds, and system sizes. Our study incorporates Finite-size effects of spectral observables to probe the transition between ergodic and localized phases. The numerical results show that localized dynamical behavior remains largely intact under sufficiently slow ramp evolution, while increasing the driving rate promotes stronger excitation generation and larger entropy growth. This trend highlights the strong dependence of nonequilibrium adiabatic dynamics in disordered interacting quantum many-body systems.