This thesis describes experiments performed on a gas of lithium atoms in an optical lattice which are relevant to the low-temperature physics of the Hubbard model. The Hubbard model is perhaps the most intensively studied model of fermionic particles on a lattice, both because it is the simplest such model which seriously challenges modern theoretical tools, and because it is a candidate model for describing high-temperature superconductivity in cuprate materials. The challenges associated with numerical simulations of the Hubbard model have motivated “quantum simulations” of the type described in this thesis, in which a well-controlled physical system is engineered to mimic the model, allowing its properties to be determined through experiments. Historically, however, such experiments could not access the low-temperature regime of the Hubbard model most relevant to open questions on materials like cuprates, where modern computational techniques cannot yet definitively characterize the physics. We developed a novel set of cooling techniques which allowed our experiment to access this regime. These techniques are based on adiabatic conversion between a band insulator, which can be prepared with very low entropies using single-atom-resolved control and detection techniques now common in ultracold atomic physics, and various quantum fluids, including Hubbard systems. We then performed the first experiments in this new regime of the Hubbard model, focusing on the anomalous “pseudogapped” metallic state. This state is predicted to exist by computational methods, and is believed to mimic the eponymous state in cuprates which gives rise to high-temperature superconductivity upon cooling. Nevertheless, its properties and mechanisms remain poorly understood. We developed several new measurements that enabled an observation of the pseudogap and which yielded novel data on its thermodynamics, linking the pseudogap to an anomaly in the equation of state. These experiments demonstrate a new mode of operating fermionic quantum gas experiments, in which single-atom-resolved controls, rather than bulk techniques like evaporative cooling, are leveraged to access low temperatures in a sample of several hundred atoms. They also signal the arrival of quantum simulators in a regime of poorly-understood emergent low-temperature phenomena that are relevant to the electronic properties of solids, concretely demonstrating the utility of quantum simulation in addressing frontier problems in quantum many-body physics.
This thesis describes experiments performed on a gas of lithium atoms in an optical lattice which are relevant to the low-temperature physics of the Hubbard model. The Hubbard model is perhaps the most intensively studied model of fermionic particles on a lattice, both because it is the simplest such model which seriously challenges modern theoretical tools, and because it is a candidate model for describing high-temperature superconductivity in cuprate materials. The challenges associated with numerical simulations of the Hubbard model have motivated “quantum simulations” of the type described in this thesis, in which a well-controlled physical system is engineered to mimic the model, allowing its properties to be determined through experiments. Historically, however, such experiments could not access the low-temperature regime of the Hubbard model most relevant to open questions on materials like cuprates, where modern computational techniques cannot yet definitively characterize the physics. We developed a novel set of cooling techniques which allowed our experiment to access this regime. These techniques are based on adiabatic conversion between a band insulator, which can be prepared with very low entropies using single-atom-resolved control and detection techniques now common in ultracold atomic physics, and various quantum fluids, including Hubbard systems. We then performed the first experiments in this new regime of the Hubbard model, focusing on the anomalous “pseudogapped” metallic state. This state is predicted to exist by computational methods, and is believed to mimic the eponymous state in cuprates which gives rise to high-temperature superconductivity upon cooling. Nevertheless, its properties and mechanisms remain poorly understood. We developed several new measurements that enabled an observation of the pseudogap and which yielded novel data on its thermodynamics, linking the pseudogap to an anomaly in the equation of state. These experiments demonstrate a new mode of operating fermionic quantum gas experiments, in which single-atom-resolved controls, rather than bulk techniques like evaporative cooling, are leveraged to access low temperatures in a sample of several hundred atoms. They also signal the arrival of quantum simulators in a regime of poorly-understood emergent low-temperature phenomena that are relevant to the electronic properties of solids, concretely demonstrating the utility of quantum simulation in addressing frontier problems in quantum many-body physics.
Studying the low-temperature phases of doped Hubbard models, including the emergence of collective order and exotic normal phases, is at the heart of strongly correlated physics. This thesis offers new theoretical insights into the low-energy physics of doped Hubbard models, with a particular emphasis on leveraging quantum simulation as a powerful investigative tool. We explore the emergence of collective order driven by many-body interactions, by predicting and observing stripe-like structures in quantum gas microscopes of engineered Hamiltonians. These fluctuating stripe patterns are developed into a theoretical framework for the pseudogap— a highly enigmatic phase in hole-doped cuprates—where we show how fluctuating domain walls yield effective toric code descriptions. Ultracold atoms in optical lattices provide a unique platform for directly probing and testing this theory. We further investigate unconventional bilayer nickelate superconductors, predicting an exceptionally stable superfluid state in a single-band Hubbard model at experimentally accessible temperatures. By proposing schemes to observe coherent pair-pair correlations, this advances the long-standing goal of realizing and observing long-range superconducting order in ultracold atomic systems in optical lattices. On non-bipartite lattices, we examine kinetic magnetism in Hubbard models and moiré heterostructures. Additionally, the low-temperature regime of doped Hubbard models with enhanced symmetries is analyzed, revealing exotic phenomena such as sub-dimensional polaronic particles. We develop machine learning techniques to extract key physical insights from many-body snapshots, offering a new avenue for understanding intricate quantum phases. Finally, we go beyond the Fermi-Hubbard paradigm and show how non-local strongly correlated models can be leveraged to address classical optimization problems through quantum annealing.
Among the most compelling targets for quantum simulation are strongly correlated phases and electronic-structure problems, where interactions organize electrons into collective states and produce intricate low-temperature phase diagrams. Cuprate superconductors are a canonical example: at low doping, superconductivity emerges from a metallic regime that exhibits a pseudogap above the critical temperature. Ultracold atoms in optical lattices provide a well-controlled, tunable setting to probe the microscopic mechanisms behind this behavior and, in fermion-native architectures, to implement electronic structure dynamics directly. This thesis demonstrates both large scale systematic analog quantum simulations in the most interesting regime of doped antiferromagnets and elementary high-fidelity gates as the basis for a digital or hybrid fermionic quantum processors in an optical lattice. Working with ultracold Lithium-6 in a two–dimensional optical superlattice, we first map the parameter space of the Fermi–Hubbard model over a wide range of temperatures and dopings. More than 40000 single-site spin- and density-resolved snapshots, spanning over a dense grid of doping levels and temperatures in a numerically challenging parameter space, reveal antiferromagnetic correlations extending across the system size and into the pseudogap regime. Temperatures down to T/t=0.2 are realized, determined by benchmarking the data to state-of-the-art numerical methods. By evaluating spin–charge correlators up to fifth order, we experimentally uncover a universal scaling behavior across the pseudogap crossover, governed by a single, doping-dependent temperature scale. These findings provide the first microscopic evidence of universal scaling in multi-point correlators and establish stringent benchmarks for future theoretical developments. On the same apparatus, we demonstrate that the controlled motion of fermions in double wells enables high-fidelity entanglement. Double-well potentials are used to engineer coherent pair-tunnelling and spin-exchange interactions that act as native SWAP^α operations. We demonstrate the collisional gate fidelities up to 99.8% and Bell-state lifetimes exceeding ten seconds, limited mainly by residual inhomogeneities in the optical potential. A composite pulse sequence isolates the pair-exchange channel, realizing the elementary building block required for fermionic quantum-chemistry algorithms. With the addition of site-selective addressing, our results mark a key step towards a fully digital, scalable fermionic quantum computer based on the controlled motion and entanglement of fermionic neutral atoms. The work establishes a hybrid analog-digital architecture: analog evolution under the Hubbard Hamiltonian yields thermodynamic states whose magnetic properties match and extend theoretical benchmarks, while digital collisional gates provide building blocks for error-corrected quantum computation and complex state preparation.
Ultracold atoms in optical lattices can provide faithful analog realizations of the Hubbard model and many of its variants. Recently, these simulators have been augmented with improved programmable controls that enable new routes to state preparation, new kinds of dynamics, and, by extension, access to new observables. I will discuss two examples where such tools have enabled simulations that challenge classical numerical techniques: first, a sampling problem involving non-interacting bosons on a lattice that is closely related to boson sampling; and second, the preparation of low temperature equilibrium states of interacting spinful fermions on a lattice. These simulations have the potential to enhance our understanding of quantum many body physics, including practical problems of interest in material design relating to strongly correlated quantum materials. However, fully harnessing their capabilities in the regime where exact classical numerical simulations are absent requires new, precision probes of the relevant microscopic or thermodynamic quantities.
An outstanding problem in physics is how to understand strongly interacting quantum many-body systems such as the quark-gluon plasma, neutron stars, superfluid 4He, and the high-temperature superconducting cuprates. The physics approach to this problem is to reduce these complex systems to minimal models that are believed to retain relevant phenomenology. For example, the Hubbard model — the focus of this thesis — describes quantum particles tunneling between sites of a lattice with on-site interactions. The Hubbard model is conjectured to describe the low-energy charge and spin properties of high-temperature superconducting cuprates. Thus far, there are no analytic solutions to the Hubbard model, and numerical calculations are difficult and even impossible in some regimes (e.g., the Fermi-Hubbard model away from half-filling). Therefore, whether the Hubbard model is a minimal model for the cuprates remains unresolved. In the face of these difficulties, a new approach has emerged — quantum simulation. The premise of quan-tum simulation is to perform experiments on a quantum system that is well-described by the model we are trying to study, has tunable parameters, and is easily probed. Ultracold atoms trapped in optical lattices are an ideal candidate for quantum simulation of the Hubbard models. This thesis describes work on two such systems: a 87Rb (boson) optical lattice experiment in the group of Brian DeMarco at the University
Over the past few decades, quantum simulation with ultracold atoms in optical lattices has proven to be a useful tool for understanding many-body Hamiltonians.One of its primary directions is the experimental study of the Fermi Hubbard Model (FHM), which aims to capture the fundamental properties of cuprate superconductors and has produced a wide range of intriguing phenomena, including the metal-to-Mott insulator transition, antiferromagnetism, strange metallicity, pseudogap physics, and is frequently used in the study of superconductivity. In recent years, new kinds of quantum simulators have emerged, and in particular we focus on two of them: 1) quantum simulators that study high symmetry fermions in optical lattices, and 2) optical tweezer arrays which provide a bottom-up approach to engineering the FHM. The first class of emerging quantum simulators utilizes alkaline-earth-like atoms (AEA’s) – which have large nuclear spin I (such as 173Yb and 87Sr) – to engineer the SU(N = 2I +1) symmetric FHM, where N = 2,..., 10. The second class corresponds to quantum simulators based on optical tweezer arrays which combine deterministic single tweezer ground state loading and tunneling coupled sites. Both classes of quantum simulators are interesting on their own. For example, the SU(N) FHM is generating a lot of interest due to its connections to multiorbital solid state systems and theoretical predictions of rich phase diagrams which exhibit exotic phases of matter, including unusual ones like spin liquids. On the other hand, optical tweezer arrays have gained a lot of attention since they can be used to study arbitrary geometries and access very low-entropy samples. This thesis reports on the numerical study of both quantum simulation architectures. Results of the SU(N) FHM at finite temperature, for different fillings and interaction strengths that span the non-interacting to the strongly interacting limit are obtained using Determinant Quantum Monte Carlo, Exact Diagonalization, and Numerical Linked Cluster Expansion. In particular, we focus on how thermodynamic and magnetic observables, such as the number of on-site pairs, energy, entropy, spin correlations, and structure factors depend on N, the interaction strength U/t, temperature T/t, and the anisotropy of tunneling rates. Our results demonstrate that in a homogeneous square lattice with one particle per site on average, thermodynamic observables as a function of temperature obey a universal scaling with N, and that short-range antiferromagnetic correlations are stronger for larger N and in lower dimensions. In addition, we present theory-experiment comparisons where possible, where we perform thermometry and provide a precise characterization of the equation of state of the SU(N) FHM. Our results emerge as a tool to perform thermometry in experiments from theory-experiment comparison and also provide guidance for future experiments with AEAs in optical lattices. Additionally, we present calculations using numerically precise discrete variable representation methods for two-dimensional stroboscopic tweezer arrays and compare the outcomes with experimental data. We quantify the effects of stroboscopic potentials on Hubbard parameters like the interaction strength U and the tunneling t in optical tweezer arrays and illustrate how heating from the stroboscopic potential relies on strobe frequency. Our calculations enable evaluation and optimization of two-dimensional tweezer array experiments. An important milestone of the results of this thesis corresponds to the record temperatures and entropies achieved with these quantum simulation architectures: 1) the coldest fermions ever created in nature in absolute temperature and in cold atoms, and 2) the lowest entropy per particle fermions ever created, albeit in a so-far small system.
Ultracold fermionic atoms in optical lattices offer pristine realizations of Hubbard models1, which are fundamental to modern condensed-matter physics2,3. Despite notable advancements4–6, the accessible temperatures in these optical lattice material analogues are still too high to address many open problems7–10. Here we demonstrate a several-fold reduction in temperature6,11–13, 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 parameters14,15, which is extremely challenging to simulate classically10. At half-filling, the long-range antiferromagnetic order is close to saturation, leading to a temperature of T/t=0.05−0.05+0.06 based on comparisons with numerically exact simulations. Doped away from half-filling, it is exceedingly challenging to realize systematically accurate and predictive numerical simulations9. 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 simulations16–18. Compared with the cuprates2,19,20, the reported temperatures correspond to a reduction from far above to below room temperature, at which physics such as the pseudogap and stripe phases may be expected3,19,21–24. 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 physics8,10. A several-fold reduction in temperature is accomplished using a neutral-atom Hubbard quantum simulator by transforming a low-entropy product state into strongly correlated states of interest via dynamic control of the model parameters.