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.

Experimental realization of an extended Fermi-Hubbard model using a 2D lattice of dopant-based quantum dots.

The Hubbard model is an essential tool for understanding many-body physics in condensed matter systems. Artificial lattices of dopants in silicon are a promising method for the analog quantum simulation of extended Fermi-Hubbard Hamiltonians in the strong interaction regime. However, complex atom-based device fabrication requirements have meant emulating a tunable two-dimensional Fermi-Hubbard Hamiltonian in silicon has not been achieved. Here, we fabricate 3 × 3 arrays of single/few-dopant quantum dots with finite disorder and demonstrate tuning of the electron ensemble using gates and probe the many-body states using quantum transport measurements. By controlling the lattice constants, we tune the hopping amplitude and long-range interactions and observe the finite-size analogue of a transition from metallic to Mott insulating behavior. We simulate thermally activated hopping and Hubbard band formation using increased temperatures. As atomically precise fabrication continues to improve, these results enable a new class of engineered artificial lattices to simulate interactive fermionic models.

Lanczos-boosted numerical linked-cluster expansion for quantum lattice models.

Numerical linked-cluster expansions allow one to calculate finite-temperature properties of quantum lattice models directly in the thermodynamic limit through exact solutions of small clusters. However, full diagonalization is often the limiting factor for these calculations. Here we show that a partial diagonalization of the largest clusters in the expansion using the Lanczos algorithm can be as useful as full diagonalization for the method while mitigating some of the time and memory issues. As test cases, we consider the frustrated Heisenberg model on the checkerboard lattice and the Fermi-Hubbard model on the square lattice. We find that our approach can surpass state of the art in performance in a parallel environment.

Machine learning in electronic-quantum-matter imaging experiments.

For centuries, the scientific discovery process has been based on systematic human observation and analysis of natural phenomena1. Today, however, automated instrumentation and large-scale data acquisition are generating datasets of such large volume and complexity as to defy conventional scientific methodology. Radically different scientific approaches are needed, and machine learning (ML) shows great promise for research fields such as materials science2-5. Given the success of ML in the analysis of synthetic data representing electronic quantum matter (EQM)6-16, the next challenge is to apply this approach to experimental data-for example, to the arrays of complex electronic-structure images17 obtained from atomic-scale visualization of EQM. Here we report the development and training of a suite of artificial neural networks (ANNs) designed to recognize different types of order hidden in such EQM image arrays. These ANNs are used to analyse an archive of experimentally derived EQM image arrays from carrier-doped copper oxide Mott insulators. In these noisy and complex data, the ANNs discover the existence of a lattice-commensurate, four-unit-cell periodic, translational-symmetry-breaking EQM state. Further, the ANNs determine that this state is unidirectional, revealing a coincident nematic EQM state. Strong-coupling theories of electronic liquid crystals18,19 are consistent with these observations.

Spin transport in a Mott insulator of ultracold fermions.

Strongly correlated materials are expected to feature unconventional transport properties, such that charge, spin, and heat conduction are potentially independent probes of the dynamics. In contrast to charge transport, the measurement of spin transport in such materials is highly challenging. We observed spin conduction and diffusion in a system of ultracold fermionic atoms that realizes the half-filled Fermi-Hubbard model. For strong interactions, spin diffusion is driven by super-exchange and doublon-hole-assisted tunneling, and strongly violates the quantum limit of charge diffusion. The technique developed in this work can be extended to finite doping, which can shed light on the complex interplay between spin and charge in the Hubbard model.

Unsupervised machine learning account of magnetic transitions in the Hubbard model.

We employ several unsupervised machine learning techniques, including autoencoders, random trees embedding, and t-distributed stochastic neighboring ensemble (t-SNE), to reduce the dimensionality of, and therefore classify, raw (auxiliary) spin configurations generated, through Monte Carlo simulations of small clusters, for the Ising and Fermi-Hubbard models at finite temperatures. Results from a convolutional autoencoder for the three-dimensional Ising model can be shown to produce the magnetization and the susceptibility as a function of temperature with a high degree of accuracy. Quantum fluctuations distort this picture and prevent us from making such connections between the output of the autoencoder and physical observables for the Hubbard model. However, we are able to define an indicator based on the output of the t-SNE algorithm that shows a near perfect agreement with the antiferromagnetic structure factor of the model in two and three spatial dimensions in the weak-coupling regime. t-SNE also predicts a transition to the canted antiferromagnetic phase for the three-dimensional model when a strong magnetic field is present. We show that these techniques cannot be expected to work away from half filling when the "sign problem" in quantum Monte Carlo simulations is present.

Spin-imbalance in a 2D Fermi-Hubbard system.

The interplay of strong interactions and magnetic fields gives rise to unusual forms of superconductivity and magnetism in quantum many-body systems. Here, we present an experimental study of the two-dimensional Fermi-Hubbard model-a paradigm for strongly correlated fermions on a lattice-in the presence of a Zeeman field and varying doping. Using site-resolved measurements, we revealed anisotropic antiferromagnetic correlations, a precursor to long-range canted order. We observed nonmonotonic behavior of the local polarization with doping for strong interactions, which we attribute to the evolution from an antiferromagnetic insulator to a metallic phase. Our results pave the way to experimentally mapping the low-temperature phase diagram of the Fermi-Hubbard model as a function of both doping and spin polarization, for which many open questions remain.

Observation of spatial charge and spin correlations in the 2D Fermi-Hubbard model.

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.

Cooling Atomic Gases With Disorder.

Cold atomic gases have proven capable of emulating a number of fundamental condensed matter phenomena including Bose-Einstein condensation, the Mott transition, Fulde-Ferrell-Larkin-Ovchinnikov pairing, and the quantum Hall effect. Cooling to a low enough temperature to explore magnetism and exotic superconductivity in lattices of fermionic atoms remains a challenge. We propose a method to produce a low temperature gas by preparing it in a disordered potential and following a constant entropy trajectory to deliver the gas into a nondisordered state which exhibits these incompletely understood phases. We show, using quantum Monte Carlo simulations, that we can approach the Néel temperature of the three-dimensional Hubbard model for experimentally achievable parameters. Recent experimental estimates suggest the randomness required lies in a regime where atom transport and equilibration are still robust.

Compressibility of a fermionic mott insulator of ultracold atoms.

We characterize the Mott insulating regime of a repulsively interacting Fermi gas of ultracold atoms in a three-dimensional optical lattice. We use in situ imaging to extract the central density of the gas and to determine its local compressibility. For intermediate to strong interactions, we observe the emergence of a plateau in the density as a function of atom number, and a reduction of the compressibility at a density of one atom per site, indicating the formation of a Mott insulator. Comparisons to state-of-the-art numerical simulations of the Hubbard model over a wide range of interactions reveal that the temperature of the gas is of the order of, or below, the tunneling energy scale. Our results hold great promise for the exploration of many-body phenomena with ultracold atoms, where the local compressibility can be a useful tool to detect signatures of different phases or phase boundaries at specific values of the filling.

Observation of antiferromagnetic correlations in the Hubbard model with ultracold atoms.

Ultracold atoms in optical lattices have great potential to contribute to a better understanding of some of the most important issues in many-body physics, such as high-temperature superconductivity. The Hubbard model--a simplified representation of fermions moving on a periodic lattice--is thought to describe the essential details of copper oxide superconductivity. This model describes many of the features shared by the copper oxides, including an interaction-driven Mott insulating state and an antiferromagnetic (AFM) state. Optical lattices filled with a two-spin-component Fermi gas of ultracold atoms can faithfully realize the Hubbard model with readily tunable parameters, and thus provide a platform for the systematic exploration of its phase diagram. Realization of strongly correlated phases, however, has been hindered by the need to cool the atoms to temperatures as low as the magnetic exchange energy, and also by the lack of reliable thermometry. Here we demonstrate spin-sensitive Bragg scattering of light to measure AFM spin correlations in a realization of the three-dimensional Hubbard model at temperatures down to 1.4 times that of the AFM phase transition. This temperature regime is beyond the range of validity of a simple high-temperature series expansion, which brings our experiment close to the limit of the capabilities of current numerical techniques, particularly at metallic densities. We reach these low temperatures using a compensated optical lattice technique, in which the confinement of each lattice beam is compensated by a blue-detuned laser beam. The temperature of the atoms in the lattice is deduced by comparing the light scattering to determinant quantum Monte Carlo simulations and numerical linked-cluster expansion calculations. Further refinement of the compensated lattice may produce even lower temperatures which, along with light scattering thermometry, would open avenues for producing and characterizing other novel quantum states of matter, such as the pseudogap regime and correlated metallic states of the two-dimensional Hubbard model.

Magnetic correlations and pairing in the 1/5-depleted square lattice Hubbard model.

We study the single-orbital Hubbard model on the 1/5-depleted square-lattice geometry, which arises in such diverse systems as the spin-gap magnetic insulator CaV4O9 and ordered-vacancy iron selenides, presenting new issues regarding the origin of both magnetic ordering and superconductivity in these materials. We find a rich phase diagram that includes a plaquette singlet phase, a dimer singlet phase, a Néel and a block-spin antiferromagnetic phase, and stripe phases. Quantum Monte Carlo simulations show that the dominant pairing correlations at half filling change character from d wave in the plaquette phase to extended s wave upon transition to the Néel phase. These findings have intriguing connections to iron-based superconductors, and suggest that some physics of multiorbital systems can be captured by a single-orbital model at different dopings.

Linked-cluster expansion for the Green's function of the infinite-U Hubbard model.

We implement a highly efficient strong-coupling expansion for the Green's function of the Hubbard model. In the limit of extreme correlations, where the onsite interaction is infinite, the evaluation of diagrams simplifies dramatically enabling us to carry out the expansion to the eighth order in powers of the hopping amplitude. We compute the finite-temperature Green's function analytically in the momentum and Matsubara frequency space as a function of the electron density. Employing Padé approximations, we study the equation of state, Kelvin thermopower, momentum distribution function, quasiparticle fraction, and quasiparticle lifetime of the system at temperatures lower than, or of the order of, the hopping amplitude. We also discuss several different approaches for obtaining the spectral functions through analytic continuation of the imaginary frequency Green's function, and show results for the system near half filling. We benchmark our results for the equation of state against those obtained from a numerical linked-cluster expansion carried out to the eleventh order.

Fluctuation-dissipation theorem in an isolated system of quantum dipolar bosons after a quench.

We examine the validity of fluctuation-dissipation relations in isolated quantum systems taken out of equilibrium by a sudden quench. We focus on the dynamics of trapped hard-core bosons in one-dimensional lattices with dipolar interactions whose strength is changed during the quench. We find indications that fluctuation-dissipation relations hold if the system is nonintegrable after the quench, as well as if it is integrable after the quench if the initial state is an equilibrium state of a nonintegrable Hamiltonian. On the other hand, we find indications that they fail if the system is integrable both before and after quenching.

Short-range correlations and cooling of ultracold fermions in the honeycomb lattice.

We use determinantal quantum Monte Carlo simulations and numerical linked-cluster expansions to study thermodynamic properties and short-range spin correlations of fermions in the honeycomb lattice. We find that, at half filling and finite temperatures, nearest-neighbor spin correlations can be stronger in this lattice than in the square lattice, even in regimes where the ground state in the former is a semimetal or a spin liquid. The honeycomb lattice also exhibits a more pronounced anomalous region in the double occupancy that leads to stronger adiabatic cooling than in the square lattice. We discuss the implications of these findings for optical lattice experiments.

Quantum quenches in disordered systems: approach to thermal equilibrium without a typical relaxation time.

We study spectral properties and the dynamics after a quench of one-dimensional spinless fermions with short-range interactions and long-range random hopping. We show that a sufficiently fast decay of the hopping term promotes localization effects at finite temperature, which prevents thermalization even if the classical motion is chaotic. For slower decays, we find that thermalization does occur. However, within this model, the latter regime falls in an unexpected universality class, namely, observables exhibit a power-law (as opposed to an exponential) approach to their thermal expectation values.