We employ a family of ancilla qubit variational wave functions [Y.-H. Zhang and S. Sachdev, Phys. Rev. Res. 2 , 023172 (2020)] to describe the polaronic correlations in the pseudogap metal phase of a hole-doped 2D Fermi–Hubbard model. Comparison to ultracold atom quantum simulator data [J. Koepsell et al ., Science 374 , 82 (2021)] reveals both qualitative and quantitative agreement with the numerical analysis from half-filling up to 80% hole-doping, capturing the crossover from the polaronic regime to the Fermi liquid observed around 40% doping.
The interplay between spin and charge degrees of freedom arising from doping a Mott insulating quantum spin liquid (QSL) has been a topic of research for several decades. Calculating properties of these fractionalized metallic states in single-band models are generally restricted to mean-field patron descriptions and small fluctuations around these states, which are insufficient for quantitative comparison of observables to measurements performed in strongly-correlated systems. In this work, we numerically study a class of correlated electronic wavefunctions which support fractionalized spin and charge excitations and which fully take into account gauge fluctuations through the enforcement of local Hilbert space constraints. By optimizing the energy of these wavefunctions against the hole-doped Fermi Hubbard Hamiltonian, we obtain a variational ansatz for describing the low-energy physics of this model. We compare measurements of hole-induced spin-spin correlation functions to measurements taken in low temperature cold-atom simulations of the Hubbard model and find quantitative agreement between the two. In particular, we demonstrate the emergence of magnetic polaron correlations in these metallic states.
The competition between antiferromagnetism and hole motion in two-dimensional Mott insulators lies at the heart of a doping-dependent transition from an anomalous metal to a conventional Fermi liquid. We observe such a crossover in Fermi-Hubbard systems on a cold-atom quantum simulator and reveal the transformation of multipoint correlations between spins and holes upon increasing doping at temperatures around the superexchange energy. Conventional observables, such as spin susceptibility, are furthermore computed from the microscopic snapshots of the system. Starting from a magnetic polaron regime, we find the system evolves into a Fermi liquid featuring incommensurate magnetic fluctuations and fundamentally altered correlations. The crossover is completed for hole dopings around 30%. Our work benchmarks theoretical approaches and discusses possible connections to lower-temperature phenomena.
Understanding quantum many-body states of correlated electrons is one main theme in modern condensed-matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, was recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Chen et al., Phys. Rev. X 8, 031082 (2018)], complemented by independent determinant quantum Monte Carlo, offers a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to an unprecedented low temperature of a few percent of the tunneling energy. We observe excellent agreement with ultracold fermion measurements at both half filling and finite doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behaviors of the antimoments.
Two of the iconic phases of the hole-doped cuprate materials are the intermediate temperature pseudogap metal and the lower temperature
d-wave superconductor. Following the suggestion of P. W. Anderson, there were early theories of these phases as doped quantum spin liquids. However, these theories have had difficulties with two prominent observations: (i) angle-dependent magnetoresistance measurements (ADMR) in the pseudogap metal, including observation of the Yamaji effect, present convincing evidence of small hole pockets which can tunnel coherently between square lattice layers (Fang et al., Nature Physics 18, 558 (2022); Chan et al., Nature Physics 21, 1753 (2025)) and (ii) the velocities of the nodal Bogoliubov quasiparticles in the
d-wave superconductor are highly anisotropic, with
vF≫vΔ (Chiao et al., Phys. Rev. B 62, 3554 (2000)). These notes review how the fractionalized Fermi Liquid (FL*) state, which dopes quantum spin liquids with gauge-neutral electron-like quasiparticles, resolves both difficulties. Theories of insulating quantum spin liquids employing fractionalization of the electron spin into bosonic or fermionic partons are discussed. Doping the bosonic parton theory leads to a holon metal theory: while not appropriate for the cuprate pseudogap, this theory is argued to apply to ultracold atom experiments on the Lieb lattice. Doping the fermionic parton theory leads to a
d-wave superconductor with nearly isotropic quasiparticle velocities. The construction of the FL* state in a single band Hubbard-type model is described using a quantum dimer model, followed by a more realistic description using the Ancilla Layer Model (ALM), which is then used to obtain the theory of the pseudogap and the
d-wave superconductor. The ALM also leads to a wavefunction for the FL* state of the Hubbard model, which agrees with ultracold atom observations.
Strong electron correlations lie at the origin of high-temperature superconductivity. Its essence is believed to be captured by the Fermi-Hubbard model of repulsively interacting fermions on a lattice. Here we report on the site-resolved observation of charge and spin correlations in the two-dimensional (2D) Fermi-Hubbard model realized with ultracold atoms. Antiferromagnetic spin correlations are maximal at half-filling and weaken monotonically upon doping. At large doping, nearest-neighbor correlations between singly charged sites are negative, revealing the formation of a correlation hole, the suppressed probability of finding two fermions near each other. As the doping is reduced, the correlations become positive, signaling strong bunching of doublons and holes, in agreement with numerical calculations. The dynamics of the doublon-hole correlations should play an important role for transport in the Fermi-Hubbard model.