One of the most interesting and challenging problems in physics is understanding strongly correlated many-body systems, where strong interactions can yield many remarkable phenomena such as superfluidity in 4He, high-temperature superconductivity, etc. In order to attack these problems, we often need to reduce the complexity of the systems to simple models in hopes of getting better insights into the properties of the systems. The Hubbard model, the focus of this dissertation, is one of the most famous examples of such model, which describes a tunneling of electrons between nearest neighbor sites of a lattice with on-site interactions. This simple model is an important concept in condensed matter physics and provides rich understandings of electronic and magnetic properties of materials. Despite its simplicity, there is no general analytical solution to the Hubbard model beyond 1D.;The discovery of ultracold atoms and optical lattices opens up the possibility of emulating the Hubbard model in experiments. Optical lattices provide an ideal realization of the Hubbard model where relevant parameters can be tuned systematically. It makes theoretical studies of the Hubbard model increasingly attractive since a direct comparison between theoretical calculations and experimental results becomes more and more possible.;In this dissertation, the ground-state properties of the repulsive Hubbard model for weak to intermediate interaction strengths in two, three dimensions and their dimensional crossover are studied within the mean field theory. We show that the system exhibits unidirectional spin-density wave (SDW) order with antiferromagnetic correlations and a long wavelength modulation. The modulating wave is along the [0011-direction at low interaction strength U/t and along the [1111-direction at higher U/t. The evolution of the wavelength of the SDW is determined as a function of U/t, the density, and t⊥/t. With an analysis of the pairing of spins based on nesting and deformation of the Fermi surface, we discuss how these results can be rationalized and how a simple, predictive model can be constructed for the properties of the SDW states.
One of the most interesting and challenging problems in physics is understanding strongly correlated many-body systems, where strong interactions can yield many remarkable phenomena such as superfluidity in 4He, high-temperature superconductivity, etc. In order to attack these problems, we often need to reduce the complexity of the systems to simple models in hopes of getting better insights into the properties of the systems. The Hubbard model, the focus of this dissertation, is one of the most famous examples of such model, which describes a tunneling of electrons between nearest neighbor sites of a lattice with on-site interactions. This simple model is an important concept in condensed matter physics and provides rich understandings of electronic and magnetic properties of materials. Despite its simplicity, there is no general analytical solution to the Hubbard model beyond 1D.;The discovery of ultracold atoms and optical lattices opens up the possibility of emulating the Hubbard model in experiments. Optical lattices provide an ideal realization of the Hubbard model where relevant parameters can be tuned systematically. It makes theoretical studies of the Hubbard model increasingly attractive since a direct comparison between theoretical calculations and experimental results becomes more and more possible.;In this dissertation, the ground-state properties of the repulsive Hubbard model for weak to intermediate interaction strengths in two, three dimensions and their dimensional crossover are studied within the mean field theory. We show that the system exhibits unidirectional spin-density wave (SDW) order with antiferromagnetic correlations and a long wavelength modulation. The modulating wave is along the [0011-direction at low interaction strength U/t and along the [1111-direction at higher U/t. The evolution of the wavelength of the SDW is determined as a function of U/t, the density, and t⊥/t. With an analysis of the pairing of spins based on nesting and deformation of the Fermi surface, we discuss how these results can be rationalized and how a simple, predictive model can be constructed for the properties of the SDW states.
Systems of fermions described by the three-dimensional (3D) repulsive Hubbard model on a cubic lattice have recently attracted considerable attention due to their possible experimental realization via cold atoms in an optical lattice. Analytical and numerical results are limited away from half-filling. We study the ground state of the doped system from weak to intermediate interaction strengths within the generalized Hartree-Fock approximation. The exact solution to the self-consistent-field equations in the thermodynamic limit is obtained and the ground state is shown to exhibit antiferromagnetic order and incommensurate spin-density waves (SDW). At low interaction strengths, the SDW has unidirectional character with a leading wavevector along the [100]-direction, and the system is metallic. As the interaction increases, the system undergoes a simultaneous structural and metal-to-insulator transition to a unidirectional SDW state along the [111]-direction, with a different wavelength. We systematically determine the real and momentum space properties of these states. The crossover from 3D to two dimensions (2D) is then studied by varying the inter-plane hopping amplitude, which can be experimentally realized by tuning the distance between a stack of square-lattice layers. Detailed comparisons are made between the exact numerical results and predictions from the pairing model, a variational ansatz based on the pairing of spins in the vicinity of the Fermi surface. Most of the numerical results can be understood quantitatively from the ansatz, which provides a simple picture for the nature of the SDW states.
This dissertation summarizes my recent work regarding systems of strongly interacting fermionic atoms in optical lattices. This work addresses the combination of two experimental techniques that have been the subject of much recent research in ultracold atom physics. One is the use of optical lattices, which provide a means to realize diverse interaction configurations within a clean, controllable system. The other is the use of magnetically tunable Feshbach resonances to control the strength of the interatomic interaction. Together, these techniques offer the possibility of an experimental realization of many important model Hamiltonians of condensed matter physics, and may also lead to the discovery of new physics. Recent study of this system has shown that strong interactions near Feshbach resonance will lead to the population of multiple lattice bands, and that collisions between atoms on neighboring sites cannot be neglected. These effects lead to a complicated Hamiltonian, but one which can be simplified to an effective single-band model equivalent to the generalized Hubbard model (GHM), which is an extension of the Hubbard model that includes correlated hopping terms. My main results concern the study of this model. The strong correlations between the particles make it difficult to definitively determine the many body physics of the GHM. As a first approach to understanding the GHM in optical lattices, I focus mainly on cases where the problem is greatly simplified by allowing interactions among only small groups of lattice sites. This restriction can be implemented in experiments using an optical superlattice potential. Our results include a proposed scheme (based on double-well superlattices) to empirically verify that the GHM describes this system and to directly measure the various parameters of this model. Other results include exact solutions on four-site square plaquettes, which demonstrate that d-wave excitations can occur in the low-energy states. By using a superlattice to give an array of weakly coupled plaquettes, one can thus produce a d-wave superfluid state. This is of relevance to the study of high-Tc superconductors, although I note certain key differences between the sort of d-wave superfluid described here and that of the superconductors.
In this thesis, I present the experimental realisation of the two-dimensional Hubbard model with ultracold fermionic Potassium-40 atoms in optical lattices. By tuning the dimensionality, kinetic and interaction in the optical lattices, we perform an analogue quantum simulation to explore the phase diagram of the Hubbard model. \n \nThe first key result is the experimental observation of particle-hole symmetry, namely a phase mapping between repulsive and attractive interactions. We compare a density-ordered, Mott-insulating phase with repulsive interaction to a spin-ordered, preformed pair phase and found excellent agreement with the theoretical prediction. The precise control and excellent detection capability of our quantum gas apparatus allow us to validate the particle-hole symmetry, and utilise it to explore quantum phases with a novel approach. \n \nNext, we investigate the spin-ordering on the repulsive side. By implementing a novel scheme based on coherent manipulation of spin correlations, we probe the anti-ferromagnetic ordering in the low-temperature phase diagram. The momentum-resolved spin correlations permit the reconstruction of spatial correlators without site-resolved imaging fidelity. \n \nFinally, we probe the attractive side of the phase diagram with a focus in pairing phenomenon, in which we draw a close analogy with BCS-BEC crossovers present in high-temperature superconducting cuprates and trapped ultracold Fermi gases. From the pair correlation function derived from thermodynamics observables, we observe the competition between the effective Pauli repulsion in fermionic systems and the on-site attraction we implemented.
Strongly correlated materials exhibit some of the most complex and technologically useful phenomena in condensed matter. We use the world line and determinant quantum Monte Carlo techniques to explore the phases of the Hubbard Hamiltonian, a model describing strong interactions, in a number of different contexts. We first consider the spin, charge, and bond order correlations of the one-dimensional extended fermion Hubbard model in the presence of a coupling to the lattice. A static alternating lattice distortion leads to enhanced charge density wave correlations at the expense of antiferromagnetic order. When the lattice degrees of freedom are dynamic, we show that a similar effect occurs even though the charge asymmetry must arise spontaneously. Although the evolution of the total energy with lattice coupling is smooth, the individual components exhibit sharp crossovers at the phase boundaries. Finally, we observe a tendency for bond order in the region between the charge and spin density wave phases. Second, we examine mixtures of bosons and fermions in one-dimensional optical lattices. We evaluate the density profiles and bosonic visibility [special characters omitted] resolving the discrepancy between theory and experiment by identifying parameter regimes where [special characters omitted] is reduced and increased. We present a simple qualitative picture of the different response to the fermion admixture in terms of the superfluid and Mott-insulating domains before and after the fermions are included. Finally, we show that [special characters omitted] exhibits kinks which are tied to the domain evolution present in the pure bosonic case, and also additional structure arising from the formation of boson-fermion molecules, a prediction for future experiments. Third, we report large scale calculations of the effective bandwidth, momentum distribution, and magnetic correlations of the square lattice fermion Hubbard Hamiltonian. The sharp Fermi surface of the non-interacting limit is significantly broadened by the electronic correlations, but retains signatures of the approach to the edges of the first Brillouin zone as the density increases. Finite size scaling of simulations on large lattices allows us to extract the interaction dependence of the antiferromagnetic order parameter, exhibiting its evolution from weak-coupling to the strong-coupling Heisenberg limit. Our lattices provide improved resolution of the momentum distribution, allowing a more quantitative comparison with time-of-flight optical lattice experiments.
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.
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