Topological phase transitions challenge conventional paradigms in many-body physics by separating phases that are locally indistinguishable yet globally distinct. Using a quantum simulator of interacting erbium atoms in an optical lattice, we observe such a transition between one-dimensional crystalline symmetry-protected topological phases (CSPTs). We detect the critical point through non-local string order parameters and reveal its connection to the transition predicted between the Mott and Haldane insulators. Moreover, we demonstrate a striking property: stacking two identical systems eliminates the transition, confirming the predicted group structure and invertibility of SPTs. Finally, while introducing symmetry-breaking disorder also removes the transition, disorder averaging restores it. Consequently, the adjacent phases realize a form of mixed-state quantum order wherein the criticality between them depends on the observer's information. Our results demonstrate how topology and information influence quantum phase transitions, opening the doors to probing novel critical phenomena in programmable quantum matter.

Understanding emergent phenomena in out-of-equilibrium interacting many-body systems is an exciting frontier in physical science. While quantum simulators represent a promising approach to this long-standing problem, in practice it can be challenging to directly realize the required interactions, measure arbitrary observables, and mitigate errors. Here we use coherent mapping between the Rydberg and hyperfine qubits in a neutral atom array simulator to engineer and probe complex quantum dynamics. We combine efficient analog dynamics with fully programmable state preparation and measurement, leverage non-destructive readout for loss information and atomic qubit reuse, and use an atom reservoir for replacing lost atoms. With this analog-digital approach, we first demonstrate dynamical engineering of ring-exchange and particle hopping dynamics via Floquet driving and measure the spectral function of single excitations by evolving initial superposition states. Extending these techniques to a 271-site kagome lattice, we employ closed-loop optimization to target an out-of-equilibrium critical quantum spin liquid of the Rokhsar-Kivelson type. We observe the key features of such a state, including the absence of local order, many-body coherences between nearly equal-amplitude dimer configurations over up to 18 sites, and universal correlations consistent with predictions from field theory. Together, these results pave the way for using dynamical control in analog-digital quantum simulators to study complex quantum many-body systems.

In quantum many-body physics, one aims to understand emergent phenomena and effects of strong interactions, ideally by developing a simple theoretical picture. Recently, progress in quantum simulators has enabled the measurement of site resolved snapshots of Fermi-Hubbard systems at finite doping on square as well as triangular lattice geometries. These experimental advances pose the quest for theorists to analyze the ensuing data in order to gain insights into these prototypical, strongly correlated many-body systems. Here we employ machine learning techniques to optimize the mean-field parameters of a resonating valence bond (RVB) state through comparison with experimental data, thus determining a possible underlying simple model that is physically motivated and fully interpretable. We find that the resulting RVB states are capable of capturing two- as well as three-point correlations measured in experiments, even when they are not specifically used in the optimization. The analysis of the mean-field parameters and their doping dependence can be used to obtain physical insights and shed light on the nature of possible underlying quantum spin liquid states. Our results show that finite temperature data from Fermi-Hubbard quantum simulators can be well captured by RVB states. This work paves the way for a new, systematic analysis of data from numerical as well as quantum simulation of strongly correlated quantum many-body systems.

Fractional quantum Hall states are a cornerstone of topological physics, hosting fractionally charged quasiparticles with exotic statistics that promise to enable topologically protected quantum information processing. Among these, the Pfaffian state introduced by Moore and Read implements a p-wave pairing structure that supports excitations with non-Abelian exchange statistics. Despite extensive study in electronic systems, direct access to its pairing structure has remained limited. Here we realize a three-particle bosonic Pfaffian state of ultracold ^87Rb atoms in an optical lattice subject to a Floquet-engineered synthetic magnetic field. Using a Bayesian-optimized adiabatic protocol, we prepare a state exhibiting Pfaffian pairing correlations. Site-resolved measurements of multi-point density correlations reveal a pronounced suppression of short-range three-body coincidences, reflecting the underlying pairing structure. We further probe the state's transport response through Hall drift measurements. Our results establish a bottom-up approach to engineering non-Abelian topological order and lay the groundwork for future explorations of anyonic braiding in synthetic matter.

Creation and manipulation of entanglement with low error is essential in quantum information systems. In practice, two-qubit entangling gates constitute a dominant error source, limiting circuit depths and performance in fault-tolerant architectures. Using a neutral-atom quantum processor, we realize entangling CZ gates with a high Rabi frequency smooth-amplitude pulse, employing state-selective readout and qubit reuse for fast calibration, and achieve state-of-the-art fidelities of 99.854(4)

Strongly correlated materials feature multiple electronic orbitals, which are crucial to accurately understanding their many-body properties. In such multiband models, quantum interference can lead to flat energy bands with large degeneracy that gives rise to itinerant magnetic phases. We report on signatures of a ferrimagnetic state realized in a Lieb lattice with ultracold fermions, characterized by antialigned magnetic moments with antiferromagnetic correlations, and concomitant with a finite spin polarization. The signatures remain robust when increasing repulsive interactions from the weakly interacting to the Heisenberg regime and emerge when continuously tuning the lattice unit cell from a square to a Lieb geometry. Our flexible approach paves the way toward exploring exotic phases, such as quantum spin liquids in kagome lattices and heavy fermion behavior in Kondo models.

Quantum coherence underlies collective quantum phenomena and emerging quantum technologies. Quantum gas microscopes have transformed quantum simulation by providing projective snapshots of many-body states with single-atom resolution, but spatially resolved measurements of off-diagonal correlations have remained elusive. Here, using the Talbot effect, we introduce a quantum coherence microscope that maps off-diagonal correlations onto site-resolved density signals with near-single-site resolution. We use this technique to locally probe the superfluid-Mott transition in a layer of a three-dimensional optical lattice and to measure coherence beyond nearest neighbors in an engineered potential landscape. By mapping off-diagonal correlations onto density signals through controlled Talbot evolution, this work opens new possibilities for accessing observables beyond the density basis through tailored matter-wave evolution and recapture.

When quantum emitters couple indistinguishably to light, they can synchronize into a collective light matter system with radiative properties profoundly different from those of independent particles. To date, the resulting collective effects have largely been confined to point like or homogeneous ensembles. Here, we open access to a qualitatively new collective regime by realizing geometrically ordered, spatially extended atom arrays with subwavelength spacing. This establishes a fundamentally new platform in which collective emission is no longer confined to a single Dicke mode but instead emerges from an ordered network of photon mediated interactions. We find that 2D atom arrays undergo strong super and subradiant emission. Despite subwavelength spacing, we achieve site resolved imaging and directly observe the buildup of spatial correlations, demonstrating the transformation of cooperative decay into a strongly correlated many-body process. We observe extensive scaling of superradiance, uncover superradiant revivals, and reveal the ferromagnetic nature of superradiance and the antiferromagnetic nature of subradiance. Our results realize a novel programmable platform for exploring and utilizing dissipative many-body quantum physics, opening new possibilities for photon capture, storage, and atom photon entanglement.

Understanding doped Mott insulators is a fundamental goal in condensed matter physics, relevant to cuprate superconductors and other quantum materials1-3. The Hubbard model minimally describes such systems and has explicated some of their complex behaviour4,5. However, many open questions remain concerning the anomalous metallic states that emerge at low temperatures and intermediate doping and which, in cuprates, give rise to high-temperature superconductivity on cooling2,6,7. Here we observe a crossover between a normal metal and a pseudogapped metal in the Hubbard model using thermodynamic and spectroscopic measurements in a cold-atom quantum simulator, leveraging a recent several-fold reduction in achievable temperatures8. On cooling, the compressibility develops a maximum at intermediate doping, signalling an inflection point in the equation of state. We track this maximum versus interaction strength, revealing a line of thermodynamic anomalies in the phase diagram separating an underdoped from an overdoped metal at large interactions. Lattice modulation spectra in the underdoped regime show a loss of low-energy response, especially pronounced in the antinodal regions of the Brillouin zone, indicating a pseudogap. We use this signal to construct a pseudogap phase diagram versus interactions and doping. Our results experimentally establish and characterize the pseudogap metal in the Hubbard model, and suggest connections to charge order that can be studied in future work. Furthermore, this work demonstrates the utility of quantum simulation in addressing frontier problems in correlated electron physics.

Fractional statistics give rise to quantum behaviors that differ fundamentally from those of bosons and fermions. While two-dimensional anyons play a major role in strongly correlated systems and topological quantum computing, the nature of their one-dimensional (1D) counterparts remains the subject of intense debate, with renewed interest fueled by recent experimental progress. Theoretically, 1D anyons are predicted to host exotic many-body phases and quantum phase transitions, yet experimental signatures have remained elusive. Using ultracold atoms in an optical lattice, we prepare two-body ground states of the 1D anyon-Hubbard model by combining Hamiltonian engineering via quasiperiodic drives and adiabatic state manipulation. We uncover the effects of statistical interactions that lead to pseudo-fermionization and to the formation of chiral bound states when particles remain close together. Our results establish a link between lattice and continuum realizations of anyon models, and mark important steps towards the precise control of 1D anyons in both equilibrium and out-of-equilibrium settings.

Programmable arrays of ultracold atoms are a leading platform for quantum computation and simulation, enabling state-of-the-art implementations of quantum error correction, and analog simulations of Hubbard models that address open problems in condensed matter physics. In these systems, all local control is mediated through precisely shaped optical fields, and so the challenge of managing many-body quantum states becomes an exercise in optical design. In particular, one wishes for fast, flexible control with low disorder and heating, and access to large arrays with many atoms. An ideal optical system therefore must generate arbitrary patterns with high spatial resolution and low disorder, and alter these patterns on a timescale that is faster than the relevant atomic dynamics. Here, we present an optical system that is comparable to previous approaches in scale, while advancing all other axes. We demonstrate arbitrary pattern generation with 10^-3 intensity resolution, a frame rate of >84 MFPS (megaframes per second), and a spatial resolution of 83 × 52 beam waists (with 11 × 52 waists accessible via a single 40 GHz electro-optic modulator). These capabilities unlock a new class of experiments. We develop and numerically validate a scheme for fully programmable Hubbard models, with time-dependent control over local chemical potentials, tunneling amplitudes, on-site interactions, and patterns of artificial magnetic flux. The same architecture performs fast, arbitrary permutations of tweezers in 2D, decoupling optical constraints from the design of high-rate error-correcting codes.

Realizing universal fault-tolerant quantum computation is a key goal in quantum information science. By encoding quantum information into logical qubits utilizing quantum error correcting codes, physical errors can be detected and corrected, enabling substantial reduction in logical error rates. However, the set of logical operations that can be easily implemented on such encoded qubits is often constrained, necessitating the use of special resource states known as 'magic states' to implement universal, classically hard circuits. A key method to prepare high-fidelity magic states is to perform 'distillation', creating them from multiple lower fidelity inputs. Here we present the experimental realization of magic state distillation with logical qubits on a neutral-atom quantum computer. Our approach makes use of a dynamically reconfigurable architecture to encode and perform quantum operations on many logical qubits in parallel. We demonstrate the distillation of magic states encoded in d=3 and d=5 color codes, observing improvements of the logical fidelity of the output magic states compared to the input logical magic states. These experiments demonstrate a key building block of universal fault-tolerant quantum computation, and represent an important step towards large-scale logical quantum processors.

Quantum error correction (QEC) is believed to be essential for the realization of large-scale quantum computers. However, due to the complexity of operating on the encoded `logical' qubits, understanding the physical principles for building fault-tolerant quantum devices and combining them into efficient architectures is an outstanding scientific challenge. Here we utilize reconfigurable arrays of up to 448 neutral atoms to implement all key elements of a universal, fault-tolerant quantum processing architecture and experimentally explore their underlying working mechanisms. We first employ surface codes to study how repeated QEC suppresses errors, demonstrating 2.14(13)x below-threshold performance in a four-round characterization circuit by leveraging atom loss detection and machine learning decoding. We then investigate logical entanglement using transversal gates and lattice surgery, and extend it to universal logic through transversal teleportation with 3D [[15,1,3]] codes, enabling arbitrary-angle synthesis with logarithmic overhead. Finally, we develop mid-circuit qubit re-use, increasing experimental cycle rates by two orders of magnitude and enabling deep-circuit protocols with dozens of logical qubits and hundreds of logical teleportations with [[7,1,3]] and high-rate [[16,6,4]] codes while maintaining constant internal entropy. Our experiments reveal key principles for efficient architecture design, involving the interplay between quantum logic and entropy removal, judiciously using physical entanglement in logic gates and magic state generation, and leveraging teleportations for universality and physical qubit reset. These results establish foundations for scalable, universal error-corrected processing and its practical implementation with neutral atom systems.

A recent experiment has examined ultracold, fermionic, spin-1/2 6Li atoms in the Lieb lattice at different Hubbard repulsion U and filling fractions v [Lebrat et al., arXiv:2404.17555]. At v = 1/2 and small U, they observe an enhanced compressibility on the px,y sites, pointing to a flat band near the Fermi energy. At v =1/2 and large U they observe an insulating ferrimagnet. Both small and large U observations at v =1/2 are consistent with theoretical expectations. Surprisingly, near v = 1/4 and large U, they again observe a large px,y compressibility, pointing to a flat px,y band of fermions across the Fermi energy. Our Hartree-Fock computations near v = 1/4 find states with canted magnetism (and related spiral states) at large U, which possess nearly flat px,y bands near the Fermi level. We employ parton theories to describe quantum fluctuations of the magnetic order found in Hartree-Fock. We find a metallic state with Z2 fractionalization possessing gapless, fermionic, spinless "chargons" carrying Z2 gauge charges which have a nearly flat px,y band near their Fermi level: this fractionalized metal is also consistent with observations. Our DMRG study does not indicate the presence of magnetic order, and so is in agreement with a fractionalized ground state. Given the conventional ferrimagnetic insulator at v = 1/2, the Z2 fractionalized metal at v = 1/4 represents a remarkable realization of doping-induced fractionalization.

When particles move through a crystal or optical lattice, their motion can sometimes become frozen by strong external forces – yet collective motion may still emerge through subtle many-body effects. In this work, we explore such constrained dynamics by realizing a dipole-conserving Bose-Hubbard model, where single atoms are immobile but pairs of particles can move cooperatively while preserving the system's center of mass, i.e. the overall dipole moment of the particle distribution. Starting from a one-dimensional chain of ultracold bosonic atoms in an optical lattice, we generate localized dipole excitations consisting of a hole and a doublon using site-resolved optical potentials and characterize their quantum walks and scattering dynamics. Our study provides a bottom-up investigation of a Hamiltonian with kinetic constraints, and paves the way for exploring low-energy phases of fractonic matter in existing experimental platforms.

High-resolution fluorescence imaging of ultracold atoms and molecules is paramount to performing quantum simulation and computation in optical lattices and tweezers. Imaging durations in these experiments typically range from a millisecond to a second, significantly limiting the cycle time. In this work, we present fast, 2.4 microseconds single-atom imaging in lattices, with 99.4% fidelity - pushing the readout duration of neutral atom quantum platforms to be close to that of superconducting qubit platforms. Additionally, we thoroughly study the performance of accordion lattices. We also demonstrate number-resolved imaging without parity projection, which will facilitate experiments such as the exploration of high-filling phases in the extended Bose-Hubbard models, multi-band or SU(N) Fermi-Hubbard models, and quantum link models.

Ultracold fermionic atoms in optical lattices offer pristine realizations of Hubbard models1, which are fundamental to modern condensed-matter physics2,3. Despite notable advancements4-6, the accessible temperatures in these optical lattice material analogues are still too high to address many open problems7-10. Here we demonstrate a several-fold reduction in temperature6,11-13, bringing large-scale quantum simulations of the Hubbard model into an entirely new regime. This is accomplished by transforming a low-entropy product state into strongly correlated states of interest via dynamic control of the model parameters14,15, which is extremely challenging to simulate classically10. At half-filling, the long-range antiferromagnetic order is close to saturation, leading to a temperature of T / t = 0.0 5 - 0.05 + 0.06 based on comparisons with numerically exact simulations. Doped away from half-filling, it is exceedingly challenging to realize systematically accurate and predictive numerical simulations9. Importantly, we are able to use quantum simulation to identify a new pathway for achieving similarly low temperatures with doping. This is confirmed by comparing short-range spin correlations to state-of-the-art, but approximate, constrained-path auxiliary-field quantum Monte Carlo simulations16-18. Compared with the cuprates2,19,20, the reported temperatures correspond to a reduction from far above to below room temperature, at which physics such as the pseudogap and stripe phases may be expected3,19,21-24. Our work opens the door to quantum simulations that solve open questions in material science, develop synergies with numerical methods and theoretical studies, and lead to discoveries of new physics8,10.

Confinement of particles into bound states is a phenomenon spanning from high-energy to condensed matter physics, which can be studied in the framework of lattice gauge theories (LGTs). Achieving a comprehensive understanding of confinement continues to pose a major challenge, in particular at finite matter density and in the presence of strong quantum fluctuations. State-of-the-art quantum simulators constitute a promising platform to address this problem. Here we study confinement in coupled chains of ℤ_2 LGTs coupled to matter fields, that can be mapped to a mixed-dimensional (mixD) XXZ model. We perform large-scale numerical matrix-product state calculations to obtain the phase diagram of this model, in which we uncover striped phases formed by the ℤ_2 charges that can be melted at finite temperature or by increasing the tunneling rate. To explore this setting experimentally, we use our quantum simulator constituted by erbium atoms with dipolar interactions in a quantum gas microscope, and observe the predicted melting of a stripe phase by increasing the particle tunneling rate. Our explorative experimental studies of thermal deconfinement of ℤ_2 charges motivate our further theoretical study of the mixD ℤ_2 LGT, in which we predict a confined meson gas at finite temperature and low magnetization where thermal fluctuations destroy stripes but enable spontaneous commensurate spin order. Overall, we demonstrate that our platform can be used to study confinement in ℤ_2 LGTs coupled to matter fields, including long-range interactions and beyond one dimension, paving the way for future research of confinement in the quantum many-body regime.