Lithium Lab

Hubbard models can help us understand and design quantum materials – materials whose exceptional properties rely on the quantum mechanics of strongly interacting electrons. One prominent example of such materials are high temperature superconductors, which can conduct electricity without any losses and whose application may revolutionize power transmission, energy generation and storage, and medical imaging. Presently we find quantum materials by trial-and-error, and our understanding of how quantum materials work is limited. To better understand such materials, theoretical physicists crafted the Hubbard model, an idealized and simplified model that is believed to capture the rich quantum physics of such materials. It turns out, however, that this simplicity is highly deceiving. Even after sixty years of intensive study, understanding the properties of the Hubbard model remains extremely challenging even with today’s most advanced supercomputers.

The Hubbard model is a quantum model. What makes it so difficult to compute this quantum model on a classical supercomputer is exactly the complexity of quantum mechanics that is believed to make quantum computers much more powerful than classical computers. It is exactly this quantumness from which the surprising and useful material properties emerge.

The idea of our line of research is to build a quantum machine that can quantum mechanically instead of classically simulate the physics of the Hubbard model. We build specialized quantum devices, so called “quantum simulators”, that are purpose-built to exactly realize a wide range of Hubbard models. In recent work we achieved a major breakthrough, where we lowered the temperatures by about a factor of five. This, for the first time, brings quantum simulations into an entirely new regime: In this regime, classical supercomputers have an extremely difficult time providing answers. But it is also a regime where finding answers will have a profound impact on the understanding of quantum materials such as high temperature superconductors.

Now that we reach a low temperature regime that had been inaccessible for decades, we are excited to investigate the low temperature phases that emerge in the two-dimensional Hubbard model. Famously, d-wave superconducting and stripe phases are expected to emerge, and we are starting to explore the “pseudogap” regime that holds significant open questions. With our new, programmable optical lattice, we will be able to explore a wide range of geometries beyond the “plain vanilla” Hubbard model, and map out uncharted landscapes in which entirely new quantum phases may emerge.

Experimental setup. Ultracold atoms are trapped in a vacuum chamber beneath a high-resolution microscope

(photo by Greiner Lab)

Lithium lab in 2025

(photo by Greiner Lab)

A neutral-atom Hubbard quantum simulator in the cryogenic regime

06/2025

Nature (2025)

M. Xu, L. Kendrick, A. Kale, Y. Gang, Chunhan Feng, Shiwei Zhang, A. Young, M. Lebrat, M. Greiner

Ultracold fermionic atoms in optical lattices offer pristine realizations of Hubbard models, which are fundamental to modern condensed matter physics and the study of strongly-correlated quantum materials. Despite significant advancements, the accessible temperatures in these optical lattice material analogs are still too high to address many open problems beyond the reach of current numerical techniques. Here, we demonstrate a several-fold reduction in temperature, bringing large-scale quantum simulations of the Hubbard model into an entirely new regime. This is accomplished by transforming a low entropy product state into strongly-correlated states of interest via dynamic control of the model parameters, which is extremely challenging to simulate classically and so explored using the quantum simulator itself. At half filling, the long-range antiferromagnetic order is close to saturated, leading to a temperature of T/t=0.05 based on comparisons to numerically exact simulations. Doped away from half-filling no unbiased numerical simulation is available. Importantly, we are able to use quantum simulation to identify a new pathway for achieving similarly low temperatures with doping. This is confirmed by comparing short-range spin correlations to state-of-the-art, but approximate, constrained-path auxiliary field quantum Monte Carlo simulations. Compared to the cuprates, the reported temperatures correspond to a reduction from far above to significantly below room temperature, where physics such as the pseudogap and stripe phases may be expected. Our work opens the door to quantum simulations that solve open questions in material science, develop synergies with numerical methods and theoretical studies, and lead to discoveries of new physics.

Ferrimagnetism of ultracold fermions in a multi-band Hubbard system

04/2024

arXiv:2404.17555

M. Lebrat*, A. Kale*, L. Kendrick, M. Xu, Y. Gang, Alexander Nikolaenko, Pietro M. Bonetti, Subir Sachdev, M. Greiner

Strongly correlated materials feature multiple electronic orbitals which are crucial to accurately understand their many-body properties, from cuprate materials to twisted bilayer graphene. In such multi-band models, quantum interference can lead to dispersionless bands whose large degeneracy gives rise to itinerant magnetism even with weak interactions. Here, we report on signatures of a ferrimagnetic state realized in a Lieb lattice at half-filling, characterized by antialigned magnetic moments with antiferromagnetic correlations, concomitant with a finite spin polarization. We demonstrate their robustness when increasing repulsive interactions from the non-interacting to the Heisenberg regime, and study their emergence when continuously tuning the lattice unit cell from a square to a Lieb geometry. Our work paves the way towards exploring exotic phases in related multi-orbital models such as quantum spin liquids in kagome lattices and heavy fermion behavior in Kondo models.

Observation of Nagaoka Polarons in a Fermi-Hubbard Quantum Simulator

08/2023

Nature 629, 317–322 (2024)

M. Lebrat*, M. Xu*, L. Kendrick, A. Kale, Y. Gang, Pranav Seetharaman, Ivan Morera, Ehsan Khatami, Eugene Demler, M. Greiner

Quantum interference can deeply alter the nature of many-body phases of matter. In the paradigmatic case of the Hubbard model, Nagaoka famously proved that introducing a single itinerant charge can transform a paramagnetic insulator into a ferromagnet through path interference. However, a microscopic observation of such 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 lattice. Using quantum gas microscopy, we reveal these polarons as extended ferromagnetic bubbles around particle dopants arising from the local interplay of coherent dopant motion and spin exchange. In contrast, kinetic frustration due to the triangular geometry promotes antiferromagnetic polarons around hole dopants, as proposed by Haerter and Shastry. 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 simulation.

Frustration- and doping-induced magnetism in a Fermi-Hubbard simulator

12/2022

Nature 620, 971-976 (2023) arXiv:2212.13983

M. Xu, L. Kendrick, A. Kale, Y. Gang, G. Ji, Richard T. Scalettar, M. Lebrat, M. Greiner

Geometrical frustration in strongly correlated systems can give rise to a plethora of novel ordered states and intriguing magnetic phases such as quantum spin liquids. Promising candidate materials for such phases can be described by the Hubbard model on an anisotropic triangular lattice, a paradigmatic model capturing the interplay between strong correlations and magnetic frustration. However, the fate of frustrated magnetism in the presence of itinerant dopants remains unclear, as well as its connection to the doped phases of the square Hubbard model. Here, we probe the local spin order of a Hubbard model with controllable frustration and doping, using ultracold fermions in anisotropic optical lattices continuously tunable from a square to a triangular geometry. At half-filling and strong interactions U/t∼9, we observe at the single-site level how frustration reduces the range of magnetic correlations and drives a transition from a collinear Néel antiferromagnet to a short-range correlated 120∘ spiral phase. Away from half-filling, magnetic correlations show a pronounced asymmetry between particle and hole doping close to triangular geometries and hint at a transition to ferromagnetic order at a particle doping above 20%. This work paves the way towards exploring possible chiral ordered or superconducting phases in triangular lattices, and realizing t-tprime square lattice Hubbard models that may be essential to describe superconductivity in cuprate materials.

Schrieffer-Wolff transformations for experiments: Dynamically suppressing virtual doublon-hole excitations in a Fermi-Hubbard simulator

07/2022

Phys. Rev. A 106, 012428 (2022) arXiv:2203.07366

A. Kale, Jakob Hendrik Huhn, M. Xu, L. Kendrick, M. Lebrat, C. S. Chiu, G. Ji, Fabian Grusdt, Annabelle Bohrdt, M. Greiner

In strongly interacting systems with a separation of energy scales, low-energy effective Hamiltonians help provide insights into the relevant physics at low temperatures. The emergent interactions in the effective model are mediated by virtual excitations of high-energy states: For example, virtual doublon-hole excitations in the Fermi-Hubbard model mediate antiferromagnetic spin-exchange interactions in the derived effective model, known as the t−J−3s model. Formally this procedure is described by performing a unitary Schrieffer-Wolff basis transformation. In the context of quantum simulation, it can be advantageous to consider the effective model to interpret experimental results. However, virtual excitations such as doublon-hole pairs can obfuscate the measurement of physical observables. Here we show that quantum simulators allow one to access the effective model even more directly by performing measurements in a rotated basis. We propose a protocol to perform a Schrieffer-Wolff transformation on Fermi-Hubbard low-energy eigenstates (or thermal states) to dynamically prepare approximate t−J−3s model states using fermionic atoms in an optical lattice. Our protocol involves performing a linear ramp of the optical lattice depth, which is slow enough to eliminate the virtual doublon-hole fluctuations but fast enough to freeze out the dynamics in the effective model. We perform a numerical study using exact diagonalization and find an optimal ramp speed for which the state after the lattice ramp has maximal overlap with the t−J−3s model state. We compare our numerics to experimental data from our Lithium-6 fermionic quantum gas microscope and show a proof-of-principle demonstration of this protocol. More generally, this protocol can be beneficial to studies of effective models by enabling the suppression of virtual excitations in a wide range of quantum simulation experiments.