It has been proposed that the states of fully many-body localized systems can be described in terms of conserved local pseudospins. Due to the multitude of ways to define these, the explicit identification of the optimally local pseudospins in specific systems is non-trivial. Given continuing intense interest in the role of disorder in strongly correlated systems, we consider the disordered Hubbard model. By studying a small system we provide concrete examples of the form of local integrals of motion in the Anderson-Hubbard model. Moreover, we are able not only to identify the most local choice but also to explore the nature of the distribution of possible choices. We track the evolution of the optimally localized pseudospins as hopping and interactions are varied to move the system away from the trivially localized atomic limit.
It has been proposed that the states of fully many-body localized systems can be described in terms of conserved local pseudospins. Due to the multitude of ways to define these, the explicit identification of the optimally local pseudospins in specific systems is non-trivial. Given continuing intense interest in the role of disorder in strongly correlated systems, we consider the disordered Hubbard model. By studying a small system we provide concrete examples of the form of local integrals of motion in the Anderson-Hubbard model. Moreover, we are able not only to identify the most local choice but also to explore the nature of the distribution of possible choices. We track the evolution of the optimally localized pseudospins as hopping and interactions are varied to move the system away from the trivially localized atomic limit.
We study the interplay of disorder localization and strong local interactions within the Anderson-Hubbard model. Taking into account local Mott-Hubbard physics and static screening of the disorder potential, the system is mapped onto an effective single-particle Anderson model, which is studied within the self-consistent theory of electron localization. For fermions, we find rich nonmonotonic behavior of the localization length ξ, particularly in two-dimensional systems, including an interaction-induced exponential enhancement of ξ for small and intermediate disorders and a strong reduction of ξ due to hopping suppression by strong interactions. In three dimensions, we identify for half filling a Mott-Hubbard-assisted Anderson localized phase existing between the metallic and the Mott-Hubbardgapped phases. For small U there is re-entrant behavior from the Anderson localized phase to the metallic phase. For bosons, the unrestricted particle occupation number per lattice site yields a monotonic enhancement of ξ as a function of decreasing interaction, which we assume to persist until the superfluid Bose- Einstein condensate phase is entered. Besides, we study cold atomic gases expanding, by a diffusion process, in a weak random potential. We show that the density-density correlation function of the expanding gas is strongly affected by disorder and we estimate the typical size of a speckle spot, i.e., a region of enhanced or depleted density. Both a Fermi gas and a Bose-Einstein condensate (in a mean-field approach) are considered.
Anderson localization is known to be inevitable in one-dimension for generic disordered models. Since localization leads to Poissonian energy level statistics, we ask if localized systems possess 'additional' integrals of motion as well, so as to enhance the analogy with quantum integrable systems. We answer this in the affirmative in the present work. We construct a set of nontrivial integrals of motion for Anderson localized models, in terms of the original creation and annihilation operators. These are found as a power series in the hopping parameter. The recently found Type-1 Hamiltonians, which are known to be quantum integrable in a precise sense, motivate our construction. We note that these models can be viewed as disordered electron models with infinite-range hopping, where a similar series truncates at the linear order. We show that despite the infinite range hopping, all states but one are localized. We also study the conservation laws for the disorder free Aubry–Andre model, where the states are either localized or extended, depending on the strength of a coupling constant. We formulate a specific procedure for averaging over disorder, in order to examine the convergence of the power series. Using this procedure in the Aubry–Andre model, we show that integrals of motion given by our construction are well-defined in localized phase, but not so in the extended phase. Finally, we also obtain the integrals of motion for a model with interactions to lowest order in the interaction.
Although many-body localization has primarily been studied in systems with a single local degree of freedom, experimental studies of many-body localization in cold atom systems motivate exploration of the disordered Hubbard model. With two coupled local degrees of freedom, it is natural to ask how localization in spin relates to disorder in charge and vice versa. Most prior work has addressed disorder in only one of these sectors and often has not used measures of localization which distinguish between charge and spin. Here, we explore localization in the Hubbard model with a wide range of independent values of charge and spin disorder using measures of localization based on charge and spin-specific integrals of motion. Our results demonstrate symmetry between the response of the spin to charge disorder and vice versa, and we find very weak disorder in one channel, so long as the disorder in the other channel is sufficiently strong, results in localization in both channels. The strength of disorder required in the less-disordered channel declines as the system size increases. Furthermore, the weaker the disorder in the less-disordered channel, the longer the timescale at which localization appears in the dynamics of this degree of freedom.
We develop a procedure which systematically generates all conserved operators in the disordered models of interacting fermions. Among these operators, we identify and count the independent and local integrals of motion (LIOM), which represent the hallmark of the many-body localization (MBL). The method is tested first on the prototype disordered chain of interacting spinless fermions. As expected for full MBL, we find for large enough disorder ${N}_{M}={2}^{M}\ensuremath{-}1$ independent and quasilocal LIOM with support on M consecutive sites. On the other hand, the study of the disordered Hubbard chain reveals that ${3}^{M}\ensuremath{-}1<{N}_{M}\ensuremath{\lesssim}{4}^{M}/2$, which is less than required for full MBL but much more than in the case of spinless fermions.
We study the ground-state phase diagram of the Anderson-Hubbard model with correlated hopping at half-filling in one dimension. The Hamiltonian has a local Coulomb repulsion U and a disorder potential with local energies randomly distributed in the interval $(\ensuremath{-}W,+W)$ with equal probability, acting on the singly occupied sites. The hopping process which modifies the number of doubly occupied sites is forbidden. The hopping between nearest-neighbor singly occupied and empty sites or between singly occupied and doubly occupied sites has the same amplitude t. We identify three different phases as functions of the disorder amplitude W and Coulomb interaction strength U>0. When U<4t the system shows a metallic phase: (i) only when no disorder is present W=0 or an Anderson-localized phase, (ii) when disorder is introduced $W\ensuremath{\ne}0$. When U>4t the Anderson-localized phase survives as long as disorder effects dominate on the interaction effects, otherwise a Mott-insulator phase (iii) arises. Phases (i) and (ii) are characterized by a finite density of doublons and a vanishing charge gap among the ground state and the excited states. Phase (iii) is characterized by the vanishing density of doublons and a finite gap for the charge excitations.
In this thesis we consider a variational wave function approach as a possible route to describe the competition between disorder and strong electron-electron interaction in two dimensions. In particular we aim to obtain a transparent and physically intuitive understanding of the competition between these two localizing forces within the simplest model where they both are active, namely the disordered Hubbard model at half filling and in a square lattice. Our approach is based on an approximate form of the ground-state wave function, which we believe contains the physically relevant ingredients for a correct description of both the Mott and the Anderson insulators, where electrons are localized by the Coulomb repulsion and by disorder, respectively. For strongly interacting fermionic systems, a standard variational wave function is constructed by a correlation term acting on a Slater determinant, the latter being an uncorrelated metallic state. Previous variational calculations showed that a long-range density-density correlation factor, so called Jastrow factor, is needed to correctly describe the Mott insulator [9]. This term, which is collective by definition, correlates spatially charge uctuations, thus preventing their free motion that would otherwise imply metallic conductance. For this reason, our variational wave function does include such a term. Anderson localization is instead mostly a matter of single-particle wave functions, hence it pertains to the uncorrelated Slater determinant which the Jastrow factor acts onto. We consider both the case in which we enforce paramagnetism in the wave function and the case in which we allow for magnetic ordering. Summarizing briefly our results, we find that, when the variational wave function is forced to be paramagnetic, the Anderson insulator to Mott insulator transition is continuous. This transition can be captured by studying several quantities. In particular, a novel one that we have identified and that is easily accessible variationally is the disconnected density-density fluctuation at long wavelength, defined by lim where ^nq is the Fourier transform of the charge density at momentum q, (...) denotes quantum average at fixed disorder and the overbar represents the average over disorder configurations. We find that Ndisc q!0 is everywhere finite in the Anderson insulator and vanishes critically at the Mott transition, staying zero in the Mott insulator. When magnetism is allowed and the hopping only connects nearest neighbor sites, upon increasing interaction the paramagnetic Anderson insulator first turns antiferromagnetic and finally the magnetic and compressible Anderson insulator gives way to an incompressible antiferromagnetic Mott insulator. The optimized uncorrelated Slater determinant is always found to be the eigenstate of a disordered non-interacting effective Hamiltonian, which suggests that the model is never metallic. Finally, when magnetism is frustrated by a next to nearest neighbor hopping, the overall sequence of phases does not change. However, the paramagnetic to magnetic transition within the Anderson insulator basin of stability turns first order. Indeed, within the magnetically ordered phase, we find many almost degenerate paramagnetic states with well defined local moments. This is suggestive of an emerging glassy behavior when the competition between disorder and strong correlation is maximum.