Quantum phases of matter on a 256-atom programmable quantum simulator.

Motivated by far-reaching applications ranging from quantum simulations of complex processes in physics and chemistry to quantum information processing1, a broad effort is currently underway to build large-scale programmable quantum systems. Such systems provide insights into strongly correlated quantum matter2-6, while at the same time enabling new methods for computation7-10 and metrology11. Here we demonstrate a programmable quantum simulator based on deterministically prepared two-dimensional arrays of neutral atoms, featuring strong interactions controlled by coherent atomic excitation into Rydberg states12. Using this approach, we realize a quantum spin model with tunable interactions for system sizes ranging from 64 to 256 qubits. We benchmark the system by characterizing high-fidelity antiferromagnetically ordered states and demonstrating quantum critical dynamics consistent with an Ising quantum phase transition in (2 + 1) dimensions13. We then create and study several new quantum phases that arise from the interplay between interactions and coherent laser excitation14, experimentally map the phase diagram and investigate the role of quantum fluctuations. Offering a new lens into the study of complex quantum matter, these observations pave the way for investigations of exotic quantum phases, non-equilibrium entanglement dynamics and hardware-efficient realization of quantum algorithms.

Quantum Kibble-Zurek mechanism and critical dynamics on a programmable Rydberg simulator.

Quantum phase transitions (QPTs) involve transformations between different states of matter that are driven by quantum fluctuations1. These fluctuations play a dominant part in the quantum critical region surrounding the transition point, where the dynamics is governed by the universal properties associated with the QPT. Although time-dependent phenomena associated with classical, thermally driven phase transitions have been extensively studied in systems ranging from the early Universe to Bose-Einstein condensates2-5, understanding critical real-time dynamics in isolated, non-equilibrium quantum systems remains a challenge6. Here we use a Rydberg atom quantum simulator with programmable interactions to study the quantum critical dynamics associated with several distinct QPTs. By studying the growth of spatial correlations when crossing the QPT, we experimentally verify the quantum Kibble-Zurek mechanism (QKZM)7-9 for an Ising-type QPT, explore scaling universality and observe corrections beyond QKZM predictions. This approach is subsequently used to measure the critical exponents associated with chiral clock models10,11, providing new insights into exotic systems that were not previously understood and opening the door to precision studies of critical phenomena, simulations of lattice gauge theories12,13 and applications to quantum optimization14,15.

Quantum phases of Rydberg atoms on a kagome lattice.

We analyze the zero-temperature phases of an array of neutral atoms on the kagome lattice, interacting via laser excitation to atomic Rydberg states. Density-matrix renormalization group calculations reveal the presence of a wide variety of complex solid phases with broken lattice symmetries. In addition, we identify a regime with dense Rydberg excitations that has a large entanglement entropy and no local order parameter associated with lattice symmetries. From a mapping to the triangular lattice quantum dimer model, and theories of quantum phase transitions out of the proximate solid phases, we argue that this regime could contain one or more phases with topological order. Our results provide the foundation for theoretical and experimental explorations of crystalline and liquid states using programmable quantum simulators based on Rydberg atom arrays.

Emergent Z_{2} Gauge Theories and Topological Excitations in Rydberg Atom Arrays.

Strongly interacting arrays of Rydberg atoms provide versatile platforms for exploring exotic many-body phases and dynamics of correlated quantum systems. Motivated by recent experimental advances, we show that the combination of Rydberg interactions and appropriate lattice geometries naturally leads to emergent Z_{2} gauge theories endowed with matter fields. Based on this mapping, we describe how Rydberg platforms could realize two distinct classes of topological Z_{2} quantum spin liquids, which differ in their patterns of translational symmetry fractionalization. We also discuss the natures of the fractionalized excitations of these Z_{2} spin liquid states using both fermionic and bosonic parton theories and illustrate their rich interplay with proximate solid phases.

Emergent Glassy Behavior in a Kagome Rydberg Atom Array.

We present large-scale quantum Monte Carlo simulation results on a realistic Hamiltonian of kagome-lattice Rydberg atom arrays. Although the system has no intrinsic disorder, intriguingly, our analyses of static and dynamic properties on large system sizes reveal emergent glassy behavior in a region of parameter space located between two valence bond solid phases. The extent of this glassy region is demarcated using the Edwards-Anderson order parameter, and its phase transitions to the two proximate valence bond solids-as well as the crossover towards a trivial paramagnetic phase-are identified. We demonstrate the intrinsically slow (imaginary) time dynamics deep inside the glassy phase and discuss experimental considerations for detecting such a quantum disordered phase with numerous nearly degenerate local minima. Our proposal paves a new route to the study of real-time glassy phenomena and highlights the potential for quantum simulation of a distinct phase of quantum matter beyond solids and liquids in current-generation Rydberg platforms.

Complex Density Wave Orders and Quantum Phase Transitions in a Model of Square-Lattice Rydberg Atom Arrays.

We describe the zero-temperature phase diagram of a model of a two-dimensional square-lattice array of neutral atoms, excited into Rydberg states and interacting via strong van der Waals interactions. Using the density-matrix renormalization group algorithm, we map out the phase diagram and obtain a rich variety of phases featuring complex density wave orderings, upon varying lattice spacing and laser detuning. While some of these phases result from the classical optimization of the van der Waals energy, we also find intrinsically quantum-ordered phases stabilized by quantum fluctuations. These phases are surrounded by novel quantum phase transitions, which we analyze by finite-size scaling numerics and Landau theories. Our work highlights Rydberg quantum simulators in higher dimensions as promising platforms to realize exotic many-body phenomena.

Triangular lattice quantum dimer model with variable dimer density.

Quantum dimer models are known to host topological quantum spin liquid phases, and it has recently become possible to simulate such models with Rydberg atoms trapped in arrays of optical tweezers. Here, we present large-scale quantum Monte Carlo simulation results on an extension of the triangular lattice quantum dimer model with terms in the Hamiltonian annihilating and creating single dimers. We find distinct odd and even [Formula: see text] spin liquids, along with several phases with no topological order: a staggered crystal, a nematic phase, and a trivial symmetric phase with no obvious broken symmetry. We also present dynamic spectra of the phases, and note implications for experiments on Rydberg atoms.

Quantum field theory for the chiral clock transition in one spatial dimension.

We describe the quantum phase transition in the N -state chiral clock model in spatial dimension d = 1 . With couplings chosen to preserve time-reversal and spatial inversion symmetries, such a model is in the universality class of recent experimental studies of the ordering of pumped Rydberg states in a one-dimensional chain of trapped ultracold alkali atoms. For such couplings and N = 3 , the clock model is expected to have a direct phase transition from a gapped phase with a broken global ℤ N symmetry, to a gapped phase with the ℤ N symmetry restored. The transition has dynamical critical exponent z ≠ 1 , and so cannot be described by a relativistic quantum field theory. We use a lattice duality transformation to map the transition onto that of a Bose gas in d = 1 , involving the onset of a single boson condensate in the background of a higher-dimensional N -boson condensate. We present a renormalization group analysis of the strongly coupled field theory for the Bose gas transition in an expansion in 2 - d , with 4 - N chosen to be of order 2 - d . At two-loop order, we find a regime of parameters with a renormalization group fixed point which can describe a direct phase transition. We also present numerical density-matrix renormalization group studies of lattice chiral clock and Bose gas models for N = 3 , finding good evidence for a direct phase transition, and obtain estimates for z and the correlation length exponent ν .