Kitaev's toric code has become one of the most studied models in physics and is highly relevant to the fields of both quantum error correction and condensed matter physics. Most notably, it is the simplest known model hosting an extended, deconfined topological bulk phase. To this day, it remains challenging to reliably and robustly probe topological phases, as many state-of-the-art order parameters are sensitive to specific models and even specific parameter regimes. With the emergence of powerful quantum simulators which are approaching the regimes of topological bulk phases, there is a timely need for experimentally accessible order parameters. Here we study the ground state physics of the parallel field toric code on the honeycomb, triangular, and cubic lattices using continuous-time quantum Monte Carlo. By extending the concept of experimentally accessible percolation-inspired order parameters (POPs) we show that electric and magnetic anyons are independently confined on the honeycomb and triangular lattices, unlike on the square lattice. Our work manifestly demonstrates that, even in the ground state, we must make a distinction between topological order and (de-)confinement. Moreover, we report multicritical points in the aforementioned confinement phase diagrams. Finally, we map out the topological phase diagrams on the honeycomb, triangular, and cubic lattices and compare the performance of the POPs with other topological order parameters. Our work paves the way for studies of confinement involving dynamical matter and the associated multicritical points in contemporary quantum simulation platforms for Z2 lattice gauge theories.

Topologically ordered quantum systems give rise to anyonic quasiparticles, whose controlled braiding operations form the foundation of topological quantum computation. Traditionally, studies of anyons have relied on edge-state interferometry, leaving the direct detection and manipulation of anyons in the bulk a major experimental challenge. Here, we propose and theoretically investigate a pathway toward this goal by demonstrating that a long-lived, optically generated interlayer exciton can bind to a quasihole in a fractional quantum Hall state, forming a novel composite excitation: the anyon-trion. Using exact diagonalization techniques, we reveal that anyon-trions exhibit millielectronvolt-scale binding energies and a linear dependence on the fractional charge of the quasihole. This scaling offers a powerful means to optically extract the quasihole's fractional charge through measurable shifts in exciton resonances. We outline a feasible experimental implementation via photoluminescence spectroscopy in a quantum twisting microscope setup, providing a promising route for the direct optical observation of anyon-trions within the bulk.

"In strongly correlated materials, interacting electrons are entangled and form collective quantum states, resulting in rich low-temperature phase diagrams. Notable examples include cuprate superconductors, in which superconductivity emerges at low doping out of an unusual ""pseudogap"" metallic state above the critical temperature. The Fermi-Hubbard model, describing a wide range of phenomena associated with strong electron correlations, still offers major computational challenges despite its simple formulation. In this context, ultracold atoms quantum simulators have provided invaluable insights into the microscopic nature of correlated quantum states. Here, we use a quantum gas microscope Fermi-Hubbard simulator to explore a wide range of dopings and temperatures in a regime where a pseudogap is known to develop. By measuring multipoint correlation functions up to fifth order, we uncover a universal scaling behavior in magnetic and higher-order spin-charge correlations characterized by a doping-dependent temperature scale. Accurate comparisons with determinant Quantum Monte Carlo and Minimally Entangled Typical Thermal States simulations confirm that this temperature scale is comparable to the pseudogap temperature T & lowast,.. Our quantitative findings reveal a qualitative behavior of magnetic properties and spin-charge correlations in an emergent pseudogap and pave the way toward the exploration of charge pairing and collective phenomena expected at lower temperatures."

Lattice gauge theories (LGTs) provide valuable insights into problems in strongly correlated many-body systems. Confinement which arises when matter is coupled to gauge fields is just one of the open problems, where LGT formalism can explain the underlying mechanism. However, coupling gauge fields to dynamical charges complicates the theoretical and experimental treatment of the problem. Developing a simplified meanfield theory is thus one of the ways to gain new insights into these complicated systems. Here we develop a mean-field theory of a paradigmatic 1+1D Z2 lattice gauge theory with superconducting pairing term, the gauged Kitaev chain, by decoupling charge and Z2 fields while enforcing the Gauss law on the mean-field level. We first determine the phase diagram of the original model in the context of confinement, which allows us to identify the symmetry-protected topological transition in the Kitaev chain as a confinement transition. We then compute the phase diagram of the effective mean-field theory, which correctly captures the main features of the original LGT. This is furthermore confirmed by the Green's function results and a direct comparison of the ground state energy. This simple LGT can be implemented in state-of-the art cold atom experiments. We thus also consider string-length histograms and the electric field polarization, which are easily accessible quantities in experimental setups and show that they reliably capture the various phases.

We address the Bose polaron problem of a mobile impurity interacting strongly with a host Bose-Einstein condensate through a Feshbach resonance. On the repulsive side at strong couplings, theoretical approaches predict two distinct polaron branches corresponding to attractive and repulsive polarons, but it remains unclear how the two are related. This is partly due to the challenges resulting from a competition of strongly attractive (destabilizing) impurity-boson interactions with weakly repulsive (stabilizing) boson-boson interactions, whose interplay is difficult to describe with contemporary theoretical methods. Here, we develop a powerful variational framework that combines Gaussian correlations among impurity-boson scattering states, including up to an infinite number of bosonic excitations, with exact non-Gaussian correlations among bosons occupying an impurity-boson bound state. This variational scheme enables a full treatment of strong nonlinearities arising in the Feshbach molecule on the repulsive side of the resonance. Within this framework, we demonstrate that the interplay of impurity-induced instability and stabilization by repulsive boson-boson interactions results in a discrete set of metastable many-body bound states at intermediate energies between the attractive and repulsive polaron branches. These states exhibit strong quantum statistical characteristics in the form of non-Gaussian quantum correlations, requiring nonperturbative beyond mean-field treatments for their characterization. Furthermore, these many-body bound states have sizable molecular spectral weights, accessible via molecular spectroscopy techniques. This work provides a unified theory of attractive and repulsive Bose polarons on the repulsive side of the Feshbach resonance.

Due to significant progress in quantum gas microscopy in recent years, there is a rapidly growing interest in real-space properties of single mobile dopands created in correlated antiferromagnetic (AFM) Mott insulators. However, a detailed numerical description remains challenging, even for simple toy models. As a consequence, previous numerical simulations for large systems were largely limited to T = 0. To provide guidance for cold-atom experiments, numerical calculations at finite temperature are required. Here, we numerically study the real-time properties of a single mobile hole in the 2D t-J model at finite temperature and draw a comparison to features observed at T = 0. We find that a three-stage process of hole motion, which was reported at T = 0, is valid even at finite temperature. However, already at low temperatures, the average hole velocity at long times is not simply proportional to the spin coupling, contrary to the T = 0 behavior. Comparing our finite-temperature numerical results with the experimental data from quantum gas microscopy we find a qualitative disagreement: in experiment, hole spreading speeds up with increasing J/t, while in our numerics it slows down. The latter is consistent with the numerical findings previously reported at T = 0.

The physics of strongly correlated fermions described by Hubbard or t-J models in the underdoped regime- relevant for high-temperature superconductivity in cuprate compounds-remains a subject of ongoing debate. In particular, the nature of charge carriers in this regime is poorly understood, in part due to the unusual properties of their spectral function. In this Letter, we present unbiased numerical results for the one-hole spectral function in a t-J model at finite temperatures. Our study provides valuable insights into the underlying physics of magnetic (or spin-) polaron formation in a doped antiferromagnet (AFM). For example, we find how the suppression of spectral weight outside the magnetic Brillouin zone-a precursor of Fermi arc formation-disappears with increasing temperature, revealing nearly-deconfined spinon excitations of the undoped AFM. The pristine setting we consider can be directly explored using quantum simulators. Our calculations demonstrate that coherent quasiparticle peaks associated with magnetic polarons can be observed up to temperatures T > J above the spin-exchange J, routinely obtained in such experiments. This paves the way for future studies of the fate of magnetic polarons in the pseudogap phase.

Since the discovery of unconventional superconductivity in cuprates, unraveling the pairing mechanism of charge carriers in doped antiferromagnets has been a long-standing challenge. Motivated by the discovery of high-Tc superconductivity in nickelate bilayer La3Ni2O7 (LNO), we study a minimal mixed dimensional t-J model supplemented with a repulsive Coulomb interaction V. When hole doped, previous numerical simulations revealed that the system exhibits strong binding energies, with a phenomenology resembling a BCS-to-BEC crossover accompanied by a Feshbach resonance between two distinct types of charge carriers. Here, we perform a mean-field analysis that enables a direct observation of the BCS-to-BEC crossover as well as microscopic insights into the crossover region and the pairing symmetry for two-dimensional bilayers. We benchmark our mean-field description by comparing it to density-matrix renormalization group simulations in quasi-one-dimensional settings and find remarkably good agreement. For the two-dimensional system relevant to LNO, our mean-field calculations predict a BCS pairing gap with an extended s-wave symmetry, directly resulting from the pairing mechanism's Feshbach origin. Our analysis hence gives insights into pairing in unconventional superconductors and, further, can be tested in currently available ultracold atom experiments.

We theoretically investigate the stationary properties of a spin-1/2 impurity immersed in a one-dimensional confined Bose gas. In particular, we consider coherently coupled spin states with an external field, where only one spin component interacts with the bath, enabling light dressing of the impurity and spin-dependent bath-impurity interactions. Through detailed comparisons with ab-initio many-body simulations, we demonstrate that the composite system is accurately described by a simplified effective Hamiltonian. The latter builds upon previously developed effective potential approaches in the absence of light dressing. It can be used to extract the impurity energy, residue, effective mass, and anharmonicity induced by the phononic dressing. Light-dressing is shown to increase the polaron residue, undressing the impurity from phononic excitations because of strong spin coupling. For strong repulsions, previously shown to trigger dynamical Bose polaron decay (a phenomenon called temporal orthogonality catastrophe), it is explained that strong light-dressing stabilizes a repulsive polaron-dressed state. Our results establish the effective Hamiltonian framework as a powerful tool for exploring strongly interacting polaronic systems and corroborating forthcoming experimental realizations.

Understanding the physics of the two-dimensional Hubbard model is widely believed to be a key step in achieving a full understanding of high-Tc cuprate superconductors. In recent years, progress has been made by large-scale numerical simulations at finite doping and, on the other hand, by microscopic theories able to capture the physics of individual charge carriers. In this work, we study single pairs of dopants in a cylindrical system using the density-matrix renormalization group algorithm. We identify two coexisting charge configurations that couple to the spin environment in different ways: a tightly bound configuration featuring (next-)nearest-neighbor pairs and a stripelike configuration of dopants on opposite sides of the cylinder, accompanied by a spin domain wall. Thus, we establish that the interplay between stripe order and uniform pairing, central to the models' phases at finite doping, has its origin at the single-pair level. By interpolating between the Hubbard and the related t-J model, we are able to quantitatively understand discrepancies in the pairing properties of the two models through the three-site hopping term usually omitted from the t-J Hamiltonian. This term is closely related to a next-nearest-neighbor tunneling t', which we observe to upset the balance between the competing stripe and pair states on the two-dopant level.

We consider the ground-state properties of a lattice Bose polaron, a quasiparticle arising from the interaction between an impurity confined to an optical lattice and a surrounding homogeneous Bose-Einstein condensate hosting phononic modes. We present an extension of Wegner's and Wilson's flow equation approach, the operator-valued-flow-equation approach, which allows us to calculate the renormalized dispersion of the polaron and assess the role of two-phonon scattering processes on the dispersion. The results obtained in this way are compared to a variational mean-field approach, which will allow us to demonstrate that the flow equations behave correctly at weak coupling. We find that in certain impurity phonon interaction regimes the shape of the dispersion is significantly altered by the inclusion of two-phonon scattering events as opposed to only single-phonon scattering events. Moreover, our results predict that a polaronic bound state may emerge, which is not present in Fr & ouml,.hlich-type models that only consider single-phonon scattering events.

One of the key features of hole-doped cuprates is the presence of an extended pseudogap phase, whose microscopic origin has been the subject of intense investigation since its discovery and is believed to be crucial for understanding high-temperature superconductivity. Various explanations have been proposed for the pseudogap, including links to symmetry-breaking orders such as stripes or pairing, and the emergence of novel fractionalized Fermi liquid (FL & lowast,.) and orthogonal metal (OM) phases. The topological nature of FL & lowast,. and OM phases has been identified as scenarios compatible with a small Fermi surface without symmetry breaking, as suggested experimentally. With recent experimental and numerical studies supporting an intricate relationship between stripe order and the pseudogap phase, we here propose an alternative scenario: an orthogonal metal with a geometric origin (GOM) driven by fluctuating domain walls. The essential mechanism behind our proposal is hidden order, where the proliferation of domain walls stabilized by charge fluctuations obscures the underlying long-range antiferromagnetic order in real space, but order is preserved in the reference frame of the background spins. As a result, well-defined fermionic quasiparticles in the form of magnetic polarons exist, which couple to Z2 topological excitations of the domain-wall string-net condensate in the ground state and constitute a small Fermi surface. At a critical doping value, we argue that hidden order is lost, driving a transition to a regular Fermi liquid at a hidden quantum critical point featuring quantum critical transport properties. Our GOM framework establishes a deep connection between the antiferromagnetic, stripe, and pseudogap phases, and suggests a possible unification of superconductivity in (electron- and hole-) doped cuprates and heavy fermion compounds.

Doping an antiferromagnetic (AFM) Mott insulator is central to our understanding of a variety of phenomena in strongly correlated electrons, including high-temperature superconductors(1,2). To describe the competition between tunnelling t of hole dopants and AFM spin interactions J, theoretical and numerical studies often focus on the paradigmatic t-J model(3) and the direct analogue quantum simulation of this model in the relevant regime of high-particle density has long been sought(4,5). Here we realize a doped quantum antiferromagnet with next-nearest-neighbour (NNN) tunnellings t (refs.(6, 7, 8, 9-10)) and hard-core bosonic holes(11) using a Rydberg tweezer platform. We use coherent dynamics between three Rydberg levels, encoding spins and holes(12), to implement a tunable bosonic t-J-V model allowing us to study previously inaccessible parameter regimes. We observe dynamical phase separation between hole and spin domains for |t/J|<< 1 and demonstrate the formation of repulsively bound hole pairs in a variety of spin backgrounds. The interference between NNN tunnellings t ' and perturbative pair tunnelling gives rise to light and heavy pairs depending on the sign of t. Using the single-site control allows us to study the dynamics of a single hole in 2D square lattice (anti)ferromagnets. The model we implement extends the toolbox of Rydberg tweezer experiments beyond spin-1/2 models(13) to a larger class of t-J and spin-1 models(14,15).

This review describes the field of Bose polarons, arising when mobile impurities are immersed into a bosonic quantum gas. The latter can be realized by a Bose-Einstein condensate of ultracold atoms, or of exciton polaritons in a semiconductor, which has led to a series of experimental observations of Bose polarons near inter-species Feshbach resonances that we survey. Following an introduction to the topic, with references to its historic roots and a presentation of the Bose polaron Hamiltonian, we summarize state-of-the-art experiments. Next we provide a detailed discussion of polaron models, starting from the ubiquitous Fr & ouml,.hlich Hamiltonian that applies at weak couplings. Already this highly simplified model allows insights into ultra-violet divergencies, logarithmic and power-law, that need to be properly regularized. To capture the physics near a Feshbach resonance, two-phonon scattering terms on the impurity as well as phonon-phonon interactions need to be included. We proceed by a survey of concurrent theoretical methods used for solving strongly interacting Bose polaron problems, ranging from Lee-Low-Pines mean-field theory, Chevy-ansatz, Gross-Pitaevskii-equation to diagrammatic Monte Carlo approaches. The subsequent sections are devoted to the large bodies of work investigating strong coupling Bose polarons, including detailed comparisons with radio-frequency spectra obtained in ultracold atom experiments,. to investigations of universal few-body and Efimov states associated with a Feshbach resonance in atomic mixtures,. to studies of quantum dynamics and polarons out of equilibrium,. Bose polarons in low-dimensional 1D and 2D quantum systems,. induced interactions among polarons and bipolaron formation,. and to Bose polarons at non-zero temperatures. We end our review by detailed discussions of closely related experimental setups and systems, including ionic impurities, systems with strong light-matter interactions, and variations and extensions of the Bose polaron concepts e.g. to baths with topological order or strong interactions relevant for correlated electrons. Finally, an outlook is presented, highlighting possible future research directions and open questions in the field as a whole.

With recent progress in quantum simulations of lattice-gauge theories, it is becoming a pressing question how to reliably protect the gauge symmetry that defines such models. Recently, an experimentally feasible gauge-protection scheme has been proposed that is based on the concept of a local pseudogenerator, which is required to act identically to the full gauge-symmetry generator in the target gauge sector, but not necessarily outside of it. The scheme has been analytically and numerically shown to reliably stabilize lattice gauge theories in the presence of perturbative errors on finite-size analog quantum-simulation devices. In this work, through uniform matrix product state calculations, we demonstrate the efficacy of this scheme for nonperturbative errors in analog quantum simulators up to all accessible evolution times in the thermodynamic limit, where it is a priori neither established nor expected that this scheme will succeed. Our results indicate the presence of an emergent gauge symmetry in an adjusted gauge theory even in the thermodynamic limit, which is beyond our analytic predictions. Additionally, we show through quantum circuit model calculations that gauge protection with local pseudogenerators also successfully suppresses gauge violations on finite quantum computers that discretize time through Trotterization. Our results firm up the robustness and feasibility of the local pseudogenerator as a viable tool for enforcing gauge invariance in modern quantum simulators and noisy intermediate-scale quantum devices.

Experimental and numerical spectroscopy have revealed rich physics in antiferromagnets, in particular in frustrated and doped systems. The Majumdar-Ghosh (MG) model has an analytically known spin-disordered ground state of dimerized singlets as a result of magnetic frustration. Here we study the single-hole angleresolved photoemission spectrum (ARPES) of a doped MG model, where we introduce a spin-hole interaction that is experimentally accessible with ultracold molecules. We report abound spinon-holon ground state and clear signatures of a spinon-holon molecule state and polarons in the ARPES spectrum at different magnetizations. Moreover, we find signatures of an emergent Feshbach resonance with tunable interactions associated with the unbinding of the spinon and the holon. Our results provide new insights into the physics of dopants in frustrated t-J models and establish the latter as a new platform for studies of emergent few-body phenomena.

The analytical study of confinement in lattice gauge theories (LGTs) remains a difficult task to this day. Taking a geometric perspective on confinement, we develop a real-space renormalization group (RG) formalism for Z(2) LGTs using percolation probability as a confinement order parameter. The RG flow we analyze is constituted by both the percolation probability and the coupling parameters. We consider a classical Z(2) LGT in two dimensions, with matter and thermal fluctuations, and analytically derive the confinement phase diagram. We find good agreement with numerical and exact benchmark results and confirm that a finite matter density enforces confinement at T < infinity in the model we consider. Our RG scheme enables future analytical studies of Z(2) LGTs with matter and quantum fluctuations and beyond.

Due to their immense representative power, neural network quantum states (NQS) have gained significant interest in current research. In recent advances in the field of NQS, it has been demonstrated that this approach can compete with state-of-the-art numerical techniques, making NQS a compelling alternative, in particular for the simulation of large, two-dimensional (2D) quantum systems. In this study, we show that recurrent neural network (RNN) wave functions can be employed to study systems relevant to current research in quantum many- body physics. Specifically, we employ a 2D tensorized gated RNN to explore the bosonic Hofstadter model with a variable local Hilbert space cutoff and long-range interactions. At first, we benchmark the RNN-NQS for the Hofstadter-Bose-Hubbard (HBH) Hamiltonian on a square lattice. We find that this method, despite the complexity of the wave function, is capable of efficiently identifying and representing most ground-state properties. Afterwards, we apply the method to an even more challenging model for current methods, namely, the Hofstadter model with long-range interactions. This model describes Rydberg-dressed atoms on a lattice subject to a synthetic magnetic field. We study systems of size up to 12 x 12 sites and identify three different regimes by tuning the interaction range and the filling fraction nu. In addition to phases known from the HBH model at short-range interaction, we observe bubble crystals and Wigner crystals for long-range interactions. Especially interesting is the evidence of a bubble crystal phase on a lattice, as this gives experiments a starting point for the search of clustered liquid phases, possibly hosting non-Abelian anyon excitations. In our work, we show that NQS are an efficient and reliable simulation method for quantum systems, which are the subject of current research. In particular, we demonstrate the ability of this method to simulate challenging systems with long-range interactions.

The search for Majorana excitations has seen tremendous efforts in recent years, ultimately aiming for their individual controllability in future topological quantum computers. A promising framework to realize such exotic Majorana fermions are topologically ordered non-Abelian phases of matter, such as certain fractional quantum Hall states. Quantum simulators provide unprecedented controllability and versatility to investigate such states, and developing experimentally feasible schemes to realize and identify them is of immediate relevance. Motivated by recent experiments, we consider bosons on coupled chains, subjected to a magnetic flux and experiencing Hubbard repulsion. At magnetic filling factor nu =1, similar systems on cylinders have been found to host the non-Abelian Moore-Read Pfaffian state in the bulk. Here, we address the question of whether more realistic few-leg ladders can host this exotic state and its chiral Majorana edge states. To this end, we perform extensive density-matrix renormalization-group simulations and determine the central charge of the ground state. While we do not find any evidence of gapless Majorana edge modes in systems of up to six legs, exact diagonalization of small systems reveals evidence for the Pfaffian state in the entanglement structure. By systematically varying the number of legs and monitoring the appearance and disappearance of this signal, our work highlights the importance of finite-size effects for the realization of exotic states in experimentally realistic systems.

Fermi surface reconstruction in cuprates can lead to an abrupt change in the Fermi momentum kF between different phases. This phenomenon remains a subject of debate and is at the heart of an ongoing discussion about the nature of the metallic state in the pseudogap regime. Here we study a minimal model of a kF changing crossover in the one-dimensional Fermi-Hubbard model, where a tuning of the on-site interaction leads to a crossover between a spin-1/2 Luttinger liquid with small Fermi momentum and a spinless chargon liquid with large Fermi momentum. We attribute this to an emergent U (1) symmetry in the strongly correlated limit, which can be used to derive a modified Luttinger sum rule recovering the large Fermi momentum. We analyze Friedel oscillations at the edge of a system to directly probe the change of Fermi momentum at zero and nonzero temperature. This paves the way for a direct experimental observation of changes of the Fermi momentum using ultracold fermions in a quantum gas microscope, with possible extensions to higher dimensional systems.

Resonant interactions associated with the emergence of a bound state constitute one of the cornerstones of modern many-body physics. Here we present a Feshbach perspective on the origin of strong pairing in Fermi-Hubbard type models. We perform a theoretical analysis of interactions between spin-polaron charge carriers in doped Mott insulators, modeled by a near-resonant two-channel scattering problem, and report evidence for Feshbach-type interactions in the dx2-y2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${d}_{{x}<^>{2}-{y}<^>{2}}$$\end{document} channel, consistent with the established phenomenology of cuprates. Existing experimental and numerical results on hole-doped cuprates lead us to conjecture the existence of a light, long-lived, low-energy excited state of two holes, which enables near-resonant interactions. To put our theory to a test we suggest to use coincidence angle-resolved photoemission spectroscopy (cARPES), pair-tunneling measurements or pump-probe experiments. The emergent Feshbach resonance among spin-polarons could also underlie superconductivity in other doped antiferromagnetic Mott insulators highlighting its potential as a unifying strong-coupling pairing mechanism rooted in quantum magnetism.

The relation between d-wave superconductivity and stripes is fundamental to the understanding of ordered phases in high-temperature cuprate superconductors1, 2, 3, 4, 5-6. These phases can be strongly influenced by anisotropic couplings, leading to higher critical temperatures, as emphasized by the recent discovery of superconductivity in nickelates7, 8, 9-10. Quantum simulators with ultracold atoms provide a versatile platform to engineer such couplings and to observe emergent structures in real space with single-particle resolution. Here we show, to our knowledge, the first signatures of individual stripes in a cold-atom Fermi-Hubbard quantum simulator using mixed-dimensional (mixD) settings. Increasing the energy scale of hole-hole attraction to the spin exchange energy, we access the interesting crossover temperature regime in which stripes begin to form11. We observe extended, attractive correlations between hole dopants and find an increased probability of forming larger structures akin to individual stripes. In the spin sector, we study correlation functions up to the third order and find results consistent with stripe formation. These observations are interpreted as a precursor to the stripe phase, which is characterized by interleaved charge and spin density wave ordering with fluctuating lines of dopants separating domains of opposite antiferromagnetic order12, 13-14.

Gauge theories constitute the basis of the Standard Model and provide useful descriptions of various phenomena in condensed matter. Realizing gauge theories on tunable tabletop quantum devices such as cold-atom quantum simulators offers the possibility to study their dynamics from first principles and to probe effects that are out of reach of dedicated particle colliders, such as deviations from gauge invariance. These quantum simulators can potentially provide insights into high-energy and nuclear physics questions, while also serving as a versatile tool for the exploration of topological phases and ergodicity-breaking mechanisms relevant to low-energy many-body physics. Recent years have seen substantial progress in the implementation of (1 + 1)D Abelian gauge theories using ultracold atoms. In this Review, we chronicle these advances, highlighting key developments in stabilizing gauge invariance and scaling up from basic building blocks to large-scale realizations where gauge-theory phenomena can be probed. We offer an outlook on future directions and the requirements for advancing this technology to the next level.

The simulation of high-temperature superconducting materials by implementing strongly correlated fermionic models in optical lattices is one of the major objectives in the field of analog quantum simulation. Here we show that local control and optical bilayer capabilities combined with spatially resolved measurements create a versatile toolbox to study fundamental properties of both nickelate and cuprate high-temperature superconductors. On the one hand, we present a scheme to implement a mixed- dimensional (mixD) bilayer model that has been proposed to capture the essential pairing physics of pressurized bilayer nickelates. This allows for the long-sought realization of a state with long-range superconducting order in current lattice quantum simulation machines. In particular, we show how coherent pairing correlations can be accessed in a partially particle-hole transformed and rotated basis. On the other hand, we demonstrate that control of local gates enables the observation of d-wave pairing order in the two-dimensional (single-layer) repulsive Fermi-Hubbard model through the simulation of a system with attractive interactions. Lastly, we introduce a scheme to measure momentum-resolved dopant densities, providing access to observables complementary to solid-state experiments-which is of particular interest for future studies of the enigmatic pseudogap phase appearing in cuprates.

We investigate the behavior of a finite-momentum impurity immersed in a weakly interacting threedimensional Bose-Einstein condensate of ultracold atoms, giving a detailed account of the dynamical quantum Cherenkov transition discussed by K. Seetharam et al. [Phys. Rev. Lett. 127, 185302 (2021)]. Using a timedependent variational approach, we identify a transition in the far-from-equilibrium dynamics of the system after the attractive short-range impurity-boson interaction is quenched on. The transition occurs as the impurity's velocity crosses an interaction-dependent critical value and manifests in the long-time behavior of the Loschmidt echo and average impurity velocity. This behavior is also reflected in the finite-momentum ground state of the system, where the group velocity of the interaction-dressed impurity loses its dependence on the total momentum of the system as the critical point is crossed. The transition we discuss should be experimentally observable via a variety of common protocols in ultracold atomic systems such as time-of-flight imaging, rf spectroscopy, Ramsey interferometry, and absorption imaging.

Lattice gauge theories (LGTs) were introduced in 1974 by Wilson to study quark confinement. These models have been shown to exhibit (de)confined phases, yet it remains challenging to define experimentally accessible order parameters. Here we propose percolation-inspired order parameters (POPs) to probe confinement of dynamical matter in Z(2) LGTs using electric field basis snapshots accessible to quantum simulators. We apply the POPs to study a classical Z(2) LGT and find a confining phase up to temperature T = infinity in two dimensions (critical T-c, i.e., finite-T phase transition, in three dimensions) for any nonzero density of Z(2) charges. Further, using quantum Monte Carlo we demonstrate that the POPs reproduce the square lattice Fradkin-Shenker phase diagram at T = infinity and explore the phase diagram at T > 0. The correlation length exponent coincides with the one of the three-dimensional Ising universality class and we determine the POP critical exponent characterizing percolation. Our proposed POPs provide a geometric perspective of confinement and are directly accessible to snapshots obtained in quantum simulators, making them suitable as a probe for quantum spin liquids.

Ever since the discovery of high-temperature superconductivity in cuprates, gaining microscopic insights into the nature of pairing in strongly correlated systems has remained one of the greatest challenges in modern condensed matter physics. Following recent experiments reporting superconductivity in the bilayer nickelate La3Ni2O7 (LNO) with remarkably high critical temperatures of Tc = 80 K, it has been argued that the low-energy physics of LNO can be described by the strongly correlated, mixed-dimensional bilayer t-J model. Here we investigate this bilayer system and utilize density matrix renormalization group techniques to establish a thorough understanding of the model and the magnetically induced pairing through comparison to the perturbative limit of dominating inter-layer spin couplings. In particular, this allows us to explain appearing finite-size effects, firmly establishing the existence of long-range superconducting order in the thermodynamic limit. By analyzing binding energies, we predict a BEC-BCS crossover as a function of the Hamiltonian parameters. We find large binding energies of the order of the inter-layer coupling that suggest strikingly high critical temperatures of the Berezinskii-Kosterlitz-Thouless transition, raising the question of whether (mixD) bilayer superconductors possibly facilitate critical temperatures above room temperature. The authors study a minimal model to describe the physics of bilayer nickelates, a novel high-temperature superconductor. They find that the model features extraordinarily high critical temperatures for superconductivity, and gain a detailed understanding of the underlying physics through an intuitive perturbative limit.

A central requirement for the faithful implementation of large-scale lattice gauge theories (LGTs) on quantum simulators is the protection of the underlying gauge symmetry. Recent advancements in the experimental realizations of large-scale LGTs have been impressive, albeit mostly restricted to Abelian gauge groups. Guided by this requirement for gauge protection, we propose an experimentally feasible approach to implement large-scale non-Abelian SU(N) and U(N) LGTs with dynamical matter in d + 1D, enabled by two-body spin-exchange interactions realizing local emergent gauge-symmetry stabilizer terms. We present two concrete proposals for 2 + 1D SU(2) and U(2) LGTs, including dynamical bosonic matter and induced plaquette terms, that can be readily implemented in current ultracold-molecule and next-generation ultracold-atom platforms. We provide numerical benchmarks showcasing experimentally accessible dynamics, and demonstrate the stability of the underlying non-Abelian gauge invariance. We develop a method to obtain the effective gauge-invariant model featuring the relevant magnetic plaquette and minimal gauge-matter coupling terms. Our approach paves the way towards near-term realizations of large-scale non-Abelian quantum link models in analog quantum simulators.

The physics of doped Mott insulators is at the heart of strongly correlated materials and is believed to constitute an essential ingredient for high-temperature superconductivity. In systems with higher SU(N) spin symmetries, even richer magnetic ground states appear at a filling of one particle per site compared to the case of SU(2) spins, but their fate upon doping remains largely unexplored. Here we address this question by studying a single hole in the SU(3) t-J model whose undoped ground state features long-range, diagonal spin stripes. By analyzing both ground-state and dynamical properties utilizing the density matrix renormalization group, we establish the appearence of magnetic polarons consisting of chargons and flavor defects whose dynamics is constrained to a single effective dimension along the ordered diagonal. We semi-analytically describe the system using geometric string theory, where paths of hole motion are the fundamental degrees of freedom. With recent advances in the realization and control of SU(N) Fermi-Hubbard models with ultracold atoms in optical lattices, our results can directly be observed in quantum gas microscopes with single-site resolution. Our work suggests the appearance of intricate ground states at finite doping constituted by emergent, coupled Luttinger liquids along diagonals, and is a first step towards exploring a wealth of physics in doped SU(N) Fermi-Hubbard models on various geometries.

A key to understanding unconventional superconductivity lies in unraveling the pairing mechanism of mobile charge carriers in doped antiferromagnets, yielding an effective attraction between charges even in the presence of strong repulsive Coulomb interactions. Here, we study pairing in a mixed-dimensional (mixD) t-J model, featuring robust binding energies-despite dominant repulsive interactions-that are strongly enhanced in the finite doping regime. The single and coupled mixD ladders we study, corresponding to bilayers of width w <= 2, feature a crossover from tightly bound pairs of holes (closed channel) at small repulsion to more spatially extended, correlated pairs of individual holes (open channel) at large repulsion. We derive an effective model for the latter, in which the attraction is mediated by the closed channel, in analogy to atomic Feshbach resonances. Using density matrix renormalization group simulations we reveal a dome of large binding energies at around 30% doping, accompanied by a change of the Fermi surface volume and a crossover from extended to tightly bound hole pairs. Our work provides a microscopic theory of pairing in the doped mixD system with dominant repulsion, closely related to bilayer, Ni-based superconductors, and our predictions can be tested in state-of-theart quantum simulators.

Microscopically understanding competing orders in strongly correlated systems is a key challenge in modern quantum many-body physics. For example, the study of magnetic polarons and their relation to pairing in the Fermi-Hubbard model in different geometries remains one of the central questions, and may help to understand the mechanism underlying unconventional superconductivity in cuprates or transition metal dichalcogenides. With recent advances in analog quantum simulation of the Fermi-Hubbard model on nonbipartite lattices, frustrated physics can now be explored experimentally in systems lacking particle-hole symmetry. Here, we study the singly doped t - J model on the triangular lattice, focusing on the competition between kinetic and magnetic frustration as a function of temperature. In doublon doped systems, we uncover a crossover between Nagaoka-type ferromagnetic (FM) correlations at high temperature and exchange mediated antiferromagnetic (AFM) order around the doublon at low temperature. For hole-doped systems, kinetic Haerter-Shastry-type AFM at high temperature as well as exchange interactions at low temperature favor 120 degrees order, strengthening magnetic correlations compared to the undoped system. In the ground state, the presence of AFM correlations throughout a wide range of interactions indicates confinement of both types of dopants. In this regime we firmly establish the presence of linear confining potentials via energy scaling arguments, supporting the picture of geometric strings in the frustrated triangular t - J model.

The combination of optical tweezer arrays with strong interactions-via dipole exchange of molecules and Van der Waals interactions of Rydberg atoms-has opened the door for the exploration of a wide variety of quantum spin models. A next significant step will be the combination of such settings with mobile dopants. This will enable one to simulate the physics believed to underlie many strongly correlated quantum materials. Here, we propose an experimental scheme to realize bosonic t-J models via encoding the local Hilbert space in a set of three internal atomic or molecular states. By engineering antiferromagnetic (AFM) couplings between spins, competition between charge motion and magnetic order similar to that in high-T-c cuprates can be realized. Since the ground states of the 2D bosonic AFM t-J model we propose to realize have not been studied extensively before, we start by analyzing the case of two dopants-the simplest instance in which their bosonic statistics plays a role-and compare our results to the fermionic case. We perform large-scale density matrix renormalization group calculations on six-legged cylinders, and find a strong tendency for bosonic holes to form stripes. This demonstrates that bosonic, AFM t-J models may contain similar physics as the collective phases in strongly correlated electrons.

Confinement is a paradigmatic phenomenon of gauge theories, and its understanding lies at the forefront of high-energy physics. Here, we study confinement in a simple one-dimensional 7L2 lattice gauge theory at finite temperature and filling, which is within the reach of current cold-atom and superconducting-qubit platforms. By employing matrix product states (MPS) calculations, we investigate the decay of the finite-temperature Green's function and uncover a smooth crossover between the confined and deconfined regimes. Furthermore, using the Friedel oscillations and string length distributions obtained from snapshots sampled from MPS, both of which are experimentally readily available, we verify that confined mesons remain well-defined at arbitrary finite temperature. This phenomenology is further supported by probing quench dynamics of mesons with exact diagonalization. Our results shed new light on confinement at finite temperature from an experimentally relevant standpoint.

Topological phase transitions go beyond Ginzburg and Landau's paradigm of spontaneous symmetry breaking and occur without an associated local order parameter. Instead, such transitions can be characterized by the emergence of nonlocal order parameters, which require measurements on extensively many particles simultaneously-an impossible venture in real materials. On the other hand, quantum simulators have demonstrated such measurements, making them prime candidates for experimental confirmation of nonlocal topological order. Here, building upon the recent advances in preparing few-particle fractional Chern insulators using ultracold atoms and photons, we propose a realistic scheme for detecting the hidden off-diagonal long-range order (HODLRO) characterizing Laughlin states. Furthermore, we demonstrate the existence of this hidden order in fractional Chern insulators, specifically for the nu = 1/2-Laughlin state in the isotropic Hofstadter-Bose-Hubbard model. This is achieved by large-scale numerical density matrix renormalization group (DMRG) simulations based on matrix product states, for which we formulate an efficient sampling procedure providing direct access to HODLRO in close analogy to the proposed experimental scheme. We confirm the characteristic power-law scaling of HODLRO, with an exponent 1/nu = 2, and show that its detection requires only a few thousand snapshots. This makes our scheme realistically achievable with current technology and paves the way for further analysis of nonlocal topological orders, e.g., in topological states with non-Abelian anyonic excitations.

"When a mobile hole is doped into an antiferromagnet, its movement will distort the surrounding magnetic order and yield a magnetic polaron. The resulting complex interplay of spin and charge degrees of freedom gives rise to very rich physics and is widely believed to be at the heart of high -temperature superconductivity in cuprates. In this paper, we develop a quantitative theoretical formalism, based on the phenomenological parton description, to describe magnetic polarons in the strong -coupling regime. We construct an effective Hamiltonian with weak coupling to the spin -wave excitations in the background, making the use of standard polaronic methods possible. Our starting point is a single hole doped into an antiferromagnet described by a ""geometric string"" capturing the strongly correlated hopping processes of charge and spin degrees of freedom, beyond linear spin -wave approximation. Subsequently, we introduce magnon excitations through a generalized 1 / S expansion and derive an effective coupling of these spin waves to the hole plus the string (the meson) to arrive at an effective polaron Hamiltonian with density -density type interactions. After making a Born-Oppenheimer-type approximation, this system is solved using the self -consistent Born approximation to extract the renormalized polaron properties. We apply our formalism (i) to calculate beyond linear spin -wave angle -resolved photoemission spectroscopy spectra, (ii) to reveal the interplay of rovibrational meson excitations, and (iii) to analyze the pseudogap expected at low doping. Moreover, our work paves the way for exploring magnetic polarons out of equilibrium or in frustrated systems, where weak -coupling approaches are desirable and going beyond linear spin -wave theory becomes necessary."

Modeling the underlying pairing mechanism of charge carriers in strongly correlated electrons, starting from a microscopic theory, is among the central challenges of condensed -matter physics. Hereby, the key task is to understand what causes the appearance of superconductivity at comparatively high temperatures upon hole doping an antiferromagnetic (AFM) Mott insulator. Recently, it has been proposed that at strong coupling and low doping, the fundamental one- and two -hole meson -type constituents-magnetic polarons and bipolaronic pairs-likely realize an emergent Feshbach resonance producing near -resonant d x 2 - y 2 interactions between charge carriers. Here, we provide detailed calculations of the proposed scenario by describing the open and closed meson scattering channels in the t - t ' - J model using a truncated basis method. After integrating out the closed channel constituted by bipolaronic pairs, we find d x 2 - y 2 attractive interactions between open channel magnetic polarons. The closed form of the derived interactions allows us analyze the resonant pairing interactions and we find enhanced (suppressed) attraction for hole (electron) doping in our model. The formalism we introduce provides a framework to analyze the implications of a possible Feshbach scenario, e.g., in the context of BEC-BCS crossover, and establishes a foundation to test quantitative aspects of the proposed Feshbach pairing mechanisms in doped antiferromagnets.

The study of fractional Chern insulators and their exotic anyonic excitations poses a major challenge in current experimental and theoretical research. Quantum simulators, in particular ultracold atoms in optical lattices, provide a promising platform to realize, manipulate, and understand such systems with a high degree of controllability. Recently, an atomic nu = 1/2 Laughlin state has been realized experimentally for a small system of two particles on 4 x 4 sites [Leonard et al., Nature (London) 619, 495 (2023)]. The next challenge concerns the preparation of Laughlin states in extended systems, ultimately giving access to anyonic braiding statistics or gapless chiral edge-states in systems with open boundaries. Here, we propose and analyze an experimentally feasible scheme to grow larger Laughlin states by connecting multiple copies of the already-existing 4 x 4 system. First, we present a minimal setting obtained by coupling two of such patches, producing an extended 8 x 4 system with four particles. Then, we analyze different preparation schemes, setting the focus on two shapes for the extended system, and discuss their respective advantages: While growing striplike lattices could give experimental access to the central charge, squarelike geometries are advantageous for creating quasihole excitations in view of braiding protocols. We highlight the robust quantization of the fractional quasihole charge upon using our preparation protocol. We benchmark the performance of our patchwork preparation scheme by comparing it to a protocol based on coupling one-dimensional chains. We find that the patchwork approach consistently gives higher target-state fidelities, especially for elongated systems. The results presented here pave the way towards near-term implementations of extended Laughlin states in quantum gas microscopes and the subsequent exploration of exotic properties of topologically ordered systems in experiments.

Since the discovery of superconductivity in cuprate materials, the minimal ingredients for high-Tc superconductivity have been an outstanding puzzle. Motivated by the recently discovered nickelate bilayer superconductor La3Ni2O7 under pressure, we study a minimal bilayer model, in which, as in La3Ni2O7, interlayer and intralayer magnetic interactions but no interlayer hopping are present: A mixed-dimensional (mixD) t-J model. In the setting of a mixD ladder, we show that the system exhibits a crossover associated with a Feshbach resonance: From a closed-channel-dominated regime of tightly bound bosonic pairs of holes to an open-channel-dominated regime of spatially more extended Cooper pairs. The crossover can be tuned by varying doping, or by a nearest-neighbor Coulomb repulsion V that we include in our model. Using density matrix renormalization group simulations and analytical descriptions of both regimes, we find that the ground state is a Luther-Emery liquid, competing with a density wave of tetraparton plaquettes at commensurate filling delta = 0.5 at large repulsion, and exhibits a pairing dome where binding is facilitated by doping. Our observations can be understood in terms of pairs of correlated spinon-chargon excitations constituting the open channel, which are subject to attractive interactions mediated by the closed channel of tightly bound chargon-chargon pairs. When the closed channel is lowered in energy by doping or tuning V, a Feshbach resonance is realized, associated with a dome in the binding energy. Our predictions can be directly tested in state-of-the art quantum simulators, and we argue that the pairing mechanism we describe may be realized in the nickelate bilayer superconductor La3Ni2O7.

A key step in unraveling the mysteries of materials exhibiting unconventional superconductivity is to understand the underlying pairing mechanism. While it is widely agreed upon that the pairing glue in many of these systems originates from antiferromagnetic spin correlations, a microscopic description of pairs of charge carriers remains lacking. Here we use state-of-the art numerical methods to probe the internal structure and dynamical properties of pairs of charge carriers in quantum antiferromagnets in four-legged cylinders. Exploiting the full momentum resolution in our simulations, we are able to distinguish two qualitatively different types of bound states: a highly mobile, meta-stable pair, which has a dispersion proportional to the hole hopping t, and a heavy pair, which can only move due to spin exchange processes and turns into a flat band in the Ising limit of the model. Understanding the pairing mechanism can on the one hand pave the way to boosting binding energies in related models, and on the other hand enable insights into the intricate competition of various phases of matter in strongly correlated electron systems.

Recently, tremendous progress has been made in the field of quantum science and technologies: different platforms for quantum simulation as well as quantum computing, ranging from superconduct-ing qubits to neutral atoms, are start-ing to reach unprecedentedly large sys-tems. In order to benchmark these sys-tems and gain physical insights, the need for efficient tools to characterize quantum states arises. The exponential growth of the Hilbert space with system size ren-ders a full reconstruction of the quantum state prohibitively demanding in terms of the number of necessary measurements. Here we propose and implement an ef-ficient scheme for quantum state tomog-raphy using active learning. Based on a few initial measurements, the active learn-ing protocol proposes the next measure-ment basis, designed to yield the max-imum information gain. We apply the active learning quantum state tomogra-phy scheme to reconstruct different multi-qubit states with varying degree of entan-glement as well as to ground states of the XXZ model in 1D and a kinetically con-strained spin chain. In all cases, we obtain a significantly improved reconstruction as compared to a reconstruction based on the exact same number of measurements and measurement configurations, but with ran-domly chosen basis configurations. Our scheme is highly relevant to gain physical insights in quantum many-body systems as well as for benchmarking and character-izing quantum devices, e.g. for quantum simulation, and paves the way for scalable adaptive protocols to probe, prepare, and manipulate quantum systems.

Unveiling the microscopic origins of quantum phases dominated by the interplay of spin and motional degrees of freedom constitutes one of the central challenges in strongly correlated many-body physics. When holes move through an antiferromagnetic spin background, they displace the positions of spins, which induces effective frustration in the magnetic environment. However, a concrete characterization of this effect in a quantum many-body system is still an unsolved problem. Here we present a Hamiltonian reconstruction scheme that allows for a precise quantification of hole-motion-induced frustration. We access non-local correlation functions through projective measurements of the many-body state, from which effective spin-Hamiltonians can be recovered after detaching the magnetic background from dominant charge fluctuations. The scheme is applied to systems of mixed dimensionality, where holes are restricted to move in one dimension, but SU(2) superexchange is two-dimensional. We demonstrate that hole motion drives the spin background into a highly frustrated regime, which can quantitatively be described by an effective J(1)-J(2)-type spin model. We exemplify the applicability of the reconstruction scheme to ultracold atom experiments by recovering effective spin-Hamiltonians of experimentally obtained 1D Fermi-Hubbard snapshots. Our method can be generalized to fully 2D systems, enabling promising microscopic perspectives on the doped Hubbard model.

Strongly interacting topological matter(1) exhibits fundamentally new phenomena with potential applications in quantum information technology(2,3). Emblematic instances are fractional quantum Hall (FQH) states(4), in which the interplay of a magnetic field and strong interactions gives rise to fractionally charged quasi-particles, long-ranged entanglement and anyonic exchange statistics. Progress in engineering synthetic magnetic fields(5-21) has raised the hope to create these exotic states in controlled quantum systems. However, except for a recent Laughlin state of light(22), preparing FQH states in engineered systems remains elusive. Here we realize a FQH state with ultracold atoms in an optical lattice. The state is a lattice version of a bosonic ? = 1/2 Laughlin state(4,23) with two particles on 16 sites. This minimal system already captures many hallmark features of Laughlin-type FQH states(24-28): we observe a suppression of two-body interactions, we find a distinctive vortex structure in the density correlations and we measure a fractional Hall conductivity of s(H)/s(0) = 0.6(2) by means of the bulk response to a magnetic perturbation. Furthermore, by tuning the magnetic field, we map out the transition point between the normal and the FQH regime through a spectroscopic investigation of the many-body gap. Our work provides a starting point for exploring highly entangled topological matter with ultracold atoms(29-33).

Gauge fields coupled to dynamical matter are ubiquitous in many disciplines of physics, ranging from particle to condensed matter physics, but their implementation in large-scale quantum simulators remains challenging. Here we propose a realistic scheme for Rydberg atom array experiments in which a Z2 gauge structure with dynamical charges emerges on experimentally relevant timescales from only local two-body interactions and one-body terms in two spatial dimensions. The scheme enables the experimental study of a variety of models, including (2+ 1)D Z2 lattice gauge theories coupled to different types of dynamical matter and quantum dimer models on the honeycomb lattice, for which we derive effective Hamiltonians. We discuss ground-state phase diagrams of the experimentally most relevant effective Z2 lattice gauge theories with dynamical matter featuring various confined and deconfined, quantum spin liquid phases. Further, we present selected probes with immediate experimental relevance, including signatures of disorder-free localization and a thermal deconfinement transition of two charges.

In strongly correlated quantum materials, the behavior of charge carriers is dominated by strong electron-electron interactions. These can lead to insulating states with spin order, and upon doping to competing ordered states including unconventional super-conductivity. The underlying pairing mechanism remains poorly understood however, even in strongly simplified theoretical models. Recent advances in quantum simulation allow to study pairing in paradigmatic settings, e.g. in the t - J and t - Jz Hamiltoni-ans. Even there, the most basic properties of paired states of only two dopants, such as their dispersion relation and excitation spectra, remain poorly studied in many cases. Here we provide new analytical insights into a possible string-based pairing mechanism of mobile holes in an antiferromagnet. We analyze an effective model of partons con-nected by a confining string and calculate the spectral properties of bound states. Our model is equally relevant for understanding Hubbard-Mott excitons consisting of a bound doublon-hole pair or confined states of dynamical matter in lattice gauge theories, which motivates our study of different parton statistics. Although an accurate semi-analytic es-timation of binding energies is challenging, our theory provides a detailed understanding of the internal structure of pairs. For example, in a range of settings we predict heavy states of immobile pairs with flat-band dispersions - including for the lowest-energy d -wave pair of fermions. Our findings shed new light on the long-standing question about the origin of pairing and competing orders in high-temperature superconductors.

Quantum many-body scarring is a paradigm of weak ergodicity breaking arising due to the presence of special nonthermal many-body eigenstates that possess low entanglement entropy, are equally spaced in energy, and concentrate in certain parts of the Hilbert space. Though scars have been shown to be intimately connected to gauge theories, their stability in such experimentally relevant models is still an open question, and it is generally considered that they exist only under fine-tuned conditions. In this work, we show through Krylov-based time-evolution methods how quantum many-body scars can be made robust in the presence of experimental errors through utilizing terms linear in the gaugesymmetry generator or a simplified pseudogenerator in U(1) and Z2 lattice gauge theories. Our findings are explained by the concept of quantum Zeno dynamics. Our experimentally feasible methods can be readily implemented in existing large-scale ultracold-atom quantum simulators and setups of Rydberg atoms with optical tweezers.

Microscopically understanding competing orders in strongly correlated systems is a key challenge in modern quantum many-body physics. For example, the origin of stripe order and its relation to pairing in the Fermi -Hubbard model remains one of the central questions, and may help to understand the origin of high-temperature superconductivity in cuprates. Here, we analyze stripe formation in the doped mixed-dimensional (mixD) variant of the t - J model, where charge carriers are restricted to move only in one direction, whereas magnetic SU(2) interactions are two-dimensional. Using the density matrix renormalization group at finite temperature, we find a stable vertical stripe phase in the absence of pairing, featuring incommensurate magnetic order and long-range charge density wave profiles over a wide range of dopings. We find high critical temperatures on the order of the magnetic coupling similar to J/2, hence being within reach of current quantum simulators. Snapshots of the many-body state, accessible to quantum simulators, reveal hidden spin correlations in the mixD setting, whereby antiferromagnetic correlations are enhanced when considering purely the magnetic background. The proposed model can be viewed as realizing a parent Hamiltonian of the stripe phase, whose hidden spin correlations lead to the predicted resilience against quantum and thermal fluctuations.

We study the phase diagram of two different mixed-dimensional t-Jz-J1 models on the square lattice, in which the hopping amplitude t is only nonzero along the x direction. In the first model, which is bosonic, the spin-exchange amplitude J1 is negative and isotropic along the x and y directions of the lattice, and Jz is isotropic and positive. The low-energy physics is characterized by spin-charge separation: the holes hop as free fermions in an easy-plane ferromagnetic background. In the second model, J1 is restricted to the x axis while Jz remains isotropic and positive. The model is agnostic to particle statistics, and shows stripe patterns with antiferromagnetic Neel order at low temperature and high hole densities, in resemblance of the mixed-dimensional t-Jz and t -J models. At lower hole density, a very strong first-order transition and hysteresis loop is seen extending to a remarkably high 14(1)% hole doping.

Strongly interacting fermionic systems host a variety of interesting quantum many-body states with exotic excitations. For instance, the interplay of strong interactions and the Pauli exclusion principle can lead to Stoner ferromagnetism, but the fate of this state remains unclear when kinetic terms are added. While in many lattice models the fermions' dispersion results in delocalization and destabilization of the ferromagnet, flat bands can restore strong interaction effects and ferromagnetic correlations. To reveal this interplay, here we propose to study the Hofstadter-Fermi-Hubbard model using ultracold atoms. We demonstrate, by performing large-scale density-matrix renormalization group simulations, that this model exhibits a lattice analog of the quantum Hall (QH) ferromagnet at magnetic filling factor nu = 1. We reveal the nature of the low energy spin-singlet states around nu asymptotic to 1 and find that they host quasi-particles and quasi-holes exhibiting spin-spin correlations reminiscent of skyrmions. Finally, we predict the breakdown of flat-band ferromagnetism at large fields. Our work paves the way towards experimental studies of lattice QH ferromagnetism, including prospects to study many-body states of interacting skyrmions and explore the relation to high- T-c superconductivity.

It is a widely accepted view that the interplay of spin and charge degrees of freedom in doped antiferromagnets (AFMs) gives rise to the rich physics of high-temperature superconductors. Nevertheless, it remains unclear how effective low-energy degrees of freedom and the corresponding field theories emerge from microscopic models, including t - J and Hubbard Hamiltonians. A promising view comprises that the charge carriers have a rich internal parton structure on intermediate scales, but the interplay of the emergent partons with collective magnon excitations of the surrounding AFM remains unexplored. Here we study a doped one-dimensional spin chain in a staggered magnetic field and demonstrate that it supports a zoo of various long-lived excitations. These include magnons, mesonic pairs of spinons and chargons along with their rovibrational excitations, and tetraparton bound states of mesons and magnons. We identify these types of quasiparticles in various spectra using density-matrix renormalization group simulations. Moreover, we introduce a strong-coupling theory describing the polaronic dressing and molecular binding of mesons to collective magnon excitations. The effective theory can be solved by standard tools developed for polaronic problems and can be extended to study similar physics in two-dimensional doped AFMs in the future. Experimentally, the doped spin-chain in a staggered field can be directly realized in quantum gas microscopes.

Coupling dynamical charges to gauge fields can result in highly non-local interactions with a linear confining potential. As a consequence, individual particles bind into mesons which, in one dimension, become the new constituents of emergent Luttinger liquids (LLs). Furthermore, at commensurate fillings, different Mott-insulating states can be stabilized by including nearest-neighbour (NN) interactions among charges. However, rich phase diagrams expected in such models have not been fully explored and still lack comprehensive theoretical explanation. Here, by combining numerical and analytical tools, we study a simple one-dimensional Z2 lattice gauge theory at half-filling, where U(1) matter is coupled to gauge fields and interacts through NN repulsion. We uncover a rich phase diagram where the local NN interaction stabilizes a Mott state of individual charges (or partons) on the one hand, and an LL of confined mesons on the other. Furthermore, at the interface between these two phases, we uncover a highly frustrated regime arising due to the competition between the local NN repulsion and the non-local confining interactions, realizing a pre-formed parton plasma. Our work is motivated by the recent progress in ultracold atom experiments, where such simple model could be readily implemented. For this reason we calculate the static structure factor which we propose as a simple probe to explore the phase diagram in an experimental setup.

Conventional superconductivity emerges from pairing of charge carriers-electrons or holes-mediated by phonons(1). In many unconventional superconductors, the pairing mechanism is conjectured to be mediated by magnetic correlations(2), as captured by models of mobile charges in doped antiferromagnets(3). However, a precise understanding of the underlying mechanism in real materials is still lacking and has been driving experimental and theoretical research for the past 40 years. Early theoretical studies predicted magnetic-mediated pairing of dopants in ladder systems(4-8), in which idealized theoretical toy models explained how pairing can emerge despite repulsive interactions(9). Here we experimentally observe this long-standing theoretical prediction, reporting hole pairing due to magnetic correlations in a quantum gas of ultracold atoms. By engineering doped antiferromagnetic ladders with mixed-dimensional couplings(10), we suppress Pauli blocking of holes at short length scales. This results in a marked increase in binding energy and decrease in pair size, enabling us to observe pairs of holes predominantly occupying the same rung of the ladder. We find a hole-hole binding energy of the order of the superexchange energy and, upon increased doping, we observe spatial structures in the pair distribution, indicating repulsion between bound hole pairs. By engineering a configuration in which binding is strongly enhanced, we delineate a strategy to increase the critical temperature for superconductivity.

Disorder-free localization in translation-invariant gauge theories presents a counterintuitive yet powerful framework of ergodicity breaking in quantum many-body physics. The fragility of this phenomenon in the presence of gauge-breaking errors has recently been addressed, but no scheme has been able to reliably stabilize disorder-free localization through all accessible evolution times while preserving the disorder-free property. Here, we introduce the concept of Stark gauge protection, which entails a linear sum in gauge-symmetry local (pseudo)generators weighted by a Stark potential. Using exact diagonalization and Krylov-based methods, we show how this scheme can stabilize or even enhance disorder-free localization against gauge-breaking errors in U(1) and Z2 gauge theories up to all accessible evolution times, without inducing bona fide Stark many-body localization. We show through a Magnus expansion that the dynamics under Stark gauge protection is described by an effective Hamiltonian where gauge-breaking terms are suppressed locally by the protection strength and additionally by the matter site index, which we argue is the main reason behind stabilizing the localization up to all accessible times. Our scheme is readily feasible in modern ultracold-atom experiments and Rydberg-atom setups with optical tweezers.

Synthetic lattice systems are powerful platforms for studying the influence of intrinsic nonlinearities on topological phenomena. Here the authors elucidate the topological transport of solitons in terms of Wannier functions displacement and they introduce a nonlinearity-induced topological transport effect that could be observed in ultracold quantum mixtures. Exploring the interplay between topological band structures and tunable nonlinearities has become possible with the development of synthetic lattice systems. In this emerging field of nonlinear topological physics, an experiment revealed the quantized motion of solitons in Thouless pumps and suggested that this phenomenon was dictated by the Chern number of the band from which solitons emanate. Here, we elucidate the origin of this nonlinear topological effect, by showing that the motion of solitons is established by the quantized displacement of the underlying Wannier functions. Our general theoretical approach, which fully clarifies the central role of the Chern number in solitonic pumps, provides a framework for describing the topological transport of nonlinear excitations in a broad class of physical systems. Exploiting this interdisciplinarity, we introduce an interaction-induced topological pump for ultracold atomic mixtures, where solitons of impurity atoms experience a quantized drift resulting from genuine interaction processes with their environment.

We investigate the formation of magnetic Bose polaron, an impurity atom dressed by spin-wave excitations, in a one-dimensional spinor Bose gas. Within an effective potential model, the impurity is strongly confined by the host excitations which can even overcome the impurity-medium repulsion leading to a self-localized quasi-particle state. The phase diagram of the attractive and self-bound repulsive magnetic polaron, repulsive non-magnetic (Frohlich-type) polaron and impurity-medium phase-separation regimes is explored with respect to the Rabi-coupling between the spin components, spin-spin interactions and impurity-medium coupling. The residue of such magnetic polarons decreases substantially in both strong attractive and repulsive branches with strong impurity-spin interactions, illustrating significant dressing of the impurity. The impurity can be used to probe and maneuver the spin polarization of the magnetic medium while suppressing ferromagnetic spin-spin correlations. It is shown that mean-field theory fails as the spinor gas approaches immiscibility since the generated spin-wave excitations are prominent. Our findings illustrate that impurities can be utilized to generate controllable spin-spin correlations and magnetic polaron states which can be realized with current cold atom setups.

The postulate of gauge invariance in nature does not lend itself directly to implementations of lattice gauge theories in modern setups of quantum synthetic matter. Unavoidable gauge-breaking errors in such devices require gauge invariance to be enforced for faithful quantum simulation of gauge-theory physics. This poses major experimental challenges, in large part due to the complexity of the gauge-symmetry generators. Here, we show that gauge invariance can be reliably stabilized by employing simplified local pseudogenerators designed such that within the physical sector they act identically to the actual local generator. Dynamically, they give rise to emergent exact gauge theories up to time scales polynomial and even exponential in the protection strength. This obviates the need for implementing often complex multibody full gauge symmetries, thereby further reducing experimental overhead in physical realizations. We showcase our method in the Z(2) lattice gauge theory, and discuss experimental considerations for its realization in modern ultracold-atom setups.

Experimental realizations of topologically ordered states of matter, such as fractional quantum Hall states, with cold atoms are now within reach. In particular, optical lattices provide a promising platform for the realization and characterization of such states, where novel detection schemes enable an unprecedented microscopic under-standing. Here we show that the central charge can be directly measured in current cold atom experiments using the number entropy as a proxy for the entanglement entropy. We perform density-matrix renormalization-group simulations of Hubbard-interacting bosons on coupled chains subject to a magnetic field with alpha = 1/4 flux quanta per plaquette. Tuning the interchain hopping, we find a transition from a trivial quasi-one-dimensional phase to the topologically ordered Laughlin state at magnetic filling factor nu = 1/2 for systems of three or more chains. We resolve the transition using the central charge, on-site correlations, momentum distributions, and the many-body Chern number. Additionally, we propose a scheme to experimentally estimate the central charge from Fock basis snapshots. The model studied here is experimentally realizable with existing cold atom techniques and the proposed observables pave the way for the detection and classification of a larger class of interacting topological states of matter.

In strongly interacting systems with a separation of energy scales, low-energy effective Hamiltonians help provide insights into the relevant physics at low temperatures. The emergent interactions in the effective model are mediated by virtual excitations of high-energy states: For example, virtual doublon-hole excitations in the Fermi-Hubbard model mediate antiferromagnetic spin-exchange interactions in the derived effective model, known as the t - J - 3s model. Formally this procedure is described by performing a unitary Schrieffer-Wolff basis transformation. In the context of quantum simulation, it can be advantageous to consider the effective model to interpret experimental results. However, virtual excitations such as doublon-hole pairs can obfuscate the measurement of physical observables. Here we show that quantum simulators allow one to access the effective model even more directly by performing measurements in a rotated basis. We propose a protocol to perform a Schrieffer-Wolff transformation on Fermi-Hubbard low-energy eigenstates (or thermal states) to dynamically prepare approximate t - J - 3s model states using fermionic atoms in an optical lattice. Our protocol involves performing a linear ramp of the optical lattice depth, which is slow enough to eliminate the virtual doublon-hole fluctuations but fast enough to freeze out the dynamics in the effective model. We perform a numerical study using exact diagonalization and find an optimal ramp speed for which the state after the lattice ramp has maximal overlap with the t - J - 3s model state. We compare our numerics to experimental data from our Lithium-6 fermionic quantum gas microscope and show a proof-of-principle demonstration of this protocol. More generally, this protocol can be beneficial to studies of effective models by enabling the suppression of virtual excitations in a wide range of quantum simulation experiments.

The mechanism underlying charge transport in strongly correlated quantum systems, such as doped antiferromagnetic Mott insulators, remains poorly understood. Here, we study the expansion dynamics of an initially localized hole inside a two-dimensional (2D) Ising antiferromagnet at variable temperature. Using a combination of classical Monte Carlo and truncated-basis methods, we reveal two dynamically distinct regimes: a spin-charge confined region below a critical temperature T*, characterized by slow spreading, and a spin-charge deconfined region above T*, characterized by an unbounded diffusive expansion. The deconfinement temperature T* ti 0.65Jz we find is around the N??el temperature TN = 0.567Jz of the Ising background in 2D, but we expect T* < TN in higher dimensions. In both regimes we find that the mobile hole does not thermalize with the Ising spin background on the considered time scales, indicating weak effective coupling of spin and charge degrees of freedom. Our results can be qualitatively understood by an effective parton model and can be tested experimentally in state-of-the-art quantum gas microscopes.

Realizing and characterizing interacting topological phases in synthetic quantum systems is a formidable challenge. Here, we propose a Floquet protocol to realize the antiferromagnetic Heisenberg model with power -law decaying interactions. Based on analytical and numerical arguments, we show that this model features a quantum phase transition from a liquid to a valence bond solid that spontaneously breaks lattice translational symmetry and is reminiscent of the Majumdar-Ghosh state. The different phases can be probed dynamically by measuring the evolution of a fully dimerized state. We moreover introduce an interferometric protocol to characterize the topological excitations and the bulk topological invariants of the interacting many-body system.

Studies of unconventional pairing mechanisms in cold atoms require ultralow temperatures. Large-scale numerics show that certain bilayer models allow for deeply bound and highly mobile pairs of charges at more accessible temperatures. Interacting many-body systems in reduced-dimensional settings, such as ladders and few-layer systems, are characterized by enhanced quantum fluctuations. Recently, two-dimensional bilayer systems have sparked considerable interest because they can host unusual phases, including unconventional superconductivity. Here we present a theoretical proposal for realizing high-temperature pairing of fermions in a class of bilayer Hubbard models. We introduce a general and highly efficient pairing mechanism for mobile charge carriers in doped antiferromagnetic Mott insulators. The pairing is caused by the energy that one charge gains when it follows the path created by another charge. We show that this mechanism leads to the formation of highly mobile but tightly bound pairs in the case of mixed-dimensional Fermi-Hubbard bilayer systems. This setting is closely related to the Fermi-Hubbard model believed to capture the physics of copper oxides, and can be realized in currently available ultracold atom experiments.

The bulk-boundary correspondence relates quantized edge states to bulk topological invariants in topological phases of matter. In one-dimensional symmetry-protected topological systems, quantized topological Thouless pumps directly reveal this principle and provide a sound mathematical foundation. Symmetry-protected higher-order topological phases of matter (HOSPTs) also feature a bulk-boundary correspondence, but its connection to quantized charge transport remains elusive. Here, we show that quantized Thouless pumps connecting C-4-symmetric HOSPTs can be described by a tuple of four Chern numbers that measure quantized bulk charge transport in a direction-dependent fashion. Moreover, this tuple of Chern numbers allows to predict the sign and value of fractional corner charges in the HOSPTs. We show that the topologically nontrivial phase can be characterized by both quadrupole and dipole configurations, shedding new light on current debates about the multipole nature of the HOSPT bulk. By employing corner-periodic boundary conditions, we generalize Restas's theory to HOSPTs. Our approach provides a simple framework for understanding topological invariants of general HOSPTs and paves the way for an in-depth description of future dynamical experiments.

Disorder-free localization is a recently discovered phenomenon of nonergodicity that can emerge in quantum many-body systems hosting gauge symmetries when the initial state is prepared in a superposition of gauge superselection sectors. Thermalization is then prevented up to all accessible evolution times despite the model being nonintegrable and translation invariant. In a recent work [Halimeh et al., arXiv:2111.02427 (2021)], it has been shown that terms linear in the gauge-symmetry generator stabilize disorder-free localization in U(1) gauge theories against gauge errors that couple different superselection sectors. Here, we show in the case of Z2 gauge theories that disorder-free localization can not only be stabilized, but also enhanced by the addition of translation-invariant terms linear in a local Z2 pseudogenerator that acts identically to the full generator in a single superselection sector, but not necessarily outside of it. We show analytically and numerically how this leads through the quantum Zeno effect to the dynamical emergence of a renormalized gauge theory with an enhanced local symmetry, which contains the Z2 gauge symmetry of the ideal model, associated with the Z2 pseudogenerator. The resulting proliferation of superselection sectors due to this dynamically emergent gauge theory creates an effective disorder greater than that in the original model, thereby enhancing disorder-free localization. We demonstrate the experimental feasibility of the Z2 pseudogenerator by providing a detailed readily implementable experimental proposal for the observation of disorder-free localization in a Rydberg setup.

Quantum simulators are attracting great interest because they promise insight into the behavior of quantum many-body systems that are prohibitive for classical simulations. The generic output of quantum simulators are snapshots, obtained by means of projective measurements. These provide new information, such as full distribution functions, that goes beyond the more commonly evaluated expectation values of observables while adding shot-noise uncertainty to the latter. Hence, a central goal of theoretical efforts must be to predict these exact same quantities that can be measured in experiments. Here, we report on a snapshot-based study of particle currents in quantum lattice models with a conserved number of particles. It is shown how the full probability distribution of locally resolved particle currents can be obtained from suitable snapshot data. Moreover, we investigate the Hall response of interacting bosonic flux ladders, exploiting snapshots drawn from matrix-product states. Flux ladders are minimal lattice models, which enable microscopic studies of the Hall response in correlated quantum phases, and they are successfully realized in current quantum-gas experiments. Using a specific pattern of unitary two-site transformations, it is shown that the Hall polarization and the Hall voltage can be faithfully computed from a realistic number of snapshots obtained in experimentally feasible quench and finite-bias simulations.

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.

The challenge of understanding the dynamics of a mobile impurity in an interacting quantum many-body medium comes from the necessity of including entanglement between the impurity and excited states of the environment in a wide range of energy scales. In this Letter, we investigate the motion of a finite mass impurity injected into a three-dimensional quantum Bose fluid as it starts shedding Bogoliubov excitations. We uncover a transition in the dynamics as the impurity's velocity crosses a critical value that depends on the strength of the interaction between the impurity and bosons as well as the impurity's recoil energy. We find that in injection experiments, the two regimes differ not only in the character of the impurity velocity abatement but also exhibit qualitative differences in the Loschmidt echo, density ripples excited in the Bose-Einstein condensate, and momentum distribution of scattered bosonic particles. The transition is a manifestation of a dynamical quantum Cherenkov effect and should be experimentally observable with ultracold atoms using Ramsey interferometry, rf spectroscopy, absorption imaging, and time-of-flight imaging.

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.

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.

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.

Kitaev's toric code is an exactly solvable model with Z(2)-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 Z(2) 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 nontrivial 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 wave function and introduce a minimal braiding protocol. Detecting a p-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.

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.

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 v = 1 and find strong indications that at alpha = 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 v = 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.

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.

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.

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 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.

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 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.

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 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 C-4 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.

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.

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.

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.

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 he 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 he utilized to implement the strongly frustrated J-Q model, a candidate for hosting a deconfined quantum critical point.

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.

Quantum simulation has the potential to investigate gauge theories in strongly interacting regimes, which are currently inaccessible through conventional numerical techniques. Here, we take a first step in this direction by implementing a Floquet-based method for studying Z(2) I 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 Z(2) I 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.

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 Z(2) 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 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.

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.

The interplay of spin and charge degrees of freedom, introduced by doping mobile holes into a Mott insulator with strong antiferromagnetic (AFM) correlations, is at the heart of strongly correlated matter such as high-T-C cuprate superconductors. Here, we capture this interplay in the strong coupling regime and propose a trial wave function of mobile holes in an AFM. Our method provides a microscopic justification for a class of theories which describe doped holes moving in an AFM environment as mesonlike bound states of spinons and chargons. We discuss a model of such bound states from the perspective of geometric strings, which describe a fluctuating lattice geometry introduced by the fast motion of the chargon, relative to the spinon. This is demonstrated to give rise to short-range hidden string order, signatures of which have recently been revealed by ultracold atom experiments at elevated temperatures. We present evidence for such short-range hidden string correlations also at zero temperature by performing numerical density-matrix renormalization-group simulations. To test our microscopic approach, we calculate the ground-state energy and dispersion relation of a hole in an AFM, as well as the magnetic polaron radius, and obtain good quantitative agreement with advanced numerical simulations at strong coupling. We discuss extensions of our analysis to systems without long-range AFM order to systems with short-range magnetic correlations.