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.
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.
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) 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.
While the eigenstate entanglement entropy has been extensively studied for fermionic systems, much less is known about bosonic systems. Here, we study the entanglement entropy of mid-spectrum eigenstates of Bose-Hubbard models, focusing on weakly disordered models with and without particle-number conservation, and contrasting them with the translationally-invariant model. We analyze the volume-law and O(1) contributions to the entanglement entropy via the averages over mid-spectrum eigenstates and the corresponding distributions. We derive the volume-law coefficient of the entanglement entropy by generalizing the mean-field approach from [Phys. Rev. Lett. 119, 220603 (2017)] to many-body systems with a tunable local bosonic cutoff, which agrees with previous analytical and numerical results from [Phys. Rev. B 110, 235154 (2024)]. We show that the volume-law contribution to the entanglement entropy does not change upon breaking translational invariance via on-site disorder. We then numerically study the role of the subleading O(1) contribution to the entanglement entropy. We find that, in the particle-number conserving case, it exhibits a nontrivial dependence on the particle-number density and the local bosonic cutoff, while without particle-number conservation, results suggest the emergence of a universal O(1) contribution beyond the random pure state predictions.
We treat the disordered Heisenberg model in one dimension as the standard model of many-body localization (MBL). Two new and independent order parameters stemming solely from the half-chain von Neumann entanglement entropy SvN are introduced to probe the eigenstate phase transition in this model. From the symmetry-endowed entropy decomposition, they are the probability distribution deviation ∣d(pn)∣ and the von Neumann entropy SvNn(Dn=max) of the maximally dimensional symmetry subdivision. The finite-size scaling reveals that the probability set pn drives the localization transition, preceded by a thermalization breakdown transition governed by SvNn. For the noninteracting case, these transitions coincide, but in the interacting circumstance, they separate. Such separability creates an intermediate phase regime and discriminates between the Anderson and MBL transitions. One obstacle whose solution eludes the community to date concerns the violation of the Harris criterion in most numerical investigations of MBL. Upon elucidating the mutually independent measures comprising SvN, it becomes clear that the previous studies may lack the resolution to pinpoint, and thus potentially overlook, the crucial internal structure of the transition. We show that after this necessary decomposition, the universal critical exponents for both transitions of ∣d(pn)∣ and SvNn(Dn=max) fulfill the Harris criterion: $\ensuremath{\nu}\ensuremath{\approx}2\phantom{\rule{4pt}{0ex}}(\ensuremath{\nu}\ensuremath{\approx}1.5)$ for quench (quasirandom) disorder. Our work puts forth symmetry combined with entanglement as an organization principle for the generic eigenstate matter and phase transition.
The space of one-dimensional disordered interacting quantum models displaying a many-body localization (MBL) transition seems sufficiently rich to produce critical points with level statistics interpolating continuously between the Poisson statistics of the localized phase and the Wigner–Dyson statistics of the delocalized phase. In this paper, we consider the strong disorder limit of the MBL transition, where the level statistics at the MBL critical point is close to the Poisson statistics. We analyze a one-dimensional quantum spin model, in order to determine the statistical properties of the rare extensive resonances that are needed to destabilize the MBL phase. At criticality, we find that the entanglement entropy can grow with an exponent 0 < α < 1 anywhere between the area law α = 0 and the volume law α = 1 , as a function of the resonances properties, while the entanglement spectrum follows the strong multifractality statistics. In the MBL phase near criticality, we obtain the simple value ν = 1 for the correlation length exponent. Independently of the strong disorder limit, we explain why, for the many-body localization transition concerning individual eigenstates, the correlation length exponent ν is not constrained by the usual Harris inequality ν ≥ 2 / d , so that there is no theoretical inconsistency with the best numerical measure ν = 0 . 8 ( 3 ) obtained by Luitz et al. (2015).