Doped Hubbard Model Research Articles

Strange metal in the doped Hubbard model via percolation

Strange metal in the doped Hubbard model via percolation

Origin and fate of the pseudogap in the doped Hubbard model

The relationship between the pseudogap and underlying ground-state phases has not yet been rigorously established. We investigated the doped two-dimensional Hubbard model at finite temperature using controlled diagrammatic Monte Carlo calculations, allowing for the computation of spectral properties in the infinite-size limit and with arbitrary momentum resolution. We found three distinct regimes as a function of doping and interaction strength: a weakly correlated metal, a correlated metal with strong interaction effects, and a pseudogap regime at low doping. We show that the pseudogap forms both at weak coupling, when the magnetic correlation length is large, and at strong coupling, when it is shorter. As the temperature goes to zero, the pseudogap regime extrapolates precisely to the ordered stripe phase found by ground-state methods.

Nagaoka ferromagnetism in doped Hubbard models in optical lattices

Nagaoka ferromagnetism in doped Hubbard models in optical lattices

Tilted stripes origin in La1.88Sr0.12CuO4 revealed by anisotropic next-nearest neighbor hopping

Spin- and charge- stripe order has been extensively studied in the superconducting cuprates, among which underdoped La2−xSrxCuO4 (LSCO) is an archetype with static spin stripes at low temperatures. An intriguing, but not completely understood, phenomenon in LSCO is that the stripes are tilted away from the high-symmetry Cu-Cu directions. Using high-resolution neutron scattering on LSCO with x = 0.12, we find two coexisting phases at low temperatures, one with static spin stripes and the other with fluctuating ones, both sharing the same tilt angle. Our numerical calculations using the doped Hubbard model elucidate the tilting’s origin, attributing it to anisotropic next-nearest neighbor hopping t′\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${t}^{{\\prime} }$$\\end{document}, consistent with the material’s slight orthorhombicity. Our results underscore the model’s success in describing specific details of the ground state of this real material and highlight the role of t′\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${t}^{{\\prime} }$$\\end{document} in the Hamiltonian, revealing the delicate interplay between stripes and superconductivity across theoretical and experimental contexts.

Dynamical t/U Expansion of the Doped Hubbard Model

We construct a new U(1) slave-spin representation for the single-band Hubbard model in the large-U limit. The mean-field theory in this representation is more amenable to describe both the spin-charge-separation physics of the Mott insulator at half-filling and the strange metal behavior at finite doping. By employing a dynamical Green’s function theory for slave spins, we calculate the single-particle spectral function of electrons. The result is comparable to that in dynamical mean field theories. We then formulate a dynamical t/U expansion for the doped Hubbard model that reproduces the mean-field results at the lowest order of expansion. To the next order of expansion, it naturally yields an effective low-energy theory of a t–J model for spinons self-consistently coupled to an XXZ model for the slave spins. We show that the superexchange J is renormalized by doping, in agreement with the Gutzwiller approximation. Surprisingly, we find a new ferromagnetic channel of exchange interactions which survives in the infinite U limit, as a manifestation of the Nagaoka ferromagnetism.

Perturbative solution of fermionic sign problem in quantum Monte Carlo computations

We have developed a strong-coupling perturbation scheme for a general doped Hubbard model around a particle-hole-symmetric reference system, which is free from the fermionic sign problem. Our approach is based on the lattice determinantal Quantum Monte Carlo (QMC) method in both continuous and discrete time versions for large periodic clusters in a fermionic bath. By considering the first-order perturbation in the shift of the chemical potential and the second-neighbor hopping, we are able to obtain an accurate electronic spectral function for a range of parameters that correspond to optimally doped cuprate systems at temperatures of up to T = 0.1t, which are challenging to access using straightforward lattice QMC calculations. We also discuss the formation of the pseudogap and the nodal-antinodal dichotomy for a doped Hubbard system with the interaction parameter U equal to the bandwidth and the optimal value of the next-nearest-neighbor hopping parameter t′\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${t}^{{\\prime} }$$\\end{document} for high-temperature superconducting cuprates.

Ground-state phase diagram and superconductivity of the doped Hubbard model on six-leg square cylinders

Ground-state phase diagram and superconductivity of the doped Hubbard model on six-leg square cylinders

Hall coefficient and resistivity in the doped bilayer Hubbard model

Hall coefficient and resistivity in the doped bilayer Hubbard model

The Wiedemann-Franz law in doped Mott insulators without quasiparticles.

Many metallic quantum materials display anomalous transport phenomena that defy a Fermi liquid description. Here, we use numerical methods to calculate thermal and charge transport in the doped Hubbard model and observe a crossover separating high- and low-temperature behaviors. Distinct from the behavior at high temperatures, the Lorenz number [Formula: see text] becomes weakly doping dependent and less sensitive to parameters at low temperatures. At the lowest numerically accessible temperatures, [Formula: see text] roughly approaches the Wiedemann-Franz constant [Formula: see text], even in a doped Mott insulator that lacks well-defined quasiparticles. Decomposing the energy current operator indicates a compensation between kinetic and potential contributions, which may help to clarify the interpretation of transport experiments beyond Boltzmann theory in strongly correlated metals.

Quantifying hole-motion-induced frustration in doped antiferromagnets by Hamiltonian reconstruction

Unveiling the microscopic origins of quantum phases dominated by the interplay of spin and motional degrees of freedom constitutes one of the central challenges in strongly correlated many-body physics. When holes move through an antiferromagnetic spin background, they displace the positions of spins, which induces effective frustration in the magnetic environment. However, a concrete characterization of this effect in a quantum many-body system is still an unsolved problem. Here we present a Hamiltonian reconstruction scheme that allows for a precise quantification of hole-motion-induced frustration. We access non-local correlation functions through projective measurements of the many-body state, from which effective spin-Hamiltonians can be recovered after detaching the magnetic background from dominant charge fluctuations. The scheme is applied to systems of mixed dimensionality, where holes are restricted to move in one dimension, but SU(2) superexchange is two-dimensional. We demonstrate that hole motion drives the spin background into a highly frustrated regime, which can quantitatively be described by an effective J1–J2-type spin model. We exemplify the applicability of the reconstruction scheme to ultracold atom experiments by recovering effective spin-Hamiltonians of experimentally obtained 1D Fermi-Hubbard snapshots. Our method can be generalized to fully 2D systems, enabling promising microscopic perspectives on the doped Hubbard model.

Competition between antiferromagnetism and superconductivity in a doped Hubbard model with anisotropic interaction

The competition between antiferromagnetism and superconductivity is one of the central questions in the research of strong correlated systems. In this work, we utilize a double layer model containing Hubbard interaction and interlayer Heisenberg interaction to reveal their competitions. This model is free of sign problem at certain conditions, and we perform projector quantum Monte Carlo simulations to extract the ground state correlations of magnetism and superconductivity. Our results shows that the superconductivity emerges when the antiferromagnetism is suppressed by tuning the filling or the anisotropy of the interlayer Heisenberg interaction. This model can be seen as an analogue of unconventional superconductors and may help us to understand the transition from an antiferromagnetic insulator to a superconductor.

Enhanced superconductivity by near-neighbor attraction in the doped extended Hubbard model

Recent experiment has unveiled an anomalously strong electron-electron attraction in one-dimensional copper-oxide chain Ba$_{2-x}$Sr$_x$CuO$_{3+\delta}$. While the near-neighbor electron attraction $V$ in the one-dimensional extended Hubbard chain has been examined recently, its effect in the Hubbard model beyond the one-dimensional chain remains unclear. We report a density-matrix renormalization group study of the extended Hubbard model on long four-leg cylinders on the square lattice. We find that the near-neighbor electron attraction $V$ can notably enhance the long-distance superconducting correlations while simultaneously suppressing the charge-density-wave correlations. Specifically, for a modestly strong electron attraction, the superconducting correlations become dominant over the CDW correlations with a Luttinger exponent $K_{sc}\sim 1$ and strong divergent superconducting susceptibility. Our results provide a promising way to realize long-range superconductivity in the doped Hubbard model in two dimensions. The relevance of our numerical results to cuprate materials is also discussed.

Finite-temperature tensor network study of the Hubbard model on an infinite square lattice

The Hubbard model is a longstanding problem in the theory of strongly correlated electrons and a very active one in the experiments with ultracold fermionic atoms. Motivated by current and prospective quantum simulations, we apply a two-dimensional tensor network, an infinite projected entangled pair state, evolved in imaginary time by the neighborhood tensor update algorithm working directly in the thermodynamic limit. With U(1)xU(1) symmetry and the bond dimensions up to 29, we generate thermal states down to the temperature of 0.17 times the hopping rate. We obtain results for spin and charge correlators, unaffected by boundary effects. The spin correlators, measurable in prospective ultracold atoms experiments attempting to approach the thermodynamic limit, provide evidence of disruption of the antiferromagnetic background with mobile holes in a slightly doped Hubbard model. The charge correlators reveal the presence of hole-doublon pairs near half filling and signatures of hole-hole repulsion on doping. We also obtain specific heat in the slightly doped regime.

Quantifying the role of antiferromagnetic fluctuations in the superconductivity of the doped Hubbard model

We study the contribution of the electron-spin fluctuation coupling to the superconducting state of the two dimensional Hubbard model within dynamical cluster approximation (DCA) using a numerical exact continuous time Monte Carlo solver. By analyzing the frequency dependence of the self energy, we show that only about half of the superconductivity can be attributed to a "pairing glue" arising from treating spin fluctuations as a pairing boson in the standard one-loop theory.

Site-resolved observables in the doped spin-imbalanced triangular Hubbard model

The suppression of antiferromagnetic ordering in geometrically frustrated Hubbard models leads to a variety of exotic quantum phases including quantum spin liquids and chiral states. Here, we focus on the Hubbard model on one of the simplest frustrated lattice geometries, a triangular lattice. Motivated by the recent realization of ultracold fermionic atoms in triangular optical lattices, we study the properties of the triangular-lattice Hubbard model through a Numerical Linked-Cluster Expansion algorithm. We investigate the Mott insulator transition finding a critical interaction $U_c/t = 7.0(2)$ and use spatial two- and three-point correlation functions to explore doped and imbalanced systems. Our results demonstrate that many interesting features occur at temperatures previously obtained for ultracold fermions in optical lattices and are accessible by upcoming experiments. Our calculations will be helpful for thermometry in ultracold atom quantum simulators and can guide experimental searches for exotic quantum phases in atomic triangular Hubbard quantum simulators.

Superconductivity from repulsion in the doped 2D electronic Hubbard model: an entanglement perspective

We employ the momentum space entanglement renormalization group (MERG) scheme developed in references (Mukherjee et al 2021 J. High Energy Phys. JHEP04(2021)148; Patra and Lal 2021 Phys. Rev. B 104 144514) for the study of various insulating, superconducting and normal phases of the doped and the undoped 2D Hubbard model on a square lattice found recently by us (Mukherjee and Lal 2020 New J. Phys. 22 063007; Mukherjee and Lal 2020 New J. Phys. 22 063008). At each MERG step, disentanglement of particular degrees of freedom, transforms the tensor network representation of the many-particle states. The MERG reveals distinct holographic entanglement features for the normal metallic, topologically ordered insulating quantum liquid and Neél antiferromagnetic symmetry-broken ground states of the 2D Hubbard model at half-filling, clarifying the essence of the entanglement phase transitions that separates the three phases. An MERG analysis of the quantum critical point of the hole-doped 2D Hubbard model reveals the evolution of the many-particle entanglement of the quantum liquid ground state with hole-doping, as well as how the collapse of Mottness is concomitant with the emergence of d-wave superconductivity.

Effect of hole doping on the 120 degree order in the triangular lattice Hubbard model: a Hartree–Fock revisit

We revisit the unrestricted Hartree Fock study on the evolution of the ground state of the Hubbard model on the triangular lattice with hole doping. At half-filling, it is known that the ground state of the Hubbard model on triangular lattice develops a 120 degree coplanar order at half-filling in the strong interaction limit, i.e., in the spin 1/2 anti-ferromagnetic Heisenberg model on the triangular lattice. The ground state property in the doped case is still in controversy even though extensive studies were performed in the past. Within Hartree Fock theory, we find that the 120 degree order persists from zero doping to about 0.3 hole doping. At 1/3 hole doping, a three-sublattice collinear order emerges in which the doped hole is concentrated on one of the three sublattices with antiferromagnetic Neel order on the remaining two sublattices, which forms a honeycomb lattice. Between the 120 degree order and 1/3 doping region, a phase separation occurs in which the 120 degree order coexists with the collinear anti-ferromagnetic order in different regions of the system. The collinear phase extends from 1/3 doping to about 0.41 doping, beyond which the ground state is paramagnetic with uniform electron density. The phase diagram from Hartree Fock could provide guidance for the future study of the doped Hubbard model on triangular lattice with more sophisticated many-body approaches.

Stripes and spin-density waves in the doped two-dimensional Hubbard model: Ground state phase diagram

We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near $1/8$ doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about $1/5$ doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.

On the Nature of Valence Charge and Spin Excitations via Multi-Orbital Hubbard Models for Infinite-Layer Nickelates

Building upon the recent progress on the intriguing underlying physics for the newly discovered infinite-layer nickelates, in this article we review an examination of valence charge and spin excitations via multi-orbital Hubbard models as way to determine the fundamental building blocks for Hamiltonians that can describe the low energy properties of infinite-layer nickelates. We summarize key results from density-functional approaches, and apply them to the study of x-ray absorption to determine the valence ground states of infinite-layer nickelates in their parent form, and show that a fundamental d9 configuration as in the cuprates is incompatible with a self-doped ground state having holes in both dx2−y2 and a rare-earth-derived axial orbital. When doped, we determine that the rare-earth-derived orbitals empty and additional holes form low spin (S = 0) d8 Ni states, which can be well-described as a doped single-band Hubbard model. Using exact diagonalization for a 2-orbital model involving Ni and rare-earth orbitals, we find clear magnons at 1/2 filling that persist when doped, albeit with larger damping, and with a dependence on the precise orbital energy separation between the Ni- and rare-earth-derived orbitals. Taken together, a full two-band model for infinite-layer nickelates can well describe the valence charge and spin excitations observed experimentally.

Effective spin-charge transport theory and spin-transfer physics in frustrated magnets within the slave-boson approach

We study the electron dynamics in magnetic conductors with frustrated interactions dominated by isotropic exchange. We present a transport theory for itinerant carriers built upon the (single-band) doped Hubbard model and the slave-boson formalism, which incorporates the spin-exchange with the magnetically frustrated background into the representation of electron operators in a clear and controllable way. We also formulate hydrodynamic equations for the itinerant charge and spin degrees of freedom, whose currents contain new contributions that depend on the spatiotemporal variations of the order parameter of the frustrated magnet, which are described by Yang-Mills fields. Furthermore, we elucidate the transfer of angular momentum from the itinerant charge fluid to the magnet (i.e., the spin-transfer torque) via reciprocity arguments. A detailed microscopic derivation of our effective theory is also provided for one of the simplest models of frustrated magnetism, namely the Heisenberg antiferromagnet on a triangular lattice. Our findings point towards the possibility of previously unanticipated Hall physics in these frustrated platforms.