Microscopic evolution of doped Mott insulators from polaronic metal to Fermi liquid

Realizing Altermagnetism in Fermi-Hubbard Models with Ultracold Atoms

Altermagnetism represents a type of collinear magnetism, that is in some aspects distinct from ferromagnetism and from conventional antiferromagnetism. In contrast to the latter, sublattices of opposite spin are related by spatial rotations and not only by translations and inversions. As a result, altermagnets have spin-split bands leading to unique experimental signatures. Here, we show theoretically how a d-wave altermagnetic phase can be realized with ultracold fermionic atoms in optical lattices. We propose an altermagnetic Hubbard model with anisotropic next-nearest neighbor hopping and obtain the Hartree-Fock phase diagram. The altermagnetic phase separates in a metallic and an insulating phase and is robust over a large parameter regime. We show that one of the defining characteristics of altermagnetism, the anisotropic spin transport, can be probed with trap-expansion experiments.

Directly imaging spin polarons in a kinetically frustrated Hubbard system

The emergence of quasiparticles in quantum many-body systems underlies the rich phenomenology in many strongly interacting materials. In the context of doped Mott insulators, magnetic polarons are quasiparticles that usually arise from an interplay between the kinetic energy of doped charge carriers and superexchange spin interactions1-8. However, in kinetically frustrated lattices, itinerant spin polarons-bound states of a dopant and a spin flip-have been theoretically predicted even in the absence of superexchange coupling9-14. Despite their important role in the theory of kinetic magnetism, a microscopic observation of these polarons is lacking. Here we directly image itinerant spin polarons in a triangular-lattice Hubbard system realized with ultracold atoms, revealing enhanced antiferromagnetic correlations in the local environment of a hole dopant. In contrast, around a charge dopant, we find ferromagnetic correlations, a manifestation of the elusive Nagaoka effect15,16. We study the evolution of these correlations with interactions and doping, and use higher-order correlation functions to further elucidate the relative contributions of superexchange and kinetic mechanisms. The robustness of itinerant spin polarons at high temperature paves the way for exploring potential mechanisms for hole pairing and superconductivity in frustrated systems10,11. Furthermore, our work provides microscopic insights into related phenomena in triangular-lattice moiré materials17-20.

Observation of Nagaoka polarons in a Fermi-Hubbard quantum simulator

Quantum interference can deeply alter the nature of many-body phases of matter1. In the case of the Hubbard model, Nagaoka proved that introducing a single itinerant charge can transform a paramagnetic insulator into a ferromagnet through path interference2-4. However, a microscopic observation of this kinetic magnetism induced by individually imaged dopants has been so far elusive. Here we demonstrate the emergence of Nagaoka polarons in a Hubbard system realized with strongly interacting fermions in a triangular optical lattice5,6. Using quantum gas microscopy, we image these polarons as extended ferromagnetic bubbles around particle dopants arising from the local interplay of coherent dopant motion and spin exchange. By contrast, kinetic frustration due to the triangular geometry promotes antiferromagnetic polarons around hole dopants7. Our work augurs the exploration of exotic quantum phases driven by charge motion in strongly correlated systems and over sizes that are challenging for numerical simulation8-10.

Equation of State and Thermometry of the 2D SU(N) Fermi-Hubbard Model

We characterize the equation of state (EoS) of the SU(N>2) Fermi-Hubbard Model (FHM) in a two-dimensional single-layer square optical lattice. We probe the density and the site occupation probabilities as functions of interaction strength and temperature for N=3, 4, and 6. Our measurements are used as a benchmark for state-of-the-art numerical methods including determinantal quantum Monte Carlo and numerical linked cluster expansion. By probing the density fluctuations, we compare temperatures determined in a model-independent way by fitting measurements to numerically calculated EoS results, making this a particularly interesting new step in the exploration and characterization of the SU(N) FHM.

Interference of holon strings in 2D Hubbard model

The 2D Hubbard model with large repulsion is an important problem in condensed matter physics. At half filling, its ground state is an antiferromagnet (AMF). The dope AMF below half filling is believed to capture the physics of highTcsuperconductors. And the fermion excitation of this dope AMF is theorized as splitting up into holons and spinons that carry charge and spin separately. It is believed that these exotic holons and spinons are the origins of the unusual properties of highTcsuperconductors. Despite the interests in holons and spinons, the direct observations of these excitations remain difficult in solid state experiments. Here, we show that with the rapid advances in the experimental techniques in cold atoms, the direct observation of holons is possible in quantum quench dynamic processes in cold atom settings. We show that the well-known holon-strings generated by the motion of a holon as well as their interferences can be detected by the measurements spin-spin correlations and demonstrate the presence of the Marshall phase associated with a holon string reflecting an underlying AMF background. Moreover, we show that the interferences of the holon strings make a holon propagate anisotropically, with a diffusion pattern clearly distinct from that of spinless fermions. At the same time, we show that these interferences lead to a large suppression in magnetic order in the region swept through by the strings (even to about 95% for some bond). We further demonstrate the Marshall phase of the holon-strings by comparing the dynamics of holon in thetJmodel with that of the so-calledσtJ-model, which is thetJmodel with the Marshall phase removed. The holons in these models propagate entirely differently.

Dynamics of a Magnetic Polaron in an Antiferromagnet

The t-J model remains an indispensable construct in high-temperature superconductivity research, bridging the gap between charge dynamics and spin interactions within antiferromagnetic matrices. This study employs the multiple Davydov Ansatz method with thermo-field dynamics to dissect the zero-temperature and finite-temperature behaviors. We uncover the nuanced dependence of hole and spin deviation dynamics on the spin-spin coupling parameter J, revealing a thermally-activated landscape where hole mobilities and spin deviations exhibit a distinct temperature-dependent relationship. This numerically accurate thermal perspective augments our understanding of charge and spin dynamics in an antiferromagnet.

Finite-Temperature Hole-Magnon Dynamics in an Antiferromagnet

Employing the numerically accurate multiple Davydov Ansatz in combination with the thermo-field dynamics approach, we delve into the interplay of the finite-temperature dynamics of holes and magnons in an antiferromagnet, which allows for scrutinizing previous predictions from the self-consistent Born approximation while offering, for the first time, accurate finite-temperature computation of detailed magnon dynamics as a response and a facilitator to the hole motion. The study also uncovers a pronounced temperature dependence of the magnon and hole populations, pointing to the feasibility of potential thermal manipulation and control of hole dynamics. Our methodology can be applied not only to the calculation of steady-state angular-resolved photoemission spectra but also to the simulation of femtosecond terahertz pump-probe and other nonlinear signals for the characterization of antiferromagnetic materials.

Commensurate and incommensurate 1D interacting quantum systems

Single-atom imaging resolution of many-body quantum systems in optical lattices is routinely achieved with quantum-gas microscopes. Key to their great versatility as quantum simulators is the ability to use engineered light potentials at the microscopic level. Here, we employ dynamically varying microscopic light potentials in a quantum-gas microscope to study commensurate and incommensurate 1D systems of interacting bosonic Rb atoms. Such incommensurate systems are analogous to doped insulating states that exhibit atom transport and compressibility. Initially, a commensurate system with unit filling and fixed atom number is prepared between two potential barriers. We deterministically create an incommensurate system by dynamically changing the position of the barriers such that the number of available lattice sites is reduced while retaining the atom number. Our systems are characterised by measuring the distribution of particles and holes as a function of the lattice filling, and interaction strength, and we probe the particle mobility by applying a bias potential. Our work provides the foundation for preparation of low-entropy states with controlled filling in optical-lattice experiments.

Direct observation of nonlocal fermion pairing in an attractive Fermi-Hubbard gas

The Hubbard model of attractively interacting fermions provides a paradigmatic setting for fermion pairing. It features a crossover between Bose-Einstein condensation of tightly bound pairs and Bardeen-Cooper-Schrieffer superfluidity of long-range Cooper pairs, and a "pseudo-gap" region where pairs form above the superfluid critical temperature. We directly observe the nonlocal nature of fermion pairing in a Hubbard lattice gas, using spin- and density-resolved imaging of [Formula: see text]1000 fermionic potassium-40 atoms under a bilayer microscope. Complete fermion pairing is revealed by the vanishing of global spin fluctuations with increasing attraction. In the strongly correlated regime, the fermion pair size is found to be on the order of the average interparticle spacing. Our study informs theories of pseudo-gap behavior in strongly correlated fermion systems.

Witnessing light-driven entanglement using time-resolved resonant inelastic X-ray scattering

Characterizing and controlling entanglement in quantum materials is crucial for the development of next-generation quantum technologies. However, defining a quantifiable figure of merit for entanglement in macroscopic solids is theoretically and experimentally challenging. At equilibrium the presence of entanglement can be diagnosed by extracting entanglement witnesses from spectroscopic observables and a nonequilibrium extension of this method could lead to the discovery of novel dynamical phenomena. Here, we propose a systematic approach to quantify the time-dependent quantum Fisher information and entanglement depth of transient states of quantum materials with time-resolved resonant inelastic x-ray scattering. Using a quarter-filled extended Hubbard model as an example, we benchmark the efficiency of this approach and predict a light-enhanced many-body entanglement due to the proximity to a phase boundary. Our work sets the stage for experimentally witnessing and controlling entanglement in light-driven quantum materials via ultrafast spectroscopic measurements.

Tangent Space Approach for Thermal Tensor Network Simulations of the 2D Hubbard Model

Accurate simulations of the two-dimensional (2D) Hubbard model constitute one of the most challenging problems in condensed matter and quantum physics. Here we develop a tangent space tensor renormalization group (tanTRG) approach for the calculations of the 2D Hubbard model at finite temperature. An optimal evolution of the density operator is achieved in tanTRG with a mild O(D^{3}) complexity, where the bond dimension D controls the accuracy. With the tanTRG approach we boost the low-temperature calculations of large-scale 2D Hubbard systems on up to a width-8 cylinder and 10×10 square lattice. For the half-filled Hubbard model, the obtained results are in excellent agreement with those of determinant quantum Monte Carlo (DQMC). Moreover, tanTRG can be used to explore the low-temperature, finite-doping regime inaccessible for DQMC. The calculated charge compressibility and Matsubara Green's function are found to reflect the strange metal and pseudogap behaviors, respectively. The superconductive pairing susceptibility is computed down to a low temperature of approximately 1/24 of the hopping energy, where we find d-wave pairing responses are most significant near the optimal doping. Equipped with the tangent-space technique, tanTRG constitutes a well-controlled, highly efficient and accurate tensor network method for strongly correlated 2D lattice models at finite temperature.

Fermionic-propagator and alternating-basis quantum Monte Carlo methods for correlated electrons on a lattice

Ultracold-atom simulations of the Hubbard model provide insights into the character of charge and spin correlations in and out of equilibrium. The corresponding numerical simulations, on the other hand, remain a significant challenge. We build on recent progress in the quantum Monte Carlo (QMC) simulation of electrons in continuous space and apply similar ideas to the square-lattice Hubbard model. We devise and benchmark two discrete-time QMC methods, namely the fermionic-propagator QMC (FPQMC) and the alternating-basis QMC (ABQMC). In FPQMC, the time evolution is represented by snapshots in real space, whereas the snapshots in ABQMC alternate between real and reciprocal space. The methods may be applied to study equilibrium properties within the grand-canonical or canonical ensemble, external field quenches, and even the evolution of pure states. Various real-space/reciprocal-space correlation functions are also within their reach. Both methods deal with matrices of size equal to the number of particles (thus independent of the number of orbitals or time slices), which allows for cheap updates. We benchmark the methods in relevant setups. In equilibrium, the FPQMC method is found to have an excellent average sign and, in some cases, yields correct results even with poor imaginary-time discretization. ABQMC has a significantly worse average sign, but also produces good results. Out of equilibrium, FPQMC suffers from a strong dynamical sign problem. On the contrary, in ABQMC, the sign problem is not time-dependent. Using ABQMC, we compute survival probabilities for several experimentally relevant pure states.

Magnetically mediated hole pairing in fermionic ladders of ultracold atoms

Conventional superconductivity emerges from pairing of charge carriers-electrons or holes-mediated by phonons¹. In many unconventional superconductors, the pairing mechanism is conjectured to be mediated by magnetic correlations², as captured by models of mobile charges in doped antiferromagnets³. However, a precise understanding of the underlying mechanism in real materials is still lacking and has been driving experimental and theoretical research for the past 40 years. Early theoretical studies predicted magnetic-mediated pairing of dopants in ladder systems^4-8, in which idealized theoretical toy models explained how pairing can emerge despite repulsive interactions⁹. Here we experimentally observe this long-standing theoretical prediction, reporting hole pairing due to magnetic correlations in a quantum gas of ultracold atoms. By engineering doped antiferromagnetic ladders with mixed-dimensional couplings¹⁰, we suppress Pauli blocking of holes at short length scales. This results in a marked increase in binding energy and decrease in pair size, enabling us to observe pairs of holes predominantly occupying the same rung of the ladder. We find a hole-hole binding energy of the order of the superexchange energy and, upon increased doping, we observe spatial structures in the pair distribution, indicating repulsion between bound hole pairs. By engineering a configuration in which binding is strongly enhanced, we delineate a strategy to increase the critical temperature for superconductivity.

Nonequilibrium Hole Dynamics in Antiferromagnets: Damped Strings and Polarons

We develop a nonperturbative theory for hole dynamics in antiferromagnetic spin lattices, as described by the t-J model. This is achieved by generalizing the self-consistent Born approximation to nonequilibrium systems, making it possible to calculate the full time-dependent many-body wave function. Our approach reveals three distinct dynamical regimes, ultimately leading to the formation of magnetic polarons. Following the initial ballistic stage of the hole dynamics, coherent formation of string excitations gives rise to characteristic oscillations in the hole density. Their damping eventually leaves behind magnetic polarons that undergo ballistic motion with a greatly reduced velocity. The developed theory provides a rigorous framework for understanding nonequilibrium physics of defects in quantum magnets and quantitatively explains recent observations from cold-atom quantum simulations in the strong coupling regime.

Rotational Resonances and Regge-like Trajectories in Lightly Doped Antiferromagnets

Understanding the nature of charge carriers in doped Mott insulators holds the key to unravelling puzzling properties of strongly correlated electron systems, including cuprate superconductors. Several theoretical models suggested that dopants can be understood as bound states of partons, the analogues of quarks in high-energy physics. However, direct signatures of spinon-chargon bound states are lacking, both in experiment and theory. Here we propose a rotational variant of angle-resolved photo-emission spectroscopy (ARPES) and calculate rotational spectra numerically using the density-matrix renormalization group. We identify long-lived rotational resonances for an individual dopant, which we interpret as a direct indicator of the microscopic structure of spinon-chargon bound states. Similar to Regge trajectories reflecting the quark structure of mesons, we establish a linear dependence of the rotational energy on the superexchange coupling. The rotational peaks we find are strongly suppressed in standard ARPES spectra, but we suggest a multiphoton extension of ARPES which allows us to access rotational spectra. Our findings suggest that multiphoton spectroscopy experiments should provide new insights into emergent universal features of strongly correlated electron systems.

Microscopic evolution of doped Mott insulators from polaronic metal to Fermi liquid

[Figure: see text].