Quantum simulation experiments have started to explore regimes that are not accessible with exact numerical methods. To probe these systems and enable new physical insights, the need for measurement protocols arises that can bridge the gap to solid-state experiments, and at the same time make optimal use of the capabilities of quantum simulation experiments. Here we propose applying time-dependent photoemission spectroscopy, an established tool in solid-state systems, in cold atom quantum simulators. Concretely, we suggest combining the method with large magnetic field gradients, unattainable in experiments on real materials, to drive Bloch oscillations of spinons, the emergent quasiparticles of spin liquids. We show in exact diagonalization simulations of the one-dimensional t-J model with a single hole that the spinons start to populate previously unoccupied states in an effective band structure, thus allowing us to visualize states invisible in the equilibrium spectrum. The dependence of the spectral function on the time after the pump pulse reveals collective interactions among spinons. In numerical simulations of small two-dimensional systems, spectral weight appears at the ground-state energy at momentum q = (pi, pi), where the equilibrium spectral response is strongly suppressed up to higher energies, indicating a possible route toward solving the mystery of the Fermi arcs in the cuprate materials.

In the last decade, quantum simulators, and in particular cold atoms in optical lattices, have emerged as a valuable tool to study strongly correlated quantum matter. These experiments are now reaching regimes that are numerically difficult or impossible to access. In particular they have started to fulfill a promise which has contributed significantly to defining and shaping the field of cold atom quantum simulations, namely the exploration of doped and frustrated quantum magnets and the search for the origins of high-temperature superconductivity in the fermionic Hubbard model. Despite many future challenges lying ahead, such as the need to further lower the experimentally accessible temperatures, remarkable studies have already been conducted. Among them, spin-charge separation in one-dimensional systems has been demonstrated, extended-range antiferromagnetism in two-dimensional systems has been observed, connections to modern day large-scale numerical simulations were made, and unprecedented comparisons with microscopic trial wavefunctions have been carried out at finite doping. In many regards, the field has acquired new realms, putting old ideas to a new test and producing new insights and inspiration for the next generation of physicists. In the first part of this paper, we review the results achieved in cold atom realizations of the Fermi-Hubbard model in recent years. We put special emphasis on the new probes available in quantum gas microscopes, such as higher-order correlation functions, full counting statistics, the ability to study far-from -equilibrium dynamics, machine learning and pattern recognition of instantaneous snapshots of the many-body wavefunction, and access to non-local correlators. Our review is written from a theoretical perspective, but aims to provide basic understanding of the experimental procedures. We cover one- dimensional systems, where the phenomenon of spin-charge separation is ubiquitous, and two-dimensional systems where we distinguish between situations with and without doping. Throughout, we focus on the strong coupling regime where the Hubbard inter-actions U dominate and connections to t - J models can be justified. In the second part of this paper, with the stage set and the current state of the field in mind, we propose a new direction for cold atoms to explore: namely mixed-dimensional bilayer systems, where the charge motion is restricted to individual layers which remain coupled through spin-exchange. These systems can be directly realized experimentally and we argue that they have a rich phase diagram, potentially including a strongly correlated BEC-to-BCS cross-over and regimes with different superconducting order parameters, as well as complex parton phases and possibly even analogs of tetraquark states. In particular, we propose a novel, strong pairing mechanism in these systems, which puts the formation of hole pairs at experimentally accessible, elevated temperatures within reach. Ultimately we propose to explore how the physics of the mixed-dimensional bilayer system can be connected to the rich phenomenology of the single-layer Hubbard model. (C) 2021 Elsevier Inc. All rights reserved.

Understanding the nature of charge carriers in doped Mott insulators holds the key to unravelling puzzling properties of strongly correlated electron systems, including cuprate superconductors. Several theoretical models suggested that dopants can be understood as bound states of partons, the analogues of quarks in high-energy physics. However, direct signatures of spinon-chargon bound states are lacking, both in experiment and theory. Here we propose a rotational variant of angle-resolved photo-emission spectroscopy (ARPES) and calculate rotational spectra numerically using the density-matrix renormalization group. We identify long-lived rotational resonances for an individual dopant, which we interpret as a direct indicator of the microscopic structure of spinon-chargon bound states. Similar to Regge trajectories reflecting the quark structure of mesons, we establish a linear dependence of the rotational energy on the superexchange coupling. The rotational peaks we find are strongly suppressed in standard ARPES spectra, but we suggest a multiphoton extension of ARPES which allows us to access rotational spectra. Our findings suggest that multiphoton spectroscopy experiments should provide new insights into emergent universal features of strongly correlated electron systems.

Confinement is an ubiquitous phenomenon when matter couples to gauge fields, which manifests itself in a linear string potential between two static charges. Although gauge fields can be integrated out in one dimension, they can mediate nonlocal interactions which in turn influence the paradigmatic Luttinger liquid properties. However, when the charges become dynamical and their densities finite, understanding confinement becomes challenging. Here we show that confinement in 1D Z(2) lattice gauge theories, with dynamical matter fields and arbitrary densities, is related to translational symmetry breaking in a nonlocal basis. The exact transformation to this string-length basis leads us to an exact mapping of Luttinger parameters reminiscent of a Luther-Emery rescaling. We include the effects of local, but beyond contact, interactions between the matter particles, and show that confined mesons can form a Mott-insulating state when the deconfined charges cannot. While the transition to the Mott state cannot be detected in the Green's function of the charges, we show that the metallic state is characterized by hidden off-diagonal quasi-long-range order. Our predictions provide new insights to the physics of confinement of dynamical charges, and can be experimentally addressed in Rydberg-dressed quantum gases in optical lattices.

The competition between antiferromagnetism and hole motion in two-dimensional Mott insulators lies at the heart of a doping-dependent transition from an anomalous metal to a conventional Fermi liquid. We observe such a crossover in Fermi-Hubbard systems on a cold-atom quantum simulator and reveal the transformation of multipoint correlations between spins and holes upon increasing doping at temperatures around the superexchange energy. Conventional observables, such as spin susceptibility, are furthermore computed from the microscopic snapshots of the system. Starting from a magnetic polaron regime, we find the system evolves into a Fermi liquid featuring incommensurate magnetic fluctuations and fundamentally altered correlations. The crossover is completed for hole dopings around 30%. Our work benchmarks theoretical approaches and discusses possible connections to lowertemperature phenomena.

Analysis of higher-order correlation functions has become a powerful tool for investigating interacting many-body systems in quantum simulators, such as quantum gas microscopes. Experimental measurements of mixed spin-charge correlation functions in the 2D Hubbard have been used to study equilibrium properties of magnetic polarons and to identify coherent and incoherent regimes of their dynamics. In this paper we consider theoretically an extension of this technique to systems which use a pinning potential to reduce the mobility of a single dopant in the Mott insulating regime of the 2D Hubbard model. We find that localization of the dopant has a dramatic effect on its magnetic dressing. The connected third order spin correlations are weakened in the case of a mobile hole but strengthened near an immobile hole. In the case of the fifth-order correlation function, we find that its bare value has opposite signs in cases of the mobile and of fully pinned dopant, whereas the connected part is similar for both cases. We study suppression of higher-order correlators by thermal fluctuations and demonstrate that they can be observed up to temperatures comparable to the spin-exchange energy J. We discuss implications of our results for understanding the interplay of spin and charge in doped Mott insulators.

Understanding the interplay between charge and spin and its effects on transport is a ubiquitous challenge in quantum many-body systems. In the Fermi-Hubbard model, this interplay is thought to give rise to magnetic polarons, whose dynamics may explain emergent properties of quantum materials such as high-temperature superconductivity. In this work, we use a cold-atom quantum simulator to directly observe the formation dynamics and subsequent spreading of individual magnetic polarons. Measuring the density- and spin-resolved evolution of a single hole in a 2D Hubbard insulator with short-range antiferromagnetic correlations reveals fast initial delocalization and a dressing of the spin background, indicating polaron formation. At long times, we find that dynamics are slowed down by the spin exchange time, and they are compatible with a polaronic model with strong density and spin coupling. Our work enables the study of out-of-equilibrium emergent phenomena in the Fermi-Hubbard model, one dopant at a time.

We investigate the crossover of the impurity-induced dynamics, in trapped one-dimensional Bose polarons subject to radio frequency (RF) pulses of varying intensity, from an adiabatic to a diabatic regime. Utilizing adiabatic pulses for either weak repulsive or attractive impurity-medium interactions, a multitude of polaronic excitations or mode-couplings of the impurity-bath interaction with the collective breathing motion of the bosonic medium are spectrally resolved. We find that for strongly repulsive impurity-bath interactions, a temporal orthogonality catastrophe manifests in resonances in the excitation spectra where impurity coherence vanishes. When two impurities are introduced, impurity–impurity correlations, for either attractive or strong repulsive couplings, induce a spectral shift of the resonances with respect to the single impurity. For a heavy impurity, the polaronic peak is accompanied by a series of equidistant side-band resonances, related to interference of the impurity spin dynamics and the sound waves of the bath. In all cases, we enter the diabatic transfer regime for an increasing bare Rabi frequency of the RF field with a Lorentzian spectral shape featuring a single polaronic resonance. The findings in this work on the effects of external trap, RF pulse and impurity–impurity interaction should have implications for the new generations of cold-atom experiments.

Topological states of matter, such as fractional quantum Hall states, are an active field of research due to their exotic excitations. In particular, ultracold atoms in optical lattices provide a highly controllable and adaptable platform to study such new types of quantum matter. However, finding a clear route to realize non-Abelian quantum Hall states in these systems remains challenging. Here we use the density-matrix renormalization-group (DMRG) method to study the Hofstadter-Bose-Hubbard model at filling factor ν=1 and find strong indications that at α=1/6 magnetic flux quanta per plaquette the ground state is a lattice analog of the continuum non-Abelian Pfaffian. We study the on-site correlations of the ground state, which indicate its paired nature at ν=1, and find an incompressible state characterized by a charge gap in the bulk. We argue that the emergence of a charge density wave on thin cylinders and the behavior of the two- and three-particle correlation functions at short distances provide evidence for the state being closely related to the continuum Pfaffian. The signatures discussed in this letter are accessible in current cold atom experiments and we show that the Pfaffian-like state is readily realizable in few-body systems using adiabatic preparation schemes.

Traditionally, one- and two-point correlation functions are used to characterize many-body systems. In strongly correlated quantum materials, such as the doped 2D Fermi-Hubbard system, these may no longer be sufficient, because higher-order correlations are crucial to understanding the character of the many-body system and can be numerically dominant. Experimentally, such higher-order correlations have recently become accessible in ultracold atom systems. Here, we reveal strong non-Gaussian correlations in doped quantum antiferromagnets and show that higher-order correlations dominate over lower-order terms. We study a single mobile hole in the t - J model using the density matrix renormalization group and reveal genuine fifth-order correlations which are directly related to the mobility of the dopant. We contrast our results to predictions using models based on doped quantum spin liquids which feature significantly reduced higher-order correlations. Our predictions can be tested at the lowest currently accessible temperatures in quantum simulators of the 2D Fermi-Hubbard model. Finally, we propose to experimentally study the same fifth-order spin-charge correlations as a function of doping. This will help to reveal the microscopic nature of charge carriers in the most debated regime of the Hubbard model, relevant for understanding high-T-c superconductivity.

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.

Kitaev's toric code is an exactly solvable model with Z2-topological order, which has potential applications in quantum computation and error correction. However, a direct experimental realization remains an open challenge. Here, we propose a building block for Z2 lattice gauge theories coupled to dynamical matter and demonstrate how it allows for an implementation of the toric-code ground state and its topological excitations. This is achieved by introducing separate matter excitations on individual plaquettes, whose motion induce the required plaquette terms. The proposed building block is realized in the second-order coupling regime and is well suited for implementations with superconducting qubits. Furthermore, we propose a pathway to prepare topologically non-trivial initial states during which a large gap on the order of the underlying coupling strength is present. This is verified by both analytical arguments and numerical studies. Moreover, we outline experimental signatures of the ground-state wavefunction and introduce a minimal braiding protocol. Detecting a π-phase shift between Ramsey fringes in this protocol reveals the anyonic excitations of the toric-code Hamiltonian in a system with only three triangular plaquettes. Our work paves the way for realizing non-Abelian anyons in analog quantum simulators.

Tremendous recent progress in the quantum simulation of the Hubbard model paves the way to controllably study doped antiferromagnetic Mott insulators. Motivated by these experimental advancements, we numerically study the real-time dynamics of a single hole created in an antiferromagnet on a square lattice, as described by the t-J model. Initially, the hole spreads ballistically with a velocity proportional to the hopping matrix element. At intermediate to long times, the dimensionality as well as the spin background determine the hole dynamics. A hole created in the ground state of a two dimensional (2D) quantum antiferromagnet propagates again ballistically at long times but with a velocity proportional to the spin exchange coupling, showing the formation of a magnetic polaron. We provide an intuitive explanation of this dynamics in terms of a parton construction, which leads to a good quantitative agreement with the numerical tensor network state simulations. In the limit of infinite temperature and no spin exchange couplings, the dynamics can be approximated by a quantum random walk on a Bethe lattice with coordination number

z

x303;

4

Adding Ising interactions corresponds to an effective disordered potential, which can dramatically slow down the hole propagation, consistent with subdiffusive dynamics. The study of the hole dynamics paves the way for understanding the microscopic constituents of this strongly correlated quantum state.

We show that ultracold fermions in an artificial magnetic field open up a new window to the physics of the spinful fractional quantum Hall (FQH) effect. We numerically study the lowest energy states of strongly interacting few-fermion systems in rapidly rotating optical microtraps. We find that skyrmion-like ground states with locally ferromagnetic, long-range spin textures emerge. To realize such states experimentally, rotating microtraps with higher-order angular momentum components may be used to prepare fermionic particles in a lowest Landau level. We find parameter regimes in which skyrmion-like ground states should be accessible in current experiments and demonstrate an adiabatic pathway for their preparation in a rapidly rotating harmonic trap. The addition of long range interactions will lead to an even richer interplay between spin textures and FQH physics.

We study higher order topology in the presence of strong interactions in a two-dimensional, experimentally accessible superlattice Bose-Hubbard model with alternating hoppings and strong on-site repulsion. We evaluate the phase diagram of the model around half-filling using the density renormalization group ansatz and find two gapped phases separated by a gapless superfluid region. We demonstrate that the gapped states realize two distinct higher order symmetry protected topological phases, which are protected by a combination of charge conservation and C4 lattice symmetry. The phases are distinguished in terms of a many-body topological invariant and a quantized, experimentally accessible fractional corner charge that is robust against arbitrary, symmetry preserving edge manipulations. We support our claims by numerically studying the full counting statistics of the corner charge, finding a sharp distribution peaked around the quantized values. Our results allow for a direct comparison with experiments and represent a confirmation of theoretically predicted higher order topology in a strongly interacting system. Experimentally, the fractional corner charge can be observed in ultracold atomic settings using state of the art quantum gas microscopy.

Angle-resolved photoemission spectroscopy (ARPES) has revealed peculiar properties of mobile dopants in correlated antiferromagnets (AFMs). But, describing them theoretically, even in simplified toy models, remains a challenge. Here, we study ARPES spectra of a single mobile hole in the t-J model. Recent progress in the microscopic description of mobile dopants allows us to use a geometric decoupling of spin and charge fluctuations at strong couplings, from which we conjecture a one-to-one relation of the one-dopant spectral function and the spectrum of a constituting spinon in the undoped parent AFM. We thoroughly test this hypothesis for a single hole doped into a two-dimensional Heisenberg AFM by comparing our semianalytical predictions to previous quantum Monte Carlo results and our large-scale time-dependent matrix product state calculations of the spectral function. Our conclusion is supported by a microscopic trial wave function describing spinon-chargon bound states, which captures the momentum and t/J dependence of the quasiparticle residue. From our conjecture we speculate that ARPES measurements in the pseudogap phase of cuprates may directly reveal the Dirac-fermion nature of the constituting spinons. Specifically, we demonstrate that our trial wave function provides a microscopic explanation for the sudden drop of spectral weight around the nodal point associated with the formation of Fermi arcs, assuming that additional frustration suppresses long-range AFM ordering. We benchmark our results by studying the crossover from two to one dimension, where spinons and chargons are confined and deconfined, respectively.

We introduce a Ramsey pulse scheme which extracts the non-Hermitian Hamiltonian associated with an arbitrary Lindblad dynamics. We propose a related protocol to measure via interferometry a generalized Loschmidt echo of a generic state evolving in time with the non-Hermitian Hamiltonian itself, and we apply the scheme to a one-dimensional weakly interacting Bose gas coupled to a stochastic atomic impurity. The Loschmidt echo is mapped into a functional integral from which we calculate the long-time decohering dynamics at arbitrary impurity strengths. For strong dissipation we uncover the phenomenology of a quantum many-body Zeno effect: Corrections to the decoherence exponent resulting from the impurity self-energy become purely imaginary, in contrast to the regime of small dissipation where they instead enhance the decay of quantum coherences. Our results illustrate the prospects for experiments employing Ramsey interferometry to study dissipative quantum impurities in condensed matter and cold-atom systems.

We investigate a quantum many-body lattice system of one-dimensional spinless fermions interacting with a dynamical Z(2) gauge field. The gauge field mediates long-range attraction between fermions resulting in their confinement into bosonic dimers. At strong coupling we develop an exactly solvable effective theory of such dimers with emergent constraints. Even at generic coupling and fermion density, the model can be rewritten as a local spin chain. Using the density matrix renormalization group the system is shown to form a Luttinger liquid, indicating the emergence of fractionalized excitations despite the confinement of lattice fermions. In a finite chain we observe the doubling of the period of Friedel oscillations which paves the way towards experimental detection of confinement in this system. We discuss the possibility of a Mott phase at the commensurate filling 2/3.

We introduce three numerical methods for characterizing the topological phases of three-dimensional multiband Hubbard models based on twisted boundary conditions, Wilson loops, as well as the local topological marker. We focus on the half-filled, three-dimensional time-reversal-invariant Hofstadter model with finite spin-orbit coupling. Besides the weak and strong topological insulator phases we find a nodal line semimetal in the parameter regime between the two three-dimensional topological insulator phases. Using dynamical mean-field theory combined with the topological Hamiltonian approach we find stabilization of these three-dimensional topological states due to the Hubbard interaction. We study surface states which exhibit an asymmetry between left and right surfaces originating from the broken parity symmetry of the system. Our results set the stage for further research on inhomogeneous three-dimensional topological systems, proximity effects, topological Mott insulators, nontrivially linked nodal line semimetals, and circuit-based quantum simulators.

Dominant multiparticle interactions can give rise to exotic physical phases with anyonic excitations and phase transitions without local order parameters. In spin systems with a global SU(N) symmetry, cyclic ring-exchange couplings constitute the first higher-order interaction in this class. In this Letter, we propose a protocol showing how SU(N)-invariant multibody interactions can be implemented in optical tweezer arrays. We utilize the flexibility to rearrange the tweezer configuration on short timescales compared to the typical lifetimes, in combination with strong nonlocal Rydberg interactions. As a specific example, we demonstrate how a chiral cyclic ring-exchange Hamiltonian can be implemented in a two-leg ladder geometry. We study its phase diagram using density-matrix renormalization group simulations and identify phases with dominant vector chirality, a ferromagnet, and an emergent spin-1 Haldane phase. We also discuss how the proposed protocol can be utilized to implement the strongly frustrated J–Q model, a candidate for hosting a deconfined quantum critical point.

Infinite projected entangled pair states (iPEPS) provide a convenient variational description of infinite, translationally-invariant two-dimensional quantum states. However, the simulation of local excitations is not directly possible due to the translationally-invariant ansatz. Furthermore, as iPEPS are either identical or orthogonal, expectation values between different states as required during the evaluation of non-equal-time correlators are ill-defined. Here, we show that by introducing auxiliary states on each site, it becomes possible to simulate both local excitations and evaluate non-equal-time correlators in an iPEPS setting under real-time evolution. We showcase the method by simulating the t - J model after a single hole has been placed in the half-filled antiferromagnetic background and evaluating both return probabilities and spin correlation functions, as accessible in quantum gas microscopes.

Elementary particles carry several quantum numbers, such as charge and spin. However, in an ensemble of strongly interacting particles, the emerging degrees of freedom can fundamentally differ from those of the individual constituents. For example, one-dimensional systems are described by independent quasiparticles carrying either spin (spinon) or charge (holon). Here, we report on the dynamical deconfinement of spin and charge excitations in real space after the removal of a particle in Fermi-Hubbard chains of ultracold atoms. Using space- and time-resolved quantum gas microscopy, we tracked the evolution of the excitations through their signatures in spin and charge correlations. By evaluating multipoint correlators, we quantified the spatial separation of the excitations in the context of fractionalization into single spinons and holons at finite temperatures.

From the standard model of particle physics to strongly correlated electrons, various physical settings are formulated in terms of matter coupled to gauge fields. Quantum simulations based on ultracold atoms in optical lattices provide a promising avenue to study these complex systems and unravel the underlying many-body physics. Here, we demonstrate how quantized dynamical gauge fields can be created in mixtures of ultracold atoms in optical lattices, using a combination of coherent lattice modulation with strong interactions. Specifically, we propose implementation of Z2 lattice gauge theories coupled to matter, reminiscent of theories previously introduced in high-temperature superconductivity. We discuss a range of settings from zero-dimensional toy models to ladders featuring transitions in the gauge sector to extended two-dimensional systems. Mastering lattice gauge theories in optical lattices constitutes a new route toward the realization of strongly correlated systems, with properties dictated by an interplay of dynamical matter and gauge fields.

We unravel the nonequilibrium correlated quantum quench dynamics of an impurity traveling through a harmonically confined Bose–Einstein condensate in one-dimension. For weak repulsive interspecies interactions the impurity oscillates within the bosonic gas. At strong repulsions and depending on its prequench position the impurity moves towards an edge of the bosonic medium and subsequently equilibrates. This equilibration being present independently of the initial velocity, the position and the mass of the impurity is inherently related to the generation of entanglement in the many-body system. Focusing on attractive interactions the impurity performs a damped oscillatory motion within the bosonic bath, a behavior that becomes more evident for stronger attractions. To elucidate our understanding of the dynamics an effective potential picture is constructed. The effective mass of the emergent quasiparticle is measured and found to be generically larger than the bare one, especially for strong attractions. In all cases, a transfer of energy from the impurity to the bosonic medium takes place. Finally, by averaging over a sample of simulated in situ single-shot images we expose how the single-particle density distributions and the two-body interspecies correlations can be probed.

Quantum simulation has the potential to investigate gauge theories in strongly-interacting regimes, which are up to now inaccessible through conventional numerical techniques. Here, we take a first step in this direction by implementing a Floquet-based method for studying Z2 lattice gauge theories using two-component ultracold atoms in a double-well potential. For resonant periodic driving at the on-site interaction strength and an appropriate choice of the modulation parameters, the effective Floquet Hamiltonian exhibits Z2 symmetry. We study the dynamics of the system for different initial states and critically contrast the observed evolution with a theoretical analysis of the full time-dependent Hamiltonian of the periodically-driven lattice model. We reveal challenges that arise due to symmetry-breaking terms and outline potential pathways to overcome these limitations. Our results provide important insights for future studies of lattice gauge theories based on Floquet techniques.

In the presence of symmetries, one-dimensional quantum systems can exhibit topological order, which in many cases can be characterized by a quantized value of the many-body geometric Zak or Berry phase. We establish that this topological Zak phase is directly related to the Zak phase of an elementary quasiparticle excitation in the system. By considering various systems, we establish this connection for a number of different interacting phases including the Su-Schrieffer-Heeger model, p-wave topological superconductors, and the Haldane chain. Crucially, in contrast to the bulk many-body Zak phase associated with the ground state of such systems, the topological invariant associated with quasiparticle excitations (above this ground state) exhibits a more natural route for direct experimental detection. To this end, we build upon recent work [F. Grusdt, et al., Nat. Commun. 7, 11994 (2016)] and demonstrate that mobile quantum impurities can be used, in combination with Ramsey interferometry and Bloch oscillations, to directly measure these quasiparticle topological invariants. Finally, a concrete experimental realization of our protocol for dimerized Mott insulators in ultracold atomic systems is discussed and analyzed.

Understanding strongly correlated quantum many-body states is one of the most difficult challenges in modern physics. For example, there remain fundamental open questions on the phase diagram of the Hubbard model, which describes strongly correlated electrons in solids. In this work, we realize the Hubbard Hamiltonian and search for specific patterns within the individual images of many realizations of strongly correlated ultracold fermions in an optical lattice. Upon doping a cold-atom antiferromagnet, we find consistency with geometric strings, entities that may explain the relationship between hole motion and spin order, in both pattern-based and conventional observables. Our results demonstrate the potential for pattern recognition to provide key insights into cold-atom quantum many-body systems.

Quantum gas microscopes for ultracold atoms can provide high-resolution real-space snapshots of complex many-body systems. We implement machine learning to analyse and classify such snapshots of ultracold atoms. Specifically, we compare the data from an experimental realization of the two-dimensional Fermi–Hubbard model to two theoretical approaches: a doped quantum spin liquid state of resonating valence bond type (1,2), and the geometric string theory (3,4), describing a state with hidden spin order. This technique considers all available information without a potential bias towards one particular theory by the choice of an observable and can therefore select the theory that is more predictive in general. Up to intermediate doping values, our algorithm tends to classify experimental snapshots as geometric-string-like, as compared to the doped spin liquid. Our results demonstrate the potential for machine learning in processing the wealth of data obtained through quantum gas microscopy for new physical insights.