Many-body localization of a disorder interacting boson system in one dimension is studied numerically at the filling factor being one-half, in terms of level statistics, local compressibility, correlation function and entanglement entropies. The von Neumann entanglement entropy is decoupled into a particle number entropy and a configuration entropy. An anomalous volume-law behavior is found for the configuration entanglement entropy to confirm a recent experimental observation [A. Lukin, M. Rispoli, R. Schittko, et al., Science 364, 256 (2019)] for sufficient strong disorder, while the particle number entropy fulfills an area-law corresponding to the total entropy for disordered spin chain. The localization length are extracted from a two-body correlation function for many-body localization states and corresponding time-evolutions states as well. A phase diagrams is established with consisting of an ergodic thermalized region and a many-body-localization region in a parameter space of the disorder strength and the energy density. Two regions are separated by a many-body mobility edge deducted from the standard deviation of the particle-number entanglement entropy, which appears consistent with that based on the localization length. Slow dynamics characterized by a logarithmic time-dependence is explicitly shown for both the particle number entropy and the configuration entropy in an intermediate regime of their time-evolutions, which does not show up in the Anderson localization case, i.e. non-interacting disorder systems.
Many-body localization of a disorder interacting boson system in one dimension is studied numerically at the filling factor being one-half, in terms of level statistics, local compressibility, correlation function and entanglement entropies. The von Neumann entanglement entropy is decoupled into a particle number entropy and a configuration entropy. An anomalous volume-law behavior is found for the configuration entanglement entropy to confirm a recent experimental observation [A. Lukin, M. Rispoli, R. Schittko, et al., Science 364, 256 (2019)] for sufficient strong disorder, while the particle number entropy fulfills an area-law corresponding to the total entropy for disordered spin chain. The localization length are extracted from a two-body correlation function for many-body localization states and corresponding time-evolutions states as well. A phase diagrams is established with consisting of an ergodic thermalized region and a many-body-localization region in a parameter space of the disorder strength and the energy density. Two regions are separated by a many-body mobility edge deducted from the standard deviation of the particle-number entanglement entropy, which appears consistent with that based on the localization length. Slow dynamics characterized by a logarithmic time-dependence is explicitly shown for both the particle number entropy and the configuration entropy in an intermediate regime of their time-evolutions, which does not show up in the Anderson localization case, i.e. non-interacting disorder systems.
Many-body localization (MBL) of a disordered interacting boson system in one dimension is studied numerically at the filling faction one-half. The von Neumann entanglement entropy SvN is commonly used to detect the MBL phase transition but remains challenging to be directly measured. Based on the U(1) symmetry from the particle number conservation, SvN can be decomposed into the particle number entropy SN and the configuration entropy SC. In light of the tendency that the eigenstate’s SC nears zero in the localized phase, we introduce a quantity describing the deviation of SN from the ideal thermalization distribution; finite-size scaling analysis illustrates that it shares the same phase transition point with SvN but displays the better critical exponents. This observation hints that the phase transition to MBL might largely be determined by SN and its fluctuations. Notably, the recent experiments [A. Lukin, et al., Science 364, 256 (2019); J. Léonard, et al., Nat. Phys. 19, 481 (2023)] demonstrated that this deviation can potentially be measured through the SN measurement. Furthermore, our investigations reveal that the thermalized states primarily occupy the low-energy section of the spectrum, as indicated by measures of localization length, gap ratio, and energy density distribution. This low-energy spectrum of the Bose model closely resembles the entire spectrum of the Fermi (or spin XXZ) model, accommodating a transition from the thermalized to the localized states. While, owing to the bosonic statistics, the high-energy spectrum of the model allows the formation of distinct clusters of bosons in the random potential background. We analyze the resulting eigenstate properties and briefly summarize the associated dynamics. To distinguish between the phase regions at the low and high energies, a probing quantity based on the structure of SvN is also devised. Our work highlights the importance of symmetry combined with entanglement in the study of MBL. In this regard, for the disordered Heisenberg XXZ chain, the recent pure eigenvalue analyses in [J. Suntajs, et al., Phys. Rev. E 102, 062144 (2020)] would appear inadequate, while methods used in [A. Morningstar, et al., Phys. Rev. B 105, 174205 (2022)] that spoil the U(1) symmetry could be misleading.
The phenomenology of many-body localization (MBL) develops mainly from tackling the one-dimensional spin or fermion systems. The situation when interacting bosons get clustered in a random potential remains less explored. Using numerical quantum quenches with special emphasis on the integration of both symmetry and energy resolutions, we comprehensively study the dynamics of symmetry-resolved entanglement in a disordered Bose-Hubbard (dBH) model, concentrating on the two types of inhomogeneous initial states to target the lower- and higher-energy sections of its dynamical phase diagram. (i) Motivated by the recent experiment [A. Lukin et al., Science 364, 256 (2019)] which focused on the lower-energy dynamical behaviors of the dBH chain, we first show that at low energies, for a thermalizing state, although the second law of thermodynamics prohibits the decrease of the total entropy over time, for part of the channel-resolved entropies, a long-term entropic reduction may arise at weak disorder. (ii) A companion channel-resolved analysis at strong disorder further hints that the prior observed double-logarithmic growth of the number entropy might not directly indicate the breakdown of MBL in spin or fermion chains, providing a refreshing perspective on this major controversy in the community. (iii) From time-evolving the line-shape low-energy product state, we subsequently reveal an abrupt formation of a novel ``entropy imbalance pattern'' across the different symmetry channels. Intriguingly, this imbalance melts in the strong-disorder limit. We conjecture that the melting of the entropic pattern, together with the freezing of a concurrent particle-density wave, embodies a dual trait inherent to MBL. (iv) Conversely, the higher-energy section of the dynamical phase diagram is where the dBH model differs most significantly from the spin or fermion systems. This parametric space was not included in the previous literature. Specifically, we find a cluster MBL regime, unique to the Bose statistics, emerging from the higher-energy section. This cluster MBL regime, realizable even at weak disorder, does not appear to suffer from the finite-size drift and is distinguished by its absence of the hallmark of MBL---the unbounded growth of the entanglement entropy. Our theoretical predictions are, by and large, testable via the present experimental facilities.
We investigate the non-equilibrium dynamics of the disordered one-dimensional Fermi-Hubbard model with a focus on many-body localization. The system is initialized in a charge-density-wave-state, and its time evolution is analyzed through sublattice imbalance (spin and charge), and bipartite entanglement entropy. A clear crossover from ergodic to non-ergodic behavior is observed with increasing disorder strength. In the weak disorder regime, rapid decay of imbalance and the fast growth of entanglement indicate efficient thermalization. In contrast, a strong disorder leads to persistent imbalance and slow dynamics, signaling the breakdown of ergodicity. The charge and spin sectors exhibit distinct relaxation behavior, providing evidence for partial decoupling between these degrees of freedom. Furthermore, in the interacting regime, the entanglement entropy shows slow logarithmic growth, reflecting the dephasing-driven dynamics characteristic of the many-body localized phase. These results highlight the interplay between disorder and interactions in determining the dynamical properties of the system and establish robust signatures of many-body localization in the Fermi-Hubbard model.
Many-body localization (MBL) has been widely investigated for both fermions and bosons, it is, however, much less explored for anyons. Here we numerically calculate several physical characteristics related to MBL of a one-dimensional disordered anyon-Hubbard model in both localized and delocalized regions. We figure out a logarithmically slow growth of the half-chain entanglement entropy and an area-law rather than volume-law obedience for the highly excited eigenstates in the MBL phase. The adjacent energy level gap-ratio parameter is calculated and is found to exhibit a Poisson-like probability distribution in the deep MBL phase. By studying a hybridization parameter, we reveal an intriguing effect that the statistics can induce localization-delocalization transition. Several physical quantities, such as the half-chain entanglement, the adjacent energy level gap-ratio parameter, the long-time limit of the particle imbalance, and the critical disorder strength, are shown to be nonmonotonically dependent on the anyon statistical angle. Furthermore, a feasible scheme based on the spectroscopy of energy levels is proposed for the experimental observation of these statistically related properties.
Many-body localization for a system of bosons trapped in a one-dimensional lattice is discussed. Two models that may be realized for cold atoms in optical lattices are considered. The model with a random on-site potential is compared with previously introduced random interactions model. While the origin and character of the disorder in both systems is different they show interesting similar properties. In particular, many-body localization appears for a sufficiently large disorder as verified by a time evolution of initial density wave states as well as using statistical properties of energy levels for small system sizes. Starting with different initial states, we observe that the localization properties are energy-dependent which reveals an inverted many-body localization edge in both systems (that finding is also verified by statistical analysis of energy spectrum). Moreover, we consider computationally challenging regime of transition between many body localized and extended phases where we observe a characteristic algebraic decay of density correlations which may be attributed to subdiffusion (and Griffiths-like regions) in the studied systems. Ergodicity breaking in the disordered Bose–Hubbard models is compared with the slowing-down of the time evolution of the clean system at large interactions.
Many-body localization occurs in isolated quantum systems when Anderson localization persists in the presence of finite interactions. Despite strong evidence for the existence of a many-body localization transition, a reliable extraction of the critical disorder strength is difficult due to a large drift with system size in the studied quantities. In this Letter, we explore two entanglement properties that are promising for the study of the many-body localization transition: the variance of the half-chain entanglement entropy of exact eigenstates and the long time change in entanglement after a local quench from an exact eigenstate. We investigate these quantities in a disordered quantum Ising chain and use them to estimate the critical disorder strength and its energy dependence. In addition, we analyze a spin-glass transition at large disorder strength and provide evidence for it being a separate transition. We, thereby, give numerical support for a recently proposed phase diagram of many-body localization with localization protected quantum order [Huse et al., Phys. Rev. B 88, 014206 (2013).