Many-body localisation is studied in a disordered quantum spin-1/2 chain with long-ranged power-law interactions, and distinct power-law exponents for interactions between longitudinal and transverse spin components. Using a self-consistent mean-field theory centring on the local propagator in Fock space and its associated self-energy, a localisation phase diagram is obtained as a function of the power-law exponents and the disorder strength of the random fields acting on longitudinal spin-components. Analytical results are corroborated using the well-studied and complementary numerical diagnostics of level statistics, entanglement entropy, and participation entropy, obtained via exact diagonalisation. We find that increasing the range of interactions between transverse spin components hinders localisation and enhances the critical disorder strength. In marked contrast, increasing the interaction range between longitudinal spin components is found to enhance localisation and lower the critical disorder.
Many-body localisation is studied in a disordered quantum spin-1/2 chain with long-ranged power-law interactions, and distinct power-law exponents for interactions between longitudinal and transverse spin components. Using a self-consistent mean-field theory centring on the local propagator in Fock space and its associated self-energy, a localisation phase diagram is obtained as a function of the power-law exponents and the disorder strength of the random fields acting on longitudinal spin-components. Analytical results are corroborated using the well-studied and complementary numerical diagnostics of level statistics, entanglement entropy, and participation entropy, obtained via exact diagonalisation. We find that increasing the range of interactions between transverse spin components hinders localisation and enhances the critical disorder strength. In marked contrast, increasing the interaction range between longitudinal spin components is found to enhance localisation and lower the critical disorder.
Many-body localisation is studied in a disordered quantum spin-1/2 chain with long-ranged power-law interactions, and distinct power-law exponents for interactions between longitudinal and transverse spin components. Using a self-consistent mean-field theory centring on the local propagator in Fock space and its associated self-energy, a localisation phase diagram is obtained as a function of the power-law exponents and the disorder strength of the random fields acting on longitudinal spin-components. Analytical results are corroborated using the well-studied and complementary numerical diagnostics of level statistics, entanglement entropy, and participation entropy, obtained via exact diagonalisation. We find that increasing the range of interactions between transverse spin components hinders localisation and enhances the critical disorder strength. In marked contrast, increasing the interaction range between longitudinal spin components is found to enhance localisation and lower the critical disorder.
We investigate the transition from the many-body localized phase to the ergodic thermalized phase at an infinite temperature in an XY spin chain with L spins, which experiences power-law decaying interactions in the form of ${V}_{ij}\ensuremath{\propto}1/{\left|i\ensuremath{-}j\right|}^{\ensuremath{\alpha}}$ (i,j=1,...,L) and a random transverse field. By performing large-scale exact diagonalization for the chain size up to L=18, we systematically analyze the energy gap statistics, the half-chain entanglement entropy, and the uncertainty of the entanglement entropy of the system at different interaction exponents $\ensuremath{\alpha}$. The finite-size critical scaling allows us to determine the critical disorder strength Wc and the critical exponent $\ensuremath{\nu}$ at the many-body-localization phase transition, as a function of the interaction exponent $\ensuremath{\alpha}$ in the limit $L\ensuremath{\rightarrow}\ensuremath{\infty}$. We find that both Wc and $\ensuremath{\nu}$ diverge when $\ensuremath{\alpha}$ decreases to a critical power of ${\ensuremath{\alpha}}_{c}\ensuremath{\simeq}1.16\ifmmode\pm\else\textpm\fi{}0.17$, indicating the absence of many-body localization for $\ensuremath{\alpha}<{\ensuremath{\alpha}}_{c}$. Our result differs from the analytical prediction of ${\ensuremath{\alpha}}_{c}=3/2$ found in a previous work based on a scaling argument [A. L. Burin, Phys. Rev. B 92, 104428 (2015)], but is close to the recent numerical result ${\ensuremath{\alpha}}_{c}\ensuremath{\approx}1$ obtained from quantum dynamics simulation [A. Safavi-Naini, M. L. Wall, O. L. Acevedo, A. M. Rey, and R. M. Nandkishore, Phys. Rev. A 99, 033610 (2019)].
We theoretically investigate the many-body localization phase transition in a one-dimensional Ising spin chain with random long-range spin-spin interactions, ${V}_{ij}\ensuremath{\propto}{\left|i\ensuremath{-}j\right|}^{\ensuremath{-}\ensuremath{\alpha}}$, where the exponent of the interaction range $\ensuremath{\alpha}$ can be tuned from zero to infinitely large. By using exact diagonalization, we calculate the half-chain entanglement entropy and the energy spectral statistics and use them to characterize the phase transition towards the many-body localization phase at infinite temperature and at sufficiently large disorder strength. We perform finite-size scaling to extract the critical disorder strength and the critical exponent of the divergent localization length. With increasing $\ensuremath{\alpha}$, the critical exponent experiences a sharp increase at about ${\ensuremath{\alpha}}_{c}\ensuremath{\simeq}1.2$ and then gradually decreases to a value found earlier in a disordered short-ranged interacting spin chain. For $\ensuremath{\alpha}<{\ensuremath{\alpha}}_{c}$, we find that the system is mostly localized and the increase in the disorder strength may drive a transition between two many-body localized phases. In contrast, for $\ensuremath{\alpha}>{\ensuremath{\alpha}}_{c}$, the transition is from a thermalized phase to the many-body localization phase. Our predictions could be experimentally tested with an ion-trap quantum emulator with programmable random long-range interactions, or with randomly distributed Rydberg atoms or polar molecules in lattices.
We investigate many-body localization in the chain of interacting spins with a transverse power-law interaction, ${J}_{0}/{r}^{\ensuremath{\alpha}}$, and random on-site potentials, ${\ensuremath{\phi}}_{i}\ensuremath{\in}\left(\ensuremath{-}W/2,W/2\right)$, in the long-range limit, $\ensuremath{\alpha}<3/2$, which has been recently examined experimentally on trapped ions. The many-body localization threshold is characterized by the critical disordering, Wc, which separates localized (W>Wc) and chaotic (W<Wc) phases. Using the analysis of the instability of localized states with respect to resonant interactions complemented by numerical finite-size scaling, we show that the critical disordering scales with the number of spins, N, as ${W}_{c}\ensuremath{\approx}[1.37{J}_{0}/(4/3\ensuremath{-}\ensuremath{\alpha})]{N}^{4/3\ensuremath{-}\ensuremath{\alpha}}lnN$ for $0<\ensuremath{\alpha}\ensuremath{\le}1$, and as ${W}_{c}\ensuremath{\approx}[{J}_{0}/(1\ensuremath{-}2\ensuremath{\alpha}/3)]{N}^{1\ensuremath{-}2\ensuremath{\alpha}/3}{ln}^{2/3}N$ for $1<\ensuremath{\alpha}<3/2$, while the transition width scales as ${\ensuremath{\sigma}}_{W}\ensuremath{\propto}{W}_{c}/N$. We use this result to predict the spin long-term evolution for a very large number of spins (N=50), inaccessible for exact diagonalization, and to suggest the rescaling of hopping interaction with the system size to attain the localization transition at finite disordering in the thermodynamic limit of infinite number of spins.
We investigate the stability of the many-body localized phase against quantum avalanche instabilities in a one-dimensional Heisenberg spin chain with long-range power-law interactions (V∝r−α). By combining exact diagonalization of static properties with Lindblad master equation simulations of open-system dynamics, we systematically map the interplay between interaction range and disorder strength. Our finite-size scaling analysis of entanglement entropy identifies a critical interaction exponent αc≈2, which separates a fragile regime, characterized by an exponentially diverging critical disorder, from a robust short-range regime. To rigorously test the system's resistance to avalanches, we couple the boundary to an infinite-temperature bath and track the propagation of the thermalization front into the localized bulk. We find that the characteristic thermalization time follows a unified scaling law, Trth∼exp[κ(α)LW] (herein, L is the system size, and W is the disorder intensity), which diverges exponentially with the product of system size and disorder strength. This suppression enables the derivation of a quantitative stability criterion, Wstab(α), representing the minimum critical disorder strength required to maintain avalanche stability. Our results confirm that the MBL phase remains asymptotically stable in the thermodynamic limit when disorder exceeds an interaction-dependent threshold, bridging theoretical debates on long-range MBL and providing a roadmap for observing these dynamics in experimental platforms such as Rydberg atom arrays.
Motivated by recent debates around the many-body localization (MBL) problem, and in particular its stability against systemwide resonances, we investigate long-distance spin-spin correlations across the phase diagram of the random-field XXZ model, with a particular focus on the strong disorder regime. Building on state-of-the-art shift-invert diagonalization techniques, we study the high-energy behavior of transverse and longitudinal correlation functions, computed at the largest possible distance, for a broad range of disorder and interaction strengths. Our results show that while transverse correlations display a fairly stable exponential decay over the entire XXZ phase diagram, longitudinal correlations exhibit markedly different behavior, revealing distinct physical regimes. More precisely, we identify an intermediate disorder region where standard observables show well-converged MBL behavior [Colbois et al., Phys. Rev. Lett. 133, 116502 (2024)] while the distributions of longitudinal correlations reveal unexpected fat tails toward large values. These rare events strongly influence the average decay of longitudinal correlations, which we find to be algebraic in a broad region inside the supposed MBL phase, whereas the typical decay remains mostly exponential. At stronger disorder and weaker interactions, this intermediate regime is replaced by a more conventional exponential decay with short correlation lengths for both typical and average correlators, as expected for standard localization. Our findings shed light on the systemwide instabilities and raise important questions about the impact of such rare but large long-range correlations on the stability of the MBL phase. Finally, we discuss the possible fate of the intermediate region in the context of recent perspectives in the field.