Non-Hermitian many-body localization (NH MBL) has emerged as a possible scenario for stable localization in open systems, as suggested by spectral indicators identifying a putative transition for finite system sizes. In this work, we shift the focus to dynamical probes, specifically the steady-state spin current, to investigate transport properties in a disordered, non-Hermitian XXZ spin chain. Through exact diagonalization for small systems and tensor-network methods for larger chains, we demonstrate that the steady-state current remains finite and decays exponentially with disorder strength, showing no evidence of a transition up to disorder values far beyond the previously claimed critical point. Our results reveal a stark discrepancy between spectral indicators, which suggest localization, and transport behavior, which indicates delocalization. This highlights the importance of dynamical observables in characterizing NH MBL and suggests that traditional spectral measures may not fully capture the physics of non-Hermitian systems. Additionally, we observe a noncommutativity of limits in system size and time, further complicating the interpretation of finite-size studies. These findings challenge the existence of NH MBL in the studied model and underscore the need for alternative approaches to understanding localization in non-Hermitian settings.
Non-Hermitian many-body localization (NH MBL) has emerged as a possible scenario for stable localization in open systems, as suggested by spectral indicators identifying a putative transition for finite system sizes. In this work, we shift the focus to dynamical probes, specifically the steady-state spin current, to investigate transport properties in a disordered, non-Hermitian XXZ spin chain. Through exact diagonalization for small systems and tensor-network methods for larger chains, we demonstrate that the steady-state current remains finite and decays exponentially with disorder strength, showing no evidence of a transition up to disorder values far beyond the previously claimed critical point. Our results reveal a stark discrepancy between spectral indicators, which suggest localization, and transport behavior, which indicates delocalization. This highlights the importance of dynamical observables in characterizing NH MBL and suggests that traditional spectral measures may not fully capture the physics of non-Hermitian systems. Additionally, we observe a noncommutativity of limits in system size and time, further complicating the interpretation of finite-size studies. These findings challenge the existence of NH MBL in the studied model and underscore the need for alternative approaches to understanding localization in non-Hermitian settings.
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.
Abstract Statistical mechanics provides a framework for describing the physics of large, complex many-body systems using only a few macroscopic parameters to determine the state of the system. For isolated quantum many-body systems, such a description is achieved via the eigenstate thermalization hypothesis (ETH), which links thermalization, ergodicity and quantum chaotic behavior. However, tendency towards thermalization is not observed at finite system sizes and evolution times in a robust many-body localization (MBL) regime found numerically and experimentally in the dynamics of interacting many-body systems at strong disorder. Although the phenomenology of the MBL regime is well-established, the central question remains unanswered: under what conditions does the MBL regime give rise to an MBL phase , in which the thermalization does not occur even in the asymptotic limit of infinite system size and evolution time? This review focuses on recent numerical investigations aiming to clarify the status of the MBL phase, and it establishes the critical open questions about the dynamics of disordered many-body systems. The last decades of research have brought an unprecedented new variety of tools and indicators to study the breakdown of ergodicity, ranging from spectral and wave function measures, matrix elements of observables, through quantities probing unitary quantum dynamics, to transport and quantum information measures. We give a comprehensive overview of these approaches and attempt to provide a unified understanding of their main features. We emphasize general trends towards ergodicity with increasing length and time scales, which exclude naive single-parameter scaling hypothesis, necessitate the use of more refined scaling procedures, and prevent unambiguous extrapolations of numerical results to the asymptotic limit. Providing a concise description of numerical methods for studying ETH and MBL, we explore various approaches to tackle the question of the MBL phase. Persistent finite size drifts towards ergodicity consistently emerge in quantities derived from eigenvalues and eigenvectors of disordered many-body systems. The drifts are related to continuous inching towards ergodicity and non-vanishing transport observed in the dynamics of many-body systems, even at strong disorder. These phenomena impede the understanding of microscopic processes at the ETH-MBL crossover. Nevertheless, the abrupt slowdown of dynamics with increasing disorder strength provides premises suggesting the proximity of the MBL phase. This review concludes that the questions about thermalization and its failure in disordered many-body systems remain a captivating area open for further explorations.
Localization-delocalization transition and PT-symmetry breaking via second-neighbor coupling in a Heisenberg spin chain
The exploration of non-Hermitian quantum systems, particularly those governed by parity-time ([Formula: see text])-symmetry, has revealed a rich landscape of unconventional phase transitions and localization phenomena. Spin chains with quasiperiodic potentials and competing Heisenberg interactions provide a versatile framework for probing these effects. In this study, we investigate the dynamical and spectral signatures of [Formula: see text]-symmetry breaking and localization-delocalization transitions in spin chains subjected to a non-Hermitian Aubry-André potential and tunable Heisenberg couplings across nearest ([Formula: see text]) and next-nearest neighbors ([Formula: see text]). Our analysis reveals that increasing [Formula: see text] monotonically shifts the critical point [Formula: see text], while [Formula: see text] induces a non-linear response, generating a distinct minimum in [Formula: see text]. Notably, [Formula: see text] leads to the emergence of partially imaginary eigenvalues beyond [Formula: see text], decoupling the onset of [Formula: see text]-symmetry breaking from the localization transition. Using inverse and normalized participation ratios, we identify three distinct phases: extended, localized, and an intermediate hybrid regime characterized by a mobility edge. To capture the dynamical fingerprints of these phases, we employ time-dependent density distribution, long-time survival probability [Formula: see text], analysis. Extended phases exhibit ballistic spreading and algebraic decay, while localized regimes show spatial confinement and exponential decay. These observables serve as robust indicators of non-Hermitian phase transitions, offering insights beyond static spectral measures. Overall, our findings underscore the critical role of dynamical metrics in characterizing phase structure in non-Hermitian spin systems. As interest in non-equilibrium quantum dynamics continues to grow, such models provide a compelling platform for understanding the interplay between [Formula: see text]-symmetry, localization, and quantum coherence in complex many-body systems.
The disordered XXZ model is a prototype model of the many-body localization (MBL) transition. Despite numerous studies of this model, the available numerical evidence of multifractality of its eigenstates is not very conclusive due to severe finite size effects. Moreover it is not clear if similarly to the case of single-particle physics, multifractal properties of the many-body eigenstates are related to anomalous transport, which is observed in this model. In this work, using a state-of-the-art, massively parallel, numerically exact method, we study systems of up to 24 spins and show that a large fraction of the delocalized phase flows towards ergodicity in the thermodynamic limit, while a region immediately preceding the MBL transition appears to be multifractal in this limit. We discuss the implication of our finding on the mechanism of subdiffusive transport.
The disordered XXZ model is a prototype model of the many-body localization transition (MBL). Despite numerous studies of this model, the available numerical evidence of multifractality of its eigenstates is not very conclusive due severe finite size effects. Moreover it is not clear if similarly to the case of single-particle physics, multifractal properties of the many-body eigenstates are related to anomalous transport, which is observed in this model. In this work, using a state-of-the-art, massively parallel, numerically exact method, we study systems of up to 24 spins and show that a large fraction of the delocalized phase flows towards ergodicity in the thermodynamic limit, while a region immediately preceding the MBL transition appears to be multifractal in this limit. We discuss the implication of our finding on the mechanism of subdiffusive transport.
We study time dynamics of 1D disordered Heisenberg spin-1/2 chain focusing on a regime of large system sizes and a long time evolution. This regime is relevant for observation of many-body localization (MBL), a phenomenon that is expected to freeze the dynamics of the system and prevent it from reaching thermal equilibrium. Performing extensive numerical simulations of the imbalance, a quantity often employed in the experimental studies of MBL, we show that the regime of a slow power-law decay of imbalance persists to disorder strengths exceeding by at least a factor of 2 the current estimates of the critical disorder strength for MBL. Even though we investigate time evolution up to few thousands tunneling times, we observe no signs of the saturation of imbalance that would suggest freezing of system dynamics and provide a smoking gun evidence of MBL. We demonstrate that the situation is qualitatively different when the disorder is replaced by a quasiperiodic potential. In this case, we observe an emergence of a pattern of oscillations of the imbalance that is stable with respect to changes in the system size. This suggests that the dynamics of quasiperiodic systems remain fully local at the longest time scales we reach provided that the quasiperiodic potential is sufficiently strong. Our study identifies challenges in an unequivocal experimental observation of the phenomenon of MBL.