We develop a numerical technique to study Anderson localization in interacting electronic systems. The ground state of the disordered system is calculated with quantum Monte Carlo simulations while the localization properties are extracted from the "Thouless conductance" g, i.e., the curvature of the energy with respect to an Aharonov-Bohm flux. We apply our method to polarized electrons in a two-dimensional system of size L. We recover the well-known universal beta(g)=dlogg/dlogL one parameter scaling function without interaction. Upon switching on the interaction, we find that beta(g) is unchanged while the system flows toward the insulating limit. We conclude that polarized electrons in two dimensions stay in an insulating state in the presence of weak to moderate electron-electron correlations.
We develop a numerical technique to study Anderson localization in interacting electronic systems. The ground state of the disordered system is calculated with quantum Monte Carlo simulations while the localization properties are extracted from the "Thouless conductance" g, i.e., the curvature of the energy with respect to an Aharonov-Bohm flux. We apply our method to polarized electrons in a two-dimensional system of size L. We recover the well-known universal beta(g)=dlogg/dlogL one parameter scaling function without interaction. Upon switching on the interaction, we find that beta(g) is unchanged while the system flows toward the insulating limit. We conclude that polarized electrons in two dimensions stay in an insulating state in the presence of weak to moderate electron-electron correlations.
We have studied the effects of interactions on persistent currents in half-filled and quarter-filled Hubbard models with weak and intermediate strength disorder. Calculations are performed using a variational Gutzwiller ansatz that describes short-range correlations near the Mott transition. We apply an Aharonov-Bohm magnetic flux, which generates a persistent current that can be related to the Thouless conductance. The magnitude of the current depends on both the strength of the screened disorder potential and the strength of electron-electron correlations, and the Anderson localization length can be extracted from the scaling of the current with system size. At half-filling, the persistent current is reduced by strong correlations when the interaction strength is large. Surprisingly, we find that the disorder potential is strongly screened in the large interaction limit, so that the localization length grows with increasing interaction strength even as the magnitude of the current is suppressed. This supports earlier dynamical mean-field-theory predictions that the elastic scattering rate is suppressed near the Mott transition.
Disorder and interactions both play crucial roles in quantum transport. Decades ago, Mott showed that electron-electron interactions can lead to insulating behavior in materials that conventional band theory predicts to be conducting. Soon thereafter, Anderson demonstrated that disorder can localize a quantum particle through the wave interference phenomenon of Anderson localization. Although interactions and disorder both separately induce insulating behavior, the interplay of these two ingredients is subtle and often leads to surprising behavior at the periphery of our current understanding. Modern experiments probe these phenomena in a variety of contexts (e.g. disordered superconductors, cold atoms, photonic waveguides, etc.); thus, theoretical and numerical advancements are urgently needed. In this thesis, we report progress on understanding two contexts in which the interplay of disorder and interactions is especially important. The first is the so-called “dirty” or random boson problem. In the past decade, a strong-disorder renormalization group (SDRG) treatment by Altman, Kafri, Polkovnikov, and Refael has raised the possibility of a new unstable fixed point governing the superfluid-insulator transition in the one-dimensional dirty boson problem. This new critical behavior may take over from the weak-disorder criticality of Giamarchi and Schulz when disorder is sufficiently strong. We analytically determine the scaling of the superfluid susceptibility at the strong-disorder fixed point and connect our analysis to recent Monte Carlo simulations by Hrahsheh and Vojta. We then shift our attention to two dimensions and use a numerical implementation of the SDRG to locate the fixed point governing the superfluid-insulator transition there. We identify several universal properties of this transition, which are fully independent of the microscopic features of the disorder. The second focus of this thesis is the interplay of localization and interactions in systems with high energy density (i.e., far from the usual low energy limit of condensed matter physics). Recent theoretical and numerical work indicates that localization can survive in this regime, provided that interactions are sufficiently weak. Stronger interactions can destroy localization, leading to a so-called many-body localization transition. This dynamical phase transition is relevant to questions of thermalization in isolated quantum systems: it separates a many-body localized phase, in which localization prevents transport and thermalization, from a conducting (“ergodic”) phase in which the usual assumptions of quantum statistical mechanics hold. Here, we present evidence that many-body localization also occurs in quasiperiodic systems that lack true disorder.
Using a combination of ground state quantum Monte Carlo and finite size scaling techniques, we perform a systematic study of the effect of Coulomb interaction on the localization length of a disordered two-dimensional electron gas. We find that correlations delocalize the 2D system. In the absence of valley degeneracy (as in GaAs heterostructures), this delocalization effect corresponds to a finite increase of the localization length. The delocalization is much more dramatic in the presence of valley degeneracy [as in Si metal-oxide-semiconductor field-effect transistors], where the localization length increases drastically. We find that a simple mechanism accounts for the main features of the metallic behavior observed in two-dimensional gases. Our findings support the claim that this behavior is a genuine effect of the presence of interactions, yet that the system is not a "true" metal in the thermodynamic sense.
We study the interplay of disorder and interaction effects including bosonic degrees of freedom in the framework of a generic one-dimensional transport model: the Anderson-Edwards model. Using the density-matrix-renormalization group technique, we extract the localization length and the renormalization of the Tomonaga-Luttinger-liquid parameter from the charge-structure factor by a elaborate sample-average finite-size scaling procedure. The properties of the Anderson localized state can be described in terms of scaling relations of the metallic phase without disorder. We analyze how disorder competes with the charge-density-wave correlations triggered by the bosons and give evidence that disorder will destroy the long-range charge-ordered state.
Disorder or sufficiently strong interactions can render a metallic state unstable, causing it to turn into an insulating one. Despite the fact that the interplay of these two routes to a vanishing conductivity has been a central research topic, a unifying picture has not emerged so far. Here, we establish that the two-dimensional Falicov-Kimball model, one of the simplest lattice models of strong electron correlation, does allow for the study of this interplay. In particular, we show that this model at particle-hole symmetry possesses three distinct thermodynamic insulating phases and exhibits Anderson localization. The previously reported metallic phase is identified as a finite-size feature due to the presence of weak localization. We characterize these phases by their electronic density of states, staggered occupation, conductivity, and the generalized inverse participation ratio. The implications of our findings for other strongly correlated systems are discussed.
We present, in this dissertation, a numerical simulation method to study interacting fermion systems. The general simulation procedures are discussed in connection with a description of a Quantum Monte Carlo simulation algorithm for interacting electrons on lattices. The algorithm presented here has been used to simulate interacting electrons on lattices, and it makes possible the study of substantially larger systems than can be studied by other numerical methods. As long as certain limits of applicability are respected, model Hamiltonians of interacting electrons can be studied without resort to uncontrolled approximations. The method then provides nearly exact solutions to model Hamiltonians of many-body systems, in the sense that the degree of error can be controlled. We also discuss some results obtained from simulations of the extended Hubbard model and and the periodic Anderson model. For the extended Hubbard model in two dimensions, different regions are identified where electron correlation produces antiferromagnetism, charge-density-waves, and singlet pairing superconducting behaviors. We also find regions where transitions from one type of correlation to another occur. For the symmetric periodic Anderson model, we observe the formation of local spin moments at high temperatures and their quenching at low temperatures, as in the single impurity Anderson model. In addition, we find antiferromagnetic interaction between f local moments at low temperatures which are not present in a single impurity problem.