We study the potential influence of the particle multi-occupations on the stability of many-body localization in the disordered Bose-Hubbard model. Within the higher-energy section of the dynamical phase diagram, we find that there is no apparent finite-size boundary drift between the thermal phase and the many-body localized regime. We substantiate this observation by introducing the Van Vleck perturbation theory into the field of many-body localization. The appropriateness of this method rests largely on the peculiar Hilbert-space structure enabled by the particles' Bose statistics. The situation is reversed in the lower-energy section of the dynamical phase diagram, where the significant finite-size boundary drift pushes the putative many-body localized regime up to the greater disorder strengths. We utilize the algebraic projection method to make a connection linking the disordered Bose-Hubbard model in the lower-energy section to an intricate disordered spin chain model. This issue of the finite-size drift could hence be analogous to what happens in the disordered Heisenberg chain. Both trends might be traced back to the particles' intrinsic or emergent single-occupancy constraint such as the spin-$1/2$, hard-core boson, or spinless fermion degrees of freedom.
We study the potential influence of the particle multi-occupations on the stability of many-body localization in the disordered Bose-Hubbard model. Within the higher-energy section of the dynamical phase diagram, we find that there is no apparent finite-size boundary drift between the thermal phase and the many-body localized regime. We substantiate this observation by introducing the Van Vleck perturbation theory into the field of many-body localization. The appropriateness of this method rests largely on the peculiar Hilbert-space structure enabled by the particles' Bose statistics. The situation is reversed in the lower-energy section of the dynamical phase diagram, where the significant finite-size boundary drift pushes the putative many-body localized regime up to the greater disorder strengths. We utilize the algebraic projection method to make a connection linking the disordered Bose-Hubbard model in the lower-energy section to an intricate disordered spin chain model. This issue of the finite-size drift could hence be analogous to what happens in the disordered Heisenberg chain. Both trends might be traced back to the particles' intrinsic or emergent single-occupancy constraint such as the spin-1/2, hard-core boson, or spinless fermion degrees of freedom.
We numerically investigate one-dimensional Bose-Hubbard chains with on-site disorder by means of exact diagonalization. A primary focus of our work is on characterizing Fock-space localization in this model from the single-particle perspective. For this purpose, we compute the one-particle density matrix (OPDM) in many-body eigenstates. We show that the natural orbitals (the eigenstates of the OPDM) are extended in the ergodic phase and real-space localized when one enters the many-body localization (MBL) phase. Furthermore, the distributions of occupations of the natural orbitals can be used as measures of Fock-space localization in the respective basis. Consistent with previous studies, we observe signatures of a transition from the ergodic to the MBL regime when increasing the disorder strength. We further demonstrate that Fock-space localization, albeit weaker, is also evidently present in the distribution of the physical densities in the MBL regime, for both soft- and hard-core bosons. Moreover, the full distribution of the densities of the physical particles provides a one-particle measure for the detection of the ergodic-MBL transition which could be directly accessed in experiments with ultracold gases.
The many-body localization (MBL) proximity effect is an intriguing phenomenon where a thermal bath localizes due to the interaction with a disordered system. The interplay of thermal and nonergodic behavior in these systems gives rise to a rich phase diagram, whose exploration is an active field of research. In this paper, we study a bosonic Hubbard model featuring two particle species representing the bath and the disordered system. Using state-of-the-art numerical techniques, we investigate the dynamics of the model in different regimes, based on which we obtain a tentative phase diagram as a function of coupling strength and bath size. When the bath is composed of a single particle, we observe clear signatures of a transition from an MBL proximity effect to a delocalized phase. Increasing the bath size, however, its thermalizing effect becomes stronger and eventually the whole system delocalizes in the range of moderate interaction strengths studied. In this regime, we characterize particle transport, revealing diffusive behavior of the originally localized bosons.
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.
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 study the $t\text{\ensuremath{-}}V$ disordered spinless fermionic chain in the strong-coupling regime, $t/V\ensuremath{\rightarrow}0$. Strong interactions highly hinder the dynamics of the model, fragmenting its Hilbert space into exponentially many blocks in system size. Macroscopically, these blocks can be characterized by the number of new degrees of freedom, which we refer to as movers. We focus on two limiting cases: blocks with only one mover and ones with a finite density of movers. The former many-particle block can be exactly mapped to a single-particle Anderson model with correlated disorder in one dimension. As a result, these eigenstates are always localized for any finite amount of disorder. The blocks with a finite density of movers, on the other side, show a many-body localized (MBL) transition that is tuned by the disorder strength. Moreover, we provide numerical evidence that its ergodic phase is diffusive at weak disorder. Approaching the MBL transition, we observe subdiffusive dynamics at finite timescales and find indications that this might be only a transient behavior before crossing over to diffusion.
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.