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.

Cold-Atom Elevator: From Edge-State Injection to the Preparation of Fractional Chern Insulators

Optical box traps offer new possibilities for quantum-gas experiments. Building on their exquisite spatial and temporal control, we propose to engineer system-reservoir configurations using box traps, in view of preparing and manipulating topological atomic states in optical lattices. First, we consider the injection of particles from the reservoir to the system: this scenario is shown to be particularly well suited to activating energy-selective chiral edge currents, but also to prepare fractional Chern insulating ground states. Then, we devise a practical evaporative-cooling scheme to effectively cool down atomic gases into topological ground states. Our open-system approach to optical-lattice settings provides a new path for the investigation of ultracold quantum matter, including strongly correlated and topological phases.

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.

Accreditation of analogue quantum simulators

We present an accreditation protocol for analogue, i.e., continuous-time, quantum simulators. For a given simulation task, it provides an upper bound on the variation distance between the probability distributions at the output of an erroneous and error-free analogue quantum simulator. As its overheads are independent of the size and nature of the simulation, the protocol is ready for immediate usage and practical for the long term. It builds on the recent theoretical advances of strongly universal Hamiltonians and quantum accreditation as well as experimental progress toward the realization of programmable hybrid analogue-digital quantum simulators.

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.

The Wiedemann-Franz law in doped Mott insulators without quasiparticles

Many metallic quantum materials display anomalous transport phenomena that defy a Fermi liquid description. Here, we use numerical methods to calculate thermal and charge transport in the doped Hubbard model and observe a crossover separating high- and low-temperature behaviors. Distinct from the behavior at high temperatures, the Lorenz number [Formula: see text] becomes weakly doping dependent and less sensitive to parameters at low temperatures. At the lowest numerically accessible temperatures, [Formula: see text] roughly approaches the Wiedemann-Franz constant [Formula: see text], even in a doped Mott insulator that lacks well-defined quasiparticles. Decomposing the energy current operator indicates a compensation between kinetic and potential contributions, which may help to clarify the interpretation of transport experiments beyond Boltzmann theory in strongly correlated metals.

L-based numerical linked cluster expansion for square lattice models

We introduce a numerical linked cluster expansion for square-lattice models whose building block is an L-shape cluster. For the spin-1/2 models studied in this work, we find that this expansion exhibits a similar or better convergence of the bare sums than that of the (larger) square-shaped clusters and can be used with resummation techniques (like the site- and bond-based expansions) to obtain results at even lower temperatures. We compare the performance of weak- and strong-embedding versions of this expansion in various spin-1/2 models and show that the strong-embedding version is preferable because of its convergence properties and lower computational cost. Finally, we show that the expansion based on the L-shape cluster can be naturally used to study properties of lattice models that smoothly connect the square and triangular lattice geometries.

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.

Two-Dimensional Programmable Tweezer Arrays of Fermions

We prepare high-filling two-component arrays of tens of fermionic ^{6}Li atoms in optical tweezers, with the atoms in the ground motional state of each tweezer. Using a stroboscopic technique, we configure the arrays in various two-dimensional geometries with negligible Floquet heating. A full spin- and density-resolved readout of individual sites allows us to postselect near-zero entropy initial states for fermionic quantum simulation. We prepare a correlated state in a two-by-two tunnel-coupled Hubbard plaquette, demonstrating all the building blocks for realizing a programmable fermionic quantum simulator.

Evaluating Second-Order Phase Transitions with Diagrammatic Monte Carlo: Néel Transition in the Doped Three-Dimensional Hubbard Model

Diagrammatic Monte Carlo-the technique for the numerically exact summation of all Feynman diagrams to high orders-offers a unique unbiased probe of continuous phase transitions. Being formulated directly in the thermodynamic limit, the diagrammatic series is bound to diverge and is not resummable at the transition due to the nonanalyticity of physical observables. This enables the detection of the transition with controlled error bars from an analysis of the series coefficients alone, avoiding the challenge of evaluating physical observables near the transition. We demonstrate this technique by the example of the Néel transition in the 3D Hubbard model. At half filling and higher temperatures, the method matches the accuracy of state-of-the-art finite-size techniques, but surpasses it at low temperatures and allows us to map the phase diagram in the doped regime, where finite-size techniques struggle from the fermion sign problem. At low temperatures and sufficient doping, the transition to an incommensurate spin density wave state is observed.

Quantum gas microscopy of Kardar-Parisi-Zhang superdiffusion

The Kardar-Parisi-Zhang (KPZ) universality class describes the coarse-grained behavior of a wealth of classical stochastic models. Surprisingly, KPZ universality was recently conjectured to also describe spin transport in the one-dimensional quantum Heisenberg model. We tested this conjecture by experimentally probing transport in a cold-atom quantum simulator via the relaxation of domain walls in spin chains of up to 50 spins. We found that domain-wall relaxation is indeed governed by the KPZ dynamical exponent z = 3/2 and that the occurrence of KPZ scaling requires both integrability and a nonabelian SU(2) symmetry. Finally, we leveraged the single-spin-sensitive detection enabled by the quantum gas microscope to measure an observable based on spin-transport statistics. Our results yield a clear signature of the nonlinearity that is a hallmark of KPZ universality.

Approximating quantum thermodynamic properties using DFT

The fabrication, utilisation, and efficiency of quantum technology devices rely on a good understanding of quantum thermodynamic properties. Many-body systems are often used as hardware for these quantum devices, but interactions between particles make the complexity of related calculations grow exponentially with the system size. Here we explore and systematically compare 'simple' and 'hybrid' approximations to the average work and entropy variation built on static density functional theory concepts. These approximations are computationally cheap and could be applied to large systems. We exemplify them considering driven one-dimensional Hubbard chains and show that, for 'simple' approximations and low to medium temperatures, it pays to consider a good estimate of the Kohn-Sham Hamiltonian to approximate the driving Hamiltonian. Our results confirm that a 'hybrid' approach, requiring a very good approximation of the initial and, for the entropy, final states of the system, provides great improvements. This approach should be particularly efficient when many-body effects are not increased by the driving Hamiltonian.

Phonon-Induced Pairing in Quantum Dot Quantum Simulator

Quantum simulations can provide new insights into the physics of strongly correlated electronic systems. A well-studied system, but still open in many regards, is the Hubbard-Holstein Hamiltonian, where electronic repulsion is in competition with attraction generated by the electron-phonon coupling. In this context, we study the behavior of four quantum dots in a suspended carbon nanotube and coupled to its flexural degrees of freedom. The system is described by a Hamiltonian of the Hubbard-Holstein class, where electrons on different sites interact with the same phonon. We find that the system presents a transition from the Mott insulating state to a polaronic state, with the appearance of pairing correlations and the breaking of the translational symmetry. These findings will motivate further theoretical and experimental efforts to employ nanoelectromechanical systems to simulate strongly correlated systems with electron-phonon interactions.

Quantum Variational Learning of the Entanglement Hamiltonian

Learning the structure of the entanglement Hamiltonian (EH) is central to characterizing quantum many-body states in analog quantum simulation. We describe a protocol where spatial deformations of the many-body Hamiltonian, physically realized on the quantum device, serve as an efficient variational ansatz for a local EH. Optimal variational parameters are determined in a feedback loop, involving quench dynamics with the deformed Hamiltonian as a quantum processing step, and classical optimization. We simulate the protocol for the ground state of Fermi-Hubbard models in quasi-1D geometries, finding excellent agreement of the EH with Bisognano-Wichmann predictions. Subsequent on-device spectroscopy enables a direct measurement of the entanglement spectrum, which we illustrate for a Fermi Hubbard model in a topological phase.

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

[Figure: see text].

Recovering Quantum Correlations in Optical Lattices from Interaction Quenches

Quantum simulations with ultracold atoms in optical lattices open up an exciting path toward understanding strongly interacting quantum systems. Atom gas microscopes are crucial for this as they offer single-site density resolution, unparalleled in other quantum many-body systems. However, currently a direct measurement of local coherent currents is out of reach. In this Letter, we show how to achieve that by measuring densities that are altered in response to quenches to noninteracting dynamics, e.g., after tilting the optical lattice. For this, we establish a data analysis method solving the closed set of equations relating tunneling currents and atom number dynamics, allowing us to reliably recover the full covariance matrix, including off-diagonal terms representing coherent currents. The signal processing builds upon semidefinite optimization, providing bona fide covariance matrices optimally matching the observed data. We demonstrate how the obtained information about noncommuting observables allows one to quantify entanglement at finite temperature, which opens up the possibility to study quantum correlations in quantum simulations going beyond classical capabilities.

Real-time dynamics of strongly correlated fermions using auxiliary field quantum Monte Carlo

Spurred by recent technological advances, there is a growing demand for computational methods that can accurately predict the dynamics of correlated electrons. Such methods can provide much-needed theoretical insights into the electron dynamics probed via time-resolved spectroscopy experiments and observed in non-equilibrium ultracold atom experiments. In this article, we develop and benchmark a numerically exact Auxiliary Field Quantum Monte Carlo (AFQMC) method for modeling the dynamics of correlated electrons in real time. AFQMC has become a powerful method for predicting the ground state and finite temperature properties of strongly correlated systems mostly by employing constraints to control the sign problem. Our initial goal in this work is to determine how well AFQMC generalizes to real-time electron dynamics problems without constraints. By modeling the repulsive Hubbard model on different lattices and with differing initial electronic configurations, we show that real-time AFQMC is capable of accurately capturing long-lived electronic coherences beyond the reach of mean field techniques. While the times to which we can meaningfully model decrease with increasing correlation strength and system size as a result of the exponential growth of the dynamical phase problem, we show that our technique can model the short-time behavior of strongly correlated systems to very high accuracy. Crucially, we find that importance sampling, combined with a novel adaptive active space sampling technique, can substantially lengthen the times to which we can simulate. These results establish real-time AFQMC as a viable technique for modeling the dynamics of correlated electron systems and serve as a basis for future sampling advances that will further mitigate the dynamical phase problem.

Entropy in the Non-Fermi-Liquid Regime of the Doped 2D Hubbard Model

We study thermodynamic properties of the doped Hubbard model on the square lattice in the regime of strong charge and spin fluctuations at low temperatures near the metal-to-insulator crossover and obtain results with controlled accuracy using the diagrammatic Monte Carlo method directly in the thermodynamic limit. The behavior of the entropy reveals a non-Fermi-liquid state at sufficiently high interactions near half filling: A maximum in the entropy at nonzero doping develops as the coupling strength is increased, along with an inflection point, evidencing a metal to non-Fermi-liquid crossover. The specific heat exhibits additional distinctive features of a non-Fermi-liquid state. Measurements of the entropy can, therefore, be used as a probe of the state of the system in quantum simulation experiments with ultracold atoms in optical lattices.

Z_{2} Parton Phases in the Mixed-Dimensional t-J_{z} Model

We study the interplay of spin and charge degrees of freedom in a doped Ising antiferromagnet, where the motion of charges is restricted to one dimension. The phase diagram of this mixed-dimensional t-J_{z} model can be understood in terms of spinless chargons coupled to a Z_{2} lattice gauge field. The antiferromagnetic couplings give rise to interactions between Z_{2} electric field lines which, in turn, lead to a robust stripe phase at low temperatures. At higher temperatures, a confined meson-gas phase is found for low doping whereas at higher doping values, a robust deconfined chargon-gas phase is seen, which features hidden antiferromagnetic order. We confirm these phases in quantum Monte Carlo simulations. Our model can be implemented and its phases detected with existing technology in ultracold atom experiments. The critical temperature for stripe formation with a sufficiently high hole concentration is around the spin-exchange energy J_{z}, i.e., well within reach of current experiments.

Thermalization of a Trapped One-Dimensional Bose Gas via Diffusion

For a decade the fate of a one-dimensional gas of interacting bosons in an external trapping potential remained mysterious. We here show that whenever the underlying integrability of the gas is broken by the presence of the external potential, the inevitable diffusive rearrangements between the quasiparticles, quantified by the diffusion constants of the gas, eventually lead the system to thermalize at late times. We show that the full thermalizing dynamics can be described by the generalized hydrodynamics with diffusion and force terms, and we compare these predictions to numerical simulations. Finally, we provide an explanation for the slow thermalization rates observed in numerical and experimental settings: the hydrodynamics of integrable models is characterized by a continuity of modes, which can have arbitrarily small diffusion coefficients. As a consequence, the approach to thermalization can display prethermal plateau and relaxation dynamics with long polynomial finite-time corrections.

Imaging the Holon string of the Hubbard model

It has been a long-sought goal of quantum simulation to find answers to outstanding questions in condensed-matter physics. A famous example is finding the ground state and the excitations of the two-dimensional (2D) Hubbard model with strong repulsion below half-filling. This system is a doped antiferromagnet and is of great interest because of its possible relation to high-[Formula: see text] superconductors. Theoretically, the fermion excitations of this model are believed to split up into holons and spinons, and a moving holon is believed to leave behind it a string of "wrong" spins that mismatch with the antiferromagnetic background. Here, we show that the properties of the ground-state wavefunction and the holon excitation of the 2D Hubbard model can be revealed in unprecedented detail by using the imaging and the interference technique in atomic physics. They allow one to reveal the Marshall sign of the doped antiferromagnet. The region of wrong Marshall sign indicates the location of the holon string.

Doublon-Hole Correlations and Fluctuation Thermometry in a Fermi-Hubbard Gas

We report on the single atom and single site-resolved detection of the total density in a cold atom realization of the 2D Fermi-Hubbard model. Fluorescence imaging of doublons is achieved by splitting each lattice site into a double well, thereby separating atom pairs. Full density readout yields a direct measurement of the equation of state, including direct thermometry via the fluctuation-dissipation theorem. Site-resolved density correlations reveal the Pauli hole at low filling, and strong doublon-hole correlations near half filling. These are shown to account for the difference between local and nonlocal density fluctuations in the Mott insulator. Our technique enables the study of atom-resolved charge transport in the Fermi-Hubbard model, the site-resolved observation of molecules, and the creation of bilayer Fermi-Hubbard systems.

Superdiffusion from Emergent Classical Solitons in Quantum Spin Chains

Finite-temperature spin transport in the quantum Heisenberg spin chain is known to be superdiffusive, and has been conjectured to lie in the Kardar-Parisi-Zhang (KPZ) universality class. Using a kinetic theory of transport, we compute the KPZ coupling strength for the Heisenberg chain as a function of temperature, directly from microscopics; the results agree well with density-matrix renormalization group simulations. We establish a rigorous quantum-classical correspondence between the "giant quasiparticles" that govern superdiffusion and solitons in the classical continuous Landau-Lifshitz ferromagnet. We conclude that KPZ universality has the same origin in classical and quantum integrable isotropic magnets: a finite-temperature gas of low-energy classical solitons.

Robust Bilayer Charge Pumping for Spin- and Density-Resolved Quantum Gas Microscopy

Quantum gas microscopy has emerged as a powerful new way to probe quantum many-body systems at the microscopic level. However, layered or efficient spin-resolved readout methods have remained scarce as they impose strong demands on the specific atomic species and constrain the simulated lattice geometry and size. Here we present a novel high-fidelity bilayer readout, which can be used for full spin- and density-resolved quantum gas microscopy of two-dimensional systems with arbitrary geometry. Our technique makes use of an initial Stern-Gerlach splitting into adjacent layers of a highly stable vertical superlattice and subsequent charge pumping to separate the layers by 21  μm. This separation enables independent high-resolution images of each layer. We benchmark our method by spin- and density-resolving two-dimensional Fermi-Hubbard systems. Our technique furthermore enables the access to advanced entropy engineering schemes, spectroscopic methods, or the realization of tunable bilayer systems.

Spin and Charge Correlations across the Metal-to-Insulator Crossover in the Half-Filled 2D Hubbard Model

The 2D Hubbard model with nearest-neighbor hopping on the square lattice and an average of one electron per site is known to undergo an extended crossover from metallic to insulating behavior driven by proliferating antiferromagnetic correlations. We study signatures of this crossover in spin and charge correlation functions and present results obtained with controlled accuracy using the diagrammatic Monte Carlo approach in the range of parameters amenable to experimental verification with ultracold atoms in optical lattices. The qualitative changes in charge and spin correlations associated with the crossover are observed at well-separated temperature scales, which encase the intermediary regime of non-Fermi-liquid character, where local magnetic moments are formed and nonlocal fluctuations in both channels are essential.

Nagaoka ferromagnetism observed in a quantum dot plaquette

Engineered, highly controllable quantum systems are promising simulators of emergent physics beyond the simulation capabilities of classical computers¹. An important problem in many-body physics is itinerant magnetism, which originates purely from long-range interactions of free electrons and whose existence in real systems has been debated for decades^2,3. Here we use a quantum simulator consisting of a four-electron-site square plaquette of quantum dots⁴ to demonstrate Nagaoka ferromagnetism⁵. This form of itinerant magnetism has been rigorously studied theoretically^6-9 but has remained unattainable in experiments. We load the plaquette with three electrons and demonstrate the predicted emergence of spontaneous ferromagnetic correlations through pairwise measurements of spin. We find that the ferromagnetic ground state is remarkably robust to engineered disorder in the on-site potentials and we can induce a transition to the low-spin state by changing the plaquette topology to an open chain. This demonstration of Nagaoka ferromagnetism highlights that quantum simulators can be used to study physical phenomena that have not yet been observed in any experimental system. The work also constitutes an important step towards large-scale quantum dot simulators of correlated electron systems.

Extended Crossover from a Fermi Liquid to a Quasiantiferromagnet in the Half-Filled 2D Hubbard Model

The ground state of the Hubbard model with nearest-neighbor hopping on the square lattice at half filling is known to be that of an antiferromagnetic (AFM) band insulator for any on-site repulsion. At finite temperature, the absence of long-range order makes the question of how the interaction-driven insulator is realized nontrivial. We address this problem with controlled accuracy in the thermodynamic limit using self-energy diagrammatic determinant Monte Carlo and dynamical cluster approximation methods and show that development of long-range AFM correlations drives an extended crossover from Fermi liquid to insulating behavior in the parameter regime that precludes a metal-to-insulator transition. The intermediate crossover state is best described as a non-Fermi liquid with a partially gapped Fermi surface.

Strange metallicity in the doped Hubbard model

Strange or bad metallic transport, defined by incompatibility with the conventional quasiparticle picture, is a theme common to many strongly correlated materials, including high-temperature superconductors. The Hubbard model represents a minimal starting point for modeling strongly correlated systems. Here we demonstrate strange metallic transport in the doped two-dimensional Hubbard model using determinantal quantum Monte Carlo calculations. Over a wide range of doping, we observe resistivities exceeding the Mott-Ioffe-Regel limit with linear temperature dependence. The temperatures of our calculations extend to as low as 1/40 of the noninteracting bandwidth, placing our findings in the degenerate regime relevant to experimental observations of strange metallicity. Our results provide a foundation for connecting theories of strange metals to models of strongly correlated materials.

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.

Conductivity Spectrum of Ultracold Atoms in an Optical Lattice

We measure the conductivity of neutral fermions in a cubic optical lattice. Using in situ fluorescence microscopy, we observe the alternating current resultant from a single-frequency uniform force applied by displacement of a weak harmonic trapping potential. In the linear response regime, a neutral-particle analog of Ohm's law gives the conductivity as the ratio of total current to force. For various lattice depths, temperatures, interaction strengths, and fillings, we measure both real and imaginary conductivity, up to a frequency sufficient to capture the transport dynamics within the lowest band. The spectral width of the real conductivity reveals the current dissipation rate in the lattice, and the integrated spectral weight is related to thermodynamic properties of the system through a sum rule. The global conductivity decreases with increased band-averaged effective mass, which at high temperatures approaches a T-linear regime. Relaxation of current is observed to require a finite lattice depth, which breaks Galilean invariance and enables damping through collisions between fermions.

Transport with strong interactions

Motion of spin and charge is explored with cold-atom quantum simulations

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.