A quantum processor based on coherent transport of entangled atom arrays

[coderman]

https://www.nature.com/articles/s41586-022-04592-6

A quantum processor based on coherent transport of entangled atom arrays

[Nature](https://www.nature.com/) volume 604, pages 451–456 (2022)[Cite this article](https://www.nature.com/articles/s41586-022-04592-6#citeas)

Abstract

The ability to engineer parallel, programmable operations between desired qubits within a quantum processor is key for building scalable quantum information systems[1](https://www.nature.com/articles/s41586-022-04592-6#ref-CR1),[2](https://www.nature.com/articles/s41586-022-04592-6#ref-CR2). In most state-of-the-art approaches, qubits interact locally, constrained by the connectivity associated with their fixed spatial layout. Here we demonstrate a quantum processor with dynamic, non-local connectivity, in which entangled qubits are coherently transported in a highly parallel manner across two spatial dimensions, between layers of single- and two-qubit operations. Our approach makes use of neutral atom arrays trapped and transported by optical tweezers; hyperfine states are used for robust quantum information storage, and excitation into Rydberg states is used for entanglement generation[3](https://www.nature.com/articles/s41586-022-04592-6#ref-CR3),[4](https://www.nature.com/articles/s41586-022-04592-6#ref-CR4),[5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5). We use this architecture to realize programmable generation of entangled graph states, such as cluster states and a seven-qubit Steane code state[6](https://www.nature.com/articles/s41586-022-04592-6#ref-CR6),[7](https://www.nature.com/articles/s41586-022-04592-6#ref-CR7). Furthermore, we shuttle entangled ancilla arrays to realize a surface code state with thirteen data and six ancillary qubits[8](https://www.nature.com/articles/s41586-022-04592-6#ref-CR8) and a toric code state on a torus with sixteen data and eight ancillary qubits[9](https://www.nature.com/articles/s41586-022-04592-6#ref-CR9). Finally, we use this architecture to realize a hybrid analogue–digital evolution[2](https://www.nature.com/articles/s41586-022-04592-6#ref-CR2) and use it for measuring entanglement entropy in quantum simulations[10](https://www.nature.com/articles/s41586-022-04592-6#ref-CR10),[11](https://www.nature.com/articles/s41586-022-04592-6#ref-CR11),[12](https://www.nature.com/articles/s41586-022-04592-6#ref-CR12), experimentally observing non-monotonic entanglement dynamics associated with quantum many-body scars[13](https://www.nature.com/articles/s41586-022-04592-6#ref-CR13),[14](https://www.nature.com/articles/s41586-022-04592-6#ref-CR14). Realizing a long-standing goal, these results provide a route towards scalable quantum processing and enable applications ranging from simulation to metrology.

Main

Quantum information systems derive their power from controllable interactions that generate quantum entanglement. However, the natural, local character of interactions limits the connectivity of quantum circuits and simulations. Non-local connectivity can be engineered via a global shared quantum data bus[15](https://www.nature.com/articles/s41586-022-04592-6#ref-CR15),[16](https://www.nature.com/articles/s41586-022-04592-6#ref-CR16),[17](https://www.nature.com/articles/s41586-022-04592-6#ref-CR17),[18](https://www.nature.com/articles/s41586-022-04592-6#ref-CR18), but in practice these approaches have been limited in either control or size. A number of visionary architectures to address this challenge have been proposed theoretically over the past two decades. On the basis of coherent, dynamical transport of quantum information using movable traps or photonic links, these techniques have been the subject of intensive experimental explorations across different platforms[1](https://www.nature.com/articles/s41586-022-04592-6#ref-CR1),[19](https://www.nature.com/articles/s41586-022-04592-6#ref-CR19),[20](https://www.nature.com/articles/s41586-022-04592-6#ref-CR20),[21](https://www.nature.com/articles/s41586-022-04592-6#ref-CR21),[22](https://www.nature.com/articles/s41586-022-04592-6#ref-CR22),[23](https://www.nature.com/articles/s41586-022-04592-6#ref-CR23),[24](https://www.nature.com/articles/s41586-022-04592-6#ref-CR24). However, progress has been limited to small-scale, few-qubit systems lacking either full connectivity, programmability or true parallelism.

Our approach to address this long-standing challenge utilizes dynamically reconfigurable arrays of entangled neutral atoms, shuttled by optical tweezers in two spatial dimensions (Fig. [1a](https://www.nature.com/articles/s41586-022-04592-6#Fig1)). Hyperfine states are used for storing and transporting quantum information between quantum operations, and excitation into Rydberg states is used for generating entanglement. Highly parallel operations are enabled via selective qubit operations in distinct zones that qubits are dynamically shuttled between. Taken together, these ingredients enable a powerful quantum information architecture, which we employ to realize applications including entangled-state generation, the creation of topological surface and toric code states, and hybrid analogue–digital quantum simulations.

Fig. 1: Quantum information architecture enabled by coherent transport of neutral atoms.
[figure 1](https://www.nature.com/articles/s41586-022-04592-6/figures/1)

a, In our approach, qubits are transported to perform entangling gates with distant qubits, enabling programmable and non-local connectivity. Atom shuttling is performed using optical tweezers, with high parallelism in two dimensions and between multiple zones allowing selective manipulations. Inset: the atomic levels used. The |0⟩, |1⟩ qubit states refer to the mF = 0 clock states of 87Rb, and |r⟩ is a Rydberg state used for generating entanglement between qubits (Extended Data Fig. [1b](https://www.nature.com/articles/s41586-022-04592-6#Fig5)). b, Atom images illustrating coherent transport of entangled qubits. Using a sequence of single-qubit and two-qubit gates, atom pairs are each prepared in the |Φ+⟩ Bell state (Methods), and are then separated by 110 μm over a span of 300 μs. c, Parity oscillations indicate that movement does not observably affect entanglement or coherence. For both the moving and the stationary measurements, qubit coherence is preserved using an XY8 dynamical decoupling sequence for 300 μs (Methods). d, Measured Bell-state fidelity as a function of separation speed over the 110 μm, showing that fidelity is unaffected for a move slower than 200 μs (average separation speed of 0.55 μm μs−1). Inset: normalizing by atom loss during the move results in constant fidelity, indicating that atom loss is the dominant error mechanism (see [Methods](https://www.nature.com/articles/s41586-022-04592-6#Sec7) for details).

[Full size image](https://www.nature.com/articles/s41586-022-04592-6/figures/1)

Entanglement transport in atom arrays

Our experiments utilize a two-dimensional (2D) atom array system described previously[25](https://www.nature.com/articles/s41586-022-04592-6#ref-CR25), with key upgrades to enable coherent transport and multiple layers of single-qubit and two-qubit gates. We store quantum information in magnetically insensitive clock states within the ground-state hyperfine manifold of 87Rb atoms[20](https://www.nature.com/articles/s41586-022-04592-6#ref-CR20), and implement robust single-qubit Raman rotations (scattering error per π pulse of about 7 × 10−5)[26](https://www.nature.com/articles/s41586-022-04592-6#ref-CR26), realized by composite pulses that are robust to pulse errors (Extended Data Fig. [3](https://www.nature.com/articles/s41586-022-04592-6#Fig7))[27](https://www.nature.com/articles/s41586-022-04592-6#ref-CR27). High-fidelity controlled-Z (CZ) entangling gates in the hyperfine basis {|0⟩, |1⟩} (Fig. [1a](https://www.nature.com/articles/s41586-022-04592-6#Fig1)) are implemented in parallel using global Rydberg excitation pulses on the |1⟩ ↔ |r⟩ Rydberg transition[5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5). For dynamic reconfiguration, we initialize atoms into two sets of traps: static traps generated by a spatial light modulator (SLM) and mobile traps generated by a crossed 2D acousto-optic deflector (AOD). To execute a specific circuit, we arrange qubits into desired pairs, perform Rydberg-mediated CZ gates on each pair simultaneously and then move all mobile traps in parallel to dynamically change the connectivity into the next desired qubit arrangement.

Figure [1](https://www.nature.com/articles/s41586-022-04592-6#Fig1) shows our ability to transport qubits across large distances while preserving entanglement and coherence[20](https://www.nature.com/articles/s41586-022-04592-6#ref-CR20). We initialize pairs at an atom–atom distance of 3 μm (Fig. [1b](https://www.nature.com/articles/s41586-022-04592-6#Fig1)) and then create a Bell state |Φ+⟩=12√(|00⟩+|11⟩)

in the hyperfine basis (Methods)[5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5). To measure the resulting entangled-state fidelity, we apply a variable single-qubit phase gate before a final π/2 pulse, resulting in oscillations of the two-atom parity ⟨σz1σz2⟩

(Fig. [1c](https://www.nature.com/articles/s41586-022-04592-6#Fig1))[5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5). We then repeat this experiment, but now move the atoms apart by 110 μm before applying the final π/2 pulse. Our transport protocol is optimized to suppress heating and loss by implementing cubic-interpolated atom trajectories (Methods), and is further accompanied by an eight-pulse XY8 robust dynamical decoupling sequence[28](https://www.nature.com/articles/s41586-022-04592-6#ref-CR28) to suppress dephasing. The resulting parity oscillations indicate that two-atom entanglement is unaffected by the transport process[20](https://www.nature.com/articles/s41586-022-04592-6#ref-CR20),[29](https://www.nature.com/articles/s41586-022-04592-6#ref-CR29). Performing this experiment as a function of movement speed[30](https://www.nature.com/articles/s41586-022-04592-6#ref-CR30) shows that the fidelity remains unchanged until the total separation speed becomes more than 0.55 μm μs−1, corresponding to the onset of atom loss (Fig. [1d](https://www.nature.com/articles/s41586-022-04592-6#Fig1)). We note that the entanglement transport in Fig. [1b](https://www.nature.com/articles/s41586-022-04592-6#Fig1) corresponds to moving quantum information across a region of space that can, in principle, host about 2,000 qubits (at an atom separation of 3 μm), on a timescale corresponding to <10−3 of the coherence time T2 (Extended Data Fig. [3](https://www.nature.com/articles/s41586-022-04592-6#Fig7)), directly enabling applications in large-scale quantum information systems.

Programmable circuits and graph states

To exemplify the ability to generate non-local connectivity between qubit arrays in parallel, we carry out the preparation of entangled graph states: a large class of useful quantum information states, with examples ranging from Greenberger–Horne–Zeilinger states and cluster states to quantum error correction (QEC) codes[31](https://www.nature.com/articles/s41586-022-04592-6#ref-CR31). Graph states are defined by initializing all qubits, located on the vertices of a geometric graph, in |+⟩≡|0⟩+|1⟩2√

and then performing CZ gates on the links between qubits (corresponding to the edges of the graph)[31](https://www.nature.com/articles/s41586-022-04592-6#ref-CR31). N-qubit graph states |G⟩ are associated with a set of N stabilizers, defined by Si=XiΠj∈uiZj

, where X and Z are the Pauli matrices, ui is the set of qubits (vertices) connected by an edge to qubit i, and Π denotes a product over qubit indices j (ref. [31](https://www.nature.com/articles/s41586-022-04592-6#ref-CR31)). The stabilizers each have +1 eigenvalue for the graph state |G⟩. Measuring these operators and their expectation values can be used to characterize the preparation of the target state.

As an example, Fig. [2a](https://www.nature.com/articles/s41586-022-04592-6#Fig2) shows the preparation of a one-dimensional (1D) cluster state, a graph state defined by a linear chain of qubits[32](https://www.nature.com/articles/s41586-022-04592-6#ref-CR32),[33](https://www.nature.com/articles/s41586-022-04592-6#ref-CR33). To realize this state, we perform one global, parallel layer of CZ gates on adjacent atom pairs, move half the atoms to form new pairs and then perform another parallel layer of CZ gates (Fig. [2a, b](https://www.nature.com/articles/s41586-022-04592-6#Fig2)). To probe the resultant 12-qubit cluster state, we measure the stabilizer set {Si} = {Zi−1XiZi+1} through readout in two measurement settings, given by a local π/2 rotation on either the odd or the even sublattice before projective measurement[34](https://www.nature.com/articles/s41586-022-04592-6#ref-CR34). We achieve the local rotation by moving one sublattice of qubits to a separate zone and then performing a rotation on the unmoved qubits with a homogeneous beam illuminating the experiment zone (Fig. [1a](https://www.nature.com/articles/s41586-022-04592-6#Fig1), Methods). We measure ⟨Si⟩ by analysing the resulting bit-string outputs and plot the resulting raw stabilizer measurements (Fig. [2c](https://www.nature.com/articles/s41586-022-04592-6#Fig2)). Across all 12 stabilizers, we find an average ⟨Si⟩ = 0.87(1) (Fig. [2c](https://www.nature.com/articles/s41586-022-04592-6#Fig2)) (accounting for state-preparation-and-measurement (SPAM) errors would yield ⟨Si⟩ = 0.91(1)), certifying biseparable entanglement in a cluster state (all ⟨Si⟩ > 0.5 (ref. [34](https://www.nature.com/articles/s41586-022-04592-6#ref-CR34))). The measured fidelities would correspond to a few-percent error per operation for a measurement-based quantum computation[32](https://www.nature.com/articles/s41586-022-04592-6#ref-CR32),[35](https://www.nature.com/articles/s41586-022-04592-6#ref-CR35).

Fig. 2: One- and two-dimensional graph states using dynamic entanglement transport.
[figure 2](https://www.nature.com/articles/s41586-022-04592-6/figures/2)

a, Generation of a 12-atom 1D cluster-state graph, created by initializing all qubits (vertices) in |+⟩ and applying CZ gates on the links (edges) between qubits. The atom images show the configuration for the first and second gate layers. b, Quantum circuit representation of the 1D cluster-state preparation and measurement. Dynamical decoupling is applied throughout all quantum circuits ([Methods](https://www.nature.com/articles/s41586-022-04592-6#Sec7)). c, Raw measured stabilizers of the resulting 1D cluster state, given by Si = Zi−1XiZi+1 (X1Z2 and Z11X12 for the edge qubits). d, Graph-state representation of the seven-qubit Steane code (colours represent stabilizer plaquettes). e, Circuit for preparing the Steane code logical |+⟩L state, performed in four parallel gate layers. f, Measured stabilizers and logical operators after preparing |+⟩L. Error detection is done by postselecting on measurements where all stabilizers are +1. For both the 1D cluster state and the Steane code, the stabilizers and logical operators are measured with two measurement settings (see text). Error bars represent 68% confidence intervals.

[Full size image](https://www.nature.com/articles/s41586-022-04592-6/figures/2)

An important class of graph states are QEC codes, where the graph-state stabilizers manifest as the stabilizers of the QEC code and can be measured to correct errors on an encoded logical qubit. As an example, we prepare the seven-qubit Steane code[6](https://www.nature.com/articles/s41586-022-04592-6#ref-CR6),[7](https://www.nature.com/articles/s41586-022-04592-6#ref-CR7), a topological colour code depicted by the graph in Fig. [2d](https://www.nature.com/articles/s41586-022-04592-6#Fig2), in the logical state |+⟩L. To prepare this state, we initialize all qubits in |+⟩, apply CZs on the links between qubits (in four parallel layers; Extended Data Fig. [5](https://www.nature.com/articles/s41586-022-04592-6#Fig9)) and then rotate either of the two sublattices for measuring stabilizers (Fig. [2e](https://www.nature.com/articles/s41586-022-04592-6#Fig2)). After sublattice rotation, six of the graph-state stabilizers transform into the six Steane code stabilizers, given by four-body products of Xi or Zi. Figure [2f](https://www.nature.com/articles/s41586-022-04592-6#Fig2) shows the raw measured expectation values of these six stabilizers. The seventh graph-state stabilizer transforms into the logical qubit operator XL and has eigenvalue +1 for the graph state |G⟩, while anticommuting with logical ZL. Accordingly, in Fig. [2f](https://www.nature.com/articles/s41586-022-04592-6#Fig2), we find ⟨XL⟩ = 0.71(2) and ⟨ZL⟩ = −0.02(3), demonstrating the preparation of the logical qubit state |+⟩L. Moreover, we perform error detection by post-selecting on measurement outcomes where all measured stabilizers yield +1 (refs. [36](https://www.nature.com/articles/s41586-022-04592-6#ref-CR36),[37](https://www.nature.com/articles/s41586-022-04592-6#ref-CR37); with 66(1)% probability of no detected errors). Using this procedure, we obtain corrected values of ⟨X¯L⟩=0.991+0.004−0.007

and, ⟨Z¯L⟩=−0.03(3)

demonstrating the error-detecting properties of the Steane code graph (see Extended Data Fig. [7](https://www.nature.com/articles/s41586-022-04592-6#Fig11) for error correction and logical operations).

Topological states with ancilla arrays

We next make use of transportable ancillary qubit arrays to mediate quantum operations between remote qubits[1](https://www.nature.com/articles/s41586-022-04592-6#ref-CR1). Owing to the ability to quickly move arrays of atoms across the entire system, the use of ancillary qubits naturally complements our movement capabilities. Specifically, we employ ancillas for state preparation by mediating entanglement between physical qubits that never directly interact, followed by projective measurement of the ancilla array (performed simultaneously with the measurement of the data qubits), a form of measurement-based quantum computation[32](https://www.nature.com/articles/s41586-022-04592-6#ref-CR32),[35](https://www.nature.com/articles/s41586-022-04592-6#ref-CR35). In particular, we prepare topological surface code and toric code states[8](https://www.nature.com/articles/s41586-022-04592-6#ref-CR8),[9](https://www.nature.com/articles/s41586-022-04592-6#ref-CR9),[38](https://www.nature.com/articles/s41586-022-04592-6#ref-CR38),[39](https://www.nature.com/articles/s41586-022-04592-6#ref-CR39), whose states are more difficult to construct by direct CZ gates between physical qubits (requiring an extensive number of layers[8](https://www.nature.com/articles/s41586-022-04592-6#ref-CR8),[40](https://www.nature.com/articles/s41586-022-04592-6#ref-CR40)). For these codes, the measured values of the ancilla qubits simply redefine the stabilizers and are handled in-software for practical QEC operation[38](https://www.nature.com/articles/s41586-022-04592-6#ref-CR38). As the redefinition is applied in-software, without physical intervention, the projective measurements on the ancillae commute with all operations on the data qubits and can be done at any time, and so we measure all qubits simultaneously.

Figure [3a](https://www.nature.com/articles/s41586-022-04592-6#Fig3) shows the preparation of a 19-qubit graph state creating the |+⟩L logical state of the surface code[8](https://www.nature.com/articles/s41586-022-04592-6#ref-CR8),[38](https://www.nature.com/articles/s41586-022-04592-6#ref-CR38). The surface code is defined by X-plaquette and Z-star stabilizers, and logical operators XL (ZL) are defined as strings of X (Z) products across the height (width) of the graph. To prepare this state, ancillas are moved to perform CZ gates with each of their four neighbours and are then measured, projecting the data qubits into the surface code state. The graph-state stabilizers now transform into the X plaquettes, the Z stars (with value ±1 for a measurement outcome of ±1 of the central ancilla) and the logical XL operator[35](https://www.nature.com/articles/s41586-022-04592-6#ref-CR35),[41](https://www.nature.com/articles/s41586-022-04592-6#ref-CR41). Remarkably, this procedure creates a topologically ordered state in a constant-depth circuit[35](https://www.nature.com/articles/s41586-022-04592-6#ref-CR35),[40](https://www.nature.com/articles/s41586-022-04592-6#ref-CR40), where measured ancilla values can be used for redefining stabilizers, which can be handled in-software for practical QEC operation[38](https://www.nature.com/articles/s41586-022-04592-6#ref-CR38). Figure [3b](https://www.nature.com/articles/s41586-022-04592-6#Fig3) shows the measured expectation values of the 12 resulting stabilizers, as well as the logical operator expectation values with and without error detection. We find a raw value of ⟨XL⟩ = 0.64(3), and a corrected value of ⟨X¯L⟩=1+0−0.01

using the measured stabilizers for error detection (with 35(1)% probability of no detected errors), demonstrating the preparation of this topological QEC state (see also Extended Data Fig. [7](https://www.nature.com/articles/s41586-022-04592-6#Fig11), showing the expected attributes for all prepared error-protected logical states).

Fig. 3: Topological surface code and toric code states using mobile ancilla qubit arrays.
[figure 3](https://www.nature.com/articles/s41586-022-04592-6/figures/3)

a, Graph state realizing the surface code. Left: the circuit depicts formation of the graph state by use of mobile ancilla qubits; each move corresponds to performing a CZ gate with a neighbouring data qubit (illustrated in box). The logical |+⟩L state is created upon projective measurement of the ancilla qubits in the X basis. Right: stabilizers and logical operators of the code. b, Measured X-plaquette and Z-star stabilizers of the resultant surface code, along with logical operators with and without error detection (implemented in postselection). c, Implementation of the toric code. Top: graph state realizing the two logical-qubit product state |+⟩(1)L|+⟩(2)L

of the toric code upon projective measurement of the ancilla qubits in the X basis. Bottom: images showing the movement steps implemented in creating and measuring the toric code state ([Supplementary Video 1](https://www.nature.com/articles/s41586-022-04592-6#MOESM2)). The blue shading in the final image represents a local rotation on the data qubit zone. d, Measured X-plaquette and Z-star stabilizers, along with logical operators for the two logical qubits with and without error detection (implemented in postselection).

[Full size image](https://www.nature.com/articles/s41586-022-04592-6/figures/3)

Although surface code states have previously been prepared with other methods, our transport capabilities allow us to use the full range of motion of ancilla qubits across the entire qubit array to enable periodic boundary conditions and realize the toric code state on a torus[9](https://www.nature.com/articles/s41586-022-04592-6#ref-CR9). To this end, we create the 24-qubit graph state shown in Fig. [3c](https://www.nature.com/articles/s41586-022-04592-6#Fig3) by performing five layers of parallel gates and moving the ancillae to their separate zone for readout in a separate basis (see also [Supplementary Video 1](https://www.nature.com/articles/s41586-022-04592-6#MOESM2) showing the full atom trajectory). The state we prepare has seven (owing to periodic boundary conditions) independent X plaquettes and seven independent Z stars. Moreover, owing to the topological properties of this graph, two independent logical qubits can be encoded with logical operators X(1)L,Z(1)L

and X(2)L,Z(2)L that wrap around the entire torus along two topologically distinct directions[9](https://www.nature.com/articles/s41586-022-04592-6#ref-CR9). Upon projective measurement of the ancilla qubits in the X basis we create the toric code state |+⟩(1)L|+⟩(2)L. State preparation is verified in Fig. [3d](https://www.nature.com/articles/s41586-022-04592-6#Fig3) by measuring the toric code stabilizers. For the two encoded logical qubits, we find raw logical qubit expectation values of ⟨X(1)L⟩=0.64(2) and ⟨X(2)L⟩=0.38(2), and error-detected values of ⟨X¯(1)L⟩=1+0−0.01 and ⟨X¯(2)L⟩=0.92+0.02−0.03 (with 20(1)% probability of no detected errors), demonstrating the preparation of the toric code. We note that the different expectation values of the corrected logical qubits originate from the aspect ratio of our torus, where X(1)L and X(2)L

are protected to code distance d = 4 and d = 2, respectively (see also Extended Data Fig. [7](https://www.nature.com/articles/s41586-022-04592-6#Fig11)). Our measured fidelities are in good agreement with numerical simulations of the circuit (Extended Data Fig. [6](https://www.nature.com/articles/s41586-022-04592-6#Fig10)), wherein each qubit experiences a per-layer error rate independent of the number of qubits or the shuttling process, indicating that errors in CZ gates (fidelity of about 97.5%; Methods[5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5)) constitute our dominant error source.

Hybrid analogue–digital circuits

Having established atom movement for realizing digital circuits, we now explore the applications to quantum simulation. In particular, we perform hybrid, modular quantum circuits composed of analogue Hamiltonian evolution, reconfiguration and digital gates (Fig. [4a](https://www.nature.com/articles/s41586-022-04592-6#Fig4)). Together, these tools open a wide variety of possibilities in quantum simulation and many-body physics. As a specific example, we measure the Renyi entanglement entropy after a quantum quench by effectively interfering two copies of a many-body system[10](https://www.nature.com/articles/s41586-022-04592-6#ref-CR10),[11](https://www.nature.com/articles/s41586-022-04592-6#ref-CR11).

Fig. 4: Dynamic reconfigurability for hybrid analogue–digital quantum simulation.
[figure 4](https://www.nature.com/articles/s41586-022-04592-6/figures/4)

a, Hybrid quantum circuit combining coherent atom transport with analogue Hamiltonian evolution and digital quantum gates. b, Measuring entanglement entropy in a many-body Rydberg system via two-copy interferometry. c, Measured half-chain Renyi entanglement entropy after many-body dynamics following quenches on two eight-atom systems. Quenching from |gggg...⟩ (|g⟩ ≡ |1⟩) results in rapid entropy growth and saturation, signifying quantum thermalization. Quenching from |rgrg...⟩ reveals a significantly slower growth of entanglement entropy. d, Measuring the mutual information at 0.5-μs quench time reveals a volume-law scaling for the thermalizing |gggg...⟩ state, and an area-law scaling for the scarring |rgrg...⟩ state. e, The single-site Renyi entropies for sites in the middle of the chain quickly increase and saturate for the |gggg...⟩ quench, but show large oscillations for the |rgrg...⟩ quench. The solid curves are results of exact numerical simulations for the isolated quantum system under HRyd with no free parameters (see [Methods](https://www.nature.com/articles/s41586-022-04592-6#Sec7) for details of data processing). Error bars represent one standard deviation.

[Full size image](https://www.nature.com/articles/s41586-022-04592-6/figures/4)

Figure [4b](https://www.nature.com/articles/s41586-022-04592-6#Fig4) illustrates the experimental procedure. After initializing both copies with all qubits in |1⟩, we independently evolve each copy under the Rydberg Hamiltonian HRyd for a time t, generating an entangled many-body state in the {|1⟩, |r⟩} basis (Methods)[13](https://www.nature.com/articles/s41586-022-04592-6#ref-CR13). Raman and Rydberg π pulses then map |1⟩ → |0⟩ and |r⟩ → |1⟩, transferring the entangled many-body state into the long-lived and non-interacting {|0⟩, |1⟩} basis[42](https://www.nature.com/articles/s41586-022-04592-6#ref-CR42). Finally, we measure entanglement entropy by rearranging the system and interfering each qubit in the first copy with its identical twin in the second copy, by use of a Bell measurement circuit. Measuring twins in the Bell basis detects occurences of the antisymmetric singlet state |Ψ−⟩=|01⟩−|10⟩2√

, the presence of which indicates that subsystems of the two copies were in different states owing to entanglement with the rest of the many-body system[10](https://www.nature.com/articles/s41586-022-04592-6#ref-CR10),[11](https://www.nature.com/articles/s41586-022-04592-6#ref-CR11). Quantitatively, analysing the number parity of observed singlets within subsystem A yields the purity Tr[ρ2A] of the reduced density matrix ρA, and thus yields the second-order Renyi entanglement entropy S2(A)=−log2Tr[ρ2A]

(Methods). This measurement circuit provides the Renyi entropy of any constituent subsystem of our whole closed quantum system, where the calculation over any desired subsystem A is simply performed in data processing[10](https://www.nature.com/articles/s41586-022-04592-6#ref-CR10),[11](https://www.nature.com/articles/s41586-022-04592-6#ref-CR11).

We use this method to probe the growth of entanglement entropy produced by many-body dynamics (see [Methods](https://www.nature.com/articles/s41586-022-04592-6#Sec7) for additional benchmarking of the technique). Specifically, we study the evolution of two eight-atom copies under the Rydberg Hamiltonian, subject to the nearest-neighbour blockade constraint[4](https://www.nature.com/articles/s41586-022-04592-6#ref-CR4),[13](https://www.nature.com/articles/s41586-022-04592-6#ref-CR13). Upon a rapid quench from an initial state with all atoms in the ground state |g⟩ ≡ |1⟩, we observe that the half-chain Renyi entanglement entropy quickly grows and saturates (Fig. [4c](https://www.nature.com/articles/s41586-022-04592-6#Fig4)), a process corresponding to quantum thermalization[12](https://www.nature.com/articles/s41586-022-04592-6#ref-CR12). By analysing the Renyi mutual information IAB = S2(A) +  S2(B) − S2(AB) between the leftmost n atoms in the chain (A) and the complement subsystem of the rightmost 8 − n atoms (B), we find a volume-law scaling in the resulting state (Fig. [4d](https://www.nature.com/articles/s41586-022-04592-6#Fig4))[11](https://www.nature.com/articles/s41586-022-04592-6#ref-CR11),[12](https://www.nature.com/articles/s41586-022-04592-6#ref-CR12).

Although such thermalizing dynamics is generically expected in strongly interacting many-body systems, remarkably, it was demonstrated previously that for certain initial states this system can evade thermalization. Underpinned by special, non-thermal eigenstates called quantum many-body scars[13](https://www.nature.com/articles/s41586-022-04592-6#ref-CR13),[14](https://www.nature.com/articles/s41586-022-04592-6#ref-CR14),[43](https://www.nature.com/articles/s41586-022-04592-6#ref-CR43), these states were theoretically predicted to feature dynamics associated with a slow, non-monotonic entanglement growth. Figure [4](https://www.nature.com/articles/s41586-022-04592-6#Fig4) shows the measurement of entanglement properties of many-body scars following a rapid quench from the initial state |Z2⟩≡|rgrg...⟩

, initialized by applying local Rydberg π pulses (Methods). We find that the rate of entropy growth for this initial state is significantly suppressed, and the mutual information reveals an area-law scaling (Fig. [4d](https://www.nature.com/articles/s41586-022-04592-6#Fig4)). Furthermore, Fig. [4e](https://www.nature.com/articles/s41586-022-04592-6#Fig4) shows the single-site entropy in the middle of the chain, demonstrating rapid growth and saturation for the thermalizing |gggg...⟩ state but large oscillations for the |Z2⟩ state[13](https://www.nature.com/articles/s41586-022-04592-6#ref-CR13),[14](https://www.nature.com/articles/s41586-022-04592-6#ref-CR14). Remarkably, the data show that when sites of one sublattice return to low entropy, the other sublattice goes to high entropy; this reveals that the scar dynamics entangle distant atoms (of the same sublattice) while disentangling nearest neighbours, even with only nearest-neighbour interactions ([Methods](https://www.nature.com/articles/s41586-022-04592-6#Sec7)). These measurements reveal non-trivial aspects of quantum many-body scars, and constitute the direct observation of exotic entanglement phenomena in a many-body system.

These observations are in excellent agreement with exact numerical simulations in the isolated system (lines plotted in Fig. [4c, e](https://www.nature.com/articles/s41586-022-04592-6#Fig4), Extended Data Fig. [10](https://www.nature.com/articles/s41586-022-04592-6#Fig14)). Moreover, whereas the single-site purity approaches that of a fully mixed state, our global purity (a 16-body observable composed of three-level systems) remains more than 100 times that of a fully mixed state (Extended Data Fig. [9](https://www.nature.com/articles/s41586-022-04592-6#Fig13)), altogether demonstrating the high accuracy and fidelity of our circuit-based technique. These results demonstrate that combining atom movement, many-body Hamiltonian evolution and digital quantum circuits yields powerful tools for simulating and probing the quantum physics of complex systems.

Discussion and outlook

Our experiments demonstrate highly parallel coherent qubit transport and entanglement enabling a powerful quantum information architecture. The present techniques can be extended along a number of directions. Local Rydberg excitation on subsets of qubit pairs would eliminate residual interactions from unintended atoms, allowing parallel, independent operations on arrays with significantly higher qubit densities. Two-qubit gate fidelity can be improved using higher Rydberg laser power or more efficient delivery methods, as well as more advanced atom cooling[44](https://www.nature.com/articles/s41586-022-04592-6#ref-CR44). These technical improvements should allow for scaling to deep quantum circuits operating on thousands of neutral atom qubits. These upgrades can be additionally supplemented by more sophisticated local single-qubit control employing, for example, parallel Raman excitation through acousto-optic modulator arrays[16](https://www.nature.com/articles/s41586-022-04592-6#ref-CR16). Mid-circuit readout can be implemented by moving ancillas into a separate zone and imaging using, for example, avalanche photodiode arrays within a few hundred microseconds[45](https://www.nature.com/articles/s41586-022-04592-6#ref-CR45).

Our method has a clear potential for realizing scalable QEC[46](https://www.nature.com/articles/s41586-022-04592-6#ref-CR46). For example, the procedure demonstrated in Fig. [3c](https://www.nature.com/articles/s41586-022-04592-6#Fig3) can be used for syndrome extraction in a practical QEC sequence, wherein ancillas are entangled with their data qubit neighbours and then moved to a separate zone for mid-circuit readout. We estimate that an entire QEC round can be implemented within a millisecond, much faster than the measured T2 > 1 s, and with projected fidelity improvements theoretically surpassing the surface code threshold (Methods). We emphasize that such a mid-circuit readout is essential for realizing scalable fault-tolerant quantum computation. Furthermore, the ability to reconfigure and interlace our arrays will allow efficient, parallel execution of transversal entangling gates between many logical qubits[38](https://www.nature.com/articles/s41586-022-04592-6#ref-CR38),[47](https://www.nature.com/articles/s41586-022-04592-6#ref-CR47). In addition, these techniques also enable implementation of higher-dimensional or non-local error-correcting codes with more favourable properties[48](https://www.nature.com/articles/s41586-022-04592-6#ref-CR48),[49](https://www.nature.com/articles/s41586-022-04592-6#ref-CR49). Together, these ingredients could enable an approach to universal, fault-tolerant quantum computing with thousands of physical qubits.

Our dynamically reconfigurable architecture also opens many opportunities for digital and analogue quantum simulations. For example, our hybrid approach can be extended to probing the entire entanglement spectrum[50](https://www.nature.com/articles/s41586-022-04592-6#ref-CR50), simulating wormhole creation[51](https://www.nature.com/articles/s41586-022-04592-6#ref-CR51), performing many-body purification[52](https://www.nature.com/articles/s41586-022-04592-6#ref-CR52) and engineering non-equilibrium states[53](https://www.nature.com/articles/s41586-022-04592-6#ref-CR53). Entanglement transport could also empower metrological applications such as creating distributed states for probing gravitational gradients[54](https://www.nature.com/articles/s41586-022-04592-6#ref-CR54). Finally, our approach can facilitate quantum networking between separated arrays, paving the way towards large-scale quantum information systems[29](https://www.nature.com/articles/s41586-022-04592-6#ref-CR29),[55](https://www.nature.com/articles/s41586-022-04592-6#ref-CR55) and distributed quantum metrology[54](https://www.nature.com/articles/s41586-022-04592-6#ref-CR54),[56](https://www.nature.com/articles/s41586-022-04592-6#ref-CR56).

Methods

Dynamic reconfiguration in 2D tweezer arrays

Our experiments utilize the same apparatus described previously in ref. [25](https://www.nature.com/articles/s41586-022-04592-6#ref-CR25). Inside our vacuum cell, 87Rb atoms are loaded from a magneto-optical trap into a backbone array of programmable optical tweezers generated by an SLM[57](https://www.nature.com/articles/s41586-022-04592-6#ref-CR57). Atoms are rearranged in parallel into defect-free target positions in this SLM backbone[57](https://www.nature.com/articles/s41586-022-04592-6#ref-CR57) by additional optical tweezers generated from a crossed 2D AOD. Following the rearrangement procedure, we transfer selected atoms from the static SLM traps back into the mobile AOD traps, and then move these mobile atoms to their starting positions in the quantum circuit. During this entire process, the atoms are cooled with polarization-gradient cooling. Before running the quantum circuit, we take a camera image of the atoms in their initial starting positions, and following the circuit we take a final camera image to detect qubit states |0⟩ (atom presence) and |1⟩ (atom loss, following resonant pushout). We postselect all data on finding perfect rearrangement of the AOD and SLM atoms before running the circuit. In all experiments here, each atom remains in a single static or single mobile trap throughout the duration of the quantum circuit[20](https://www.nature.com/articles/s41586-022-04592-6#ref-CR20),[58](https://www.nature.com/articles/s41586-022-04592-6#ref-CR58),[59](https://www.nature.com/articles/s41586-022-04592-6#ref-CR59).

The crossed AOD system is composed of two independently controlled AODs (AA Opto Electronic DTSX-400) for x and y control of the beam positions. Both AODs are driven by independent arbitrary waveforms, which are generated by a dual-channel arbitrary waveform generator (M4i.6631-x8 by Spectrum Instrumentation) and then amplified through independent MW amplifiers (Minicircuits ZHL-5W-1). The time-domain arbitrary waveforms are composed of multiple frequency tones corresponding to the x and y positions of columns and rows, which are independently changed as a function of time for steering around the AOD-trapped atoms dynamically; the full x and y waveforms are calculated by adding together the time-domain profile of all frequency components with a given amplitude and phase for each component. For running quantum circuits, we program the positions of the AOD atoms at each gate location and then smoothly interpolate (with a cubic profile) the AOD frequencies as a function of time between gate positions. The cubic profile enacts a constant jerk onto the atoms, which allows us to move roughly five-to-ten-times faster (without heating and loss) than if we move at a constant velocity (linear profile). In our movement protocol, we only do stretches, compressions and translations of the AOD trap array: that is, the AOD rows and columns never cross each other to avoid atom loss and heating associated with two frequency components crossing each other.

We homogenize the AOD tweezer intensity throughout the whole atom trajectory to minimize dephasing induced by a time-varying magnitude of differential light shifts. To this end, we use a reference camera in the image plane to gauge the intensity of each AOD tweezer at each gate location and homogenize by varying the amplitude of each frequency component; during motion between two locations, we interpolate the amplitude of each individual frequency component.

The SLM tweezer light (830 nm) and the AOD tweezer light (828 nm) are generated by two separate, free-running titanium:sapphire lasers (M Squared, 18-W pump). Projected through a 0.5 numerical aperture objective, the SLM tweezers have a waist of roughly 900 nm (roughly 1,000 nm for AODs). When loading the atoms, the trap depths are about 2π × 16 MHz, with radial trap frequencies of about 2π × 80 kHz, and when running quantum circuits the trap depths are about 2π × 4 MHz, with radial trap frequencies of about 2π × 40 kHz.

Raman laser system

Fast, high-fidelity single-qubit manipulations are critical ingredients of the quantum circuits demonstrated in this work. To this end, we use a high-power 795-nm Raman laser system for driving global single-qubit rotations between magnetic sublevel mF = 0 clock states. This Raman laser system is based on dispersive optics, developed and described in ref. [26](https://www.nature.com/articles/s41586-022-04592-6#ref-CR26). The 795-nm light (Toptica TA pro, 1.8 W) is phase-modulated by an electro-optic modulator (Qubig), which is driven by microwaves at 3.4 GHz (Stanford Research Systems SRS SG384) that are doubled to 6.8 GHz and amplified. The laser phase modulation is converted to amplitude modulation for driving Raman transitions through use of a chirped Bragg grating (Optigrate)[26](https://www.nature.com/articles/s41586-022-04592-6#ref-CR26). IQ (in-phase and quadrature) control of the SG384 is used for frequency and phase control of the microwaves, which are imprinted onto the laser amplitude modulation and thus give us direct frequency and phase control over the hyperfine qubit drive.

The Raman laser illuminates the atom plane from the side in a circularly polarized elliptical beam with waists of 40 μm and 560 μm on the thin axis and the tall axis, respectively, with a total average optical power of 150 mW on the atoms. The large vertical extent ensures <1% inhomogeneity across the atoms, and shot-to-shot fluctuations in the laser intensity are also <1%. For Figs. [1](https://www.nature.com/articles/s41586-022-04592-6#Fig1)–[3](https://www.nature.com/articles/s41586-022-04592-6#Fig3), we operate our Raman laser at a blue-detuned intermediate-state detuning of 180 GHz, resulting in two-photon Rabi frequencies of 1 MHz and an estimated scattering error per π pulse of 7 × 10−5 (that is, 1 scattering event per 15,000 π pulses)[26](https://www.nature.com/articles/s41586-022-04592-6#ref-CR26). For Fig. [4](https://www.nature.com/articles/s41586-022-04592-6#Fig4), to shorten the duration of the coherent mapping pulse sequence, we increase the Raman laser power and operate at a smaller blue-detuned intermediate-state detuning of 63 GHz, with a corresponding two-photon Rabi frequency of 3.2 MHz and an estimated scattering error per π pulse of 2 × 10−4.

Robust single-qubit rotations

For almost all single-qubit rotations in this work (other than XY8 and XY16 self-correcting sequences), we implement robust single-qubit rotations in the form of composite pulse sequences. These composite pulse sequences are well known in the NMR community[27](https://www.nature.com/articles/s41586-022-04592-6#ref-CR27),[60](https://www.nature.com/articles/s41586-022-04592-6#ref-CR60) and can be highly insensitive to pulse errors such as amplitude or detuning miscalibrations. Our dominant source of coherent single-qubit errors arise from ≲1% amplitude drifts and inhomogeneity across the array; as such, we primarily use the ‘BB1’ (broadband 1) pulse sequence, which is a sequence of four pulses that implements an arbitrary rotation on the Bloch sphere while being insensitive to amplitude errors to sixth order[27](https://www.nature.com/articles/s41586-022-04592-6#ref-CR27),[60](https://www.nature.com/articles/s41586-022-04592-6#ref-CR60). We benchmark the performance of these robust pulses in Extended Data Fig. [3a](https://www.nature.com/articles/s41586-022-04592-6#Fig7). Furthermore, by applying a train of BB1 pulses, we find an accumulated error consistent with the estimated scattering limit (not plotted here), suggesting that the scattering limit roughly represents our single-qubit rotation infidelities (about 3 × 10−4 error per BB1 pulse owing to the increased length of the composite pulse sequence). Randomized benchmarking[61](https://www.nature.com/articles/s41586-022-04592-6#ref-CR61) can be applied in future studies to further study single-qubit rotation fidelity.

Qubit coherence and dynamical decoupling

In our 830-nm traps, hyperfine qubit coherence is characterized by inhomogeneous dephasing time T∗2

 = 4 ms (not plotted here), T2 = 1.5 s (XY16 with 128 total π pulses) and relaxation time T1 = 4 s (including atom loss) (Extended Data Fig. [3b, c](https://www.nature.com/articles/s41586-022-04592-6#Fig7)). All of our experiments in this work are performed in a d.c. magnetic field of 8.5 G. Coherence can be further improved by using further-detuned optical tweezers (with trap depth held constant, the tweezer differential lightshifts decrease as 1/Δ and 1/T1 decreases as 1/Δ3 (ref. [62](https://www.nature.com/articles/s41586-022-04592-6#ref-CR62)), where Δ is the detuning of the trap wavelength) and shielding against magnetic field fluctuations. For practical QEC operation, atom loss can be detected in a hardware-efficient manner[46](https://www.nature.com/articles/s41586-022-04592-6#ref-CR46) and the atom then replaced from a reservoir, which could in principle be continuously reloaded by a magneto-optical trap for reaching arbitrarily deep circuits.

All of our transport sequences[20](https://www.nature.com/articles/s41586-022-04592-6#ref-CR20),[58](https://www.nature.com/articles/s41586-022-04592-6#ref-CR58),[59](https://www.nature.com/articles/s41586-022-04592-6#ref-CR59) are accompanied by dynamical decoupling sequences. The number of pulses we use is a trade-off between preserving qubit coherence while minimizing pulse errors. We interchange between two types of dynamical decoupling sequence: XY8 and XY16 sequences, composed of phase-alternated individual π pulses that are self-correcting for amplitude and detuning errors[28](https://www.nature.com/articles/s41586-022-04592-6#ref-CR28),[63](https://www.nature.com/articles/s41586-022-04592-6#ref-CR63), and Carr–Purcell–Meiboom–Gill (CPMG)-type dynamical decoupling sequences composed of robust BB1 pulses. The CPMG-BB1 sequence is more robust to amplitude errors but incurs more scattering error. We empirically optimize for any given experiment by choosing between these different sequences and a variable number of decoupling π pulses, optimizing on either single-qubit coherence (including the movement) or the final signal. Typically, our decoupling sequences are composed of a total 12–18 π pulses.

Movement effects on atom heating and loss

We study here the effects of movement on atom loss and heating in the harmonic oscillator potential given by the tweezer trap. Motion of the trap potential is equivalent to the non-inertial frame of reference where the harmonic oscillator potential is stationary but the atom experiences a fictitious force given by F(t) = −ma(t), where m is the mass of the particle and a(t) is the acceleration of the trap as a function of time[64](https://www.nature.com/articles/s41586-022-04592-6#ref-CR64),[65](https://www.nature.com/articles/s41586-022-04592-6#ref-CR65). By following ref. [66](https://www.nature.com/articles/s41586-022-04592-6#ref-CR66) (equation 5.4), we find the average vibrational quantum number increase ΔN is given by

ΔN=|a~(ω0)|2(2xzpfω0)2,

(1)

where a~(ω0)

is the Fourier transform of a(t) evaluated at the trap frequency ω0, and the zero-point size of the particle xzpf≡ℏ/(2mω0)−−−−−−−−√, where ℏ is the reduced Planck constant. ΔN is the same for all initial levels of the oscillator[66](https://www.nature.com/articles/s41586-022-04592-6#ref-CR66). Experimentally, we apply an acceleration profile a(t) = jt to the atom, from time −T/2 to +T/2 to move a distance D with constant jerk j. We calculate |a~(ω)|2

, simplify using ω0T ≫ 1, and assume a small range of trap frequencies to average the oscillatory terms, resulting in

ΔN=12(6D/xzpfω20T2)2.

(2)

Several relevant insights can be gleaned from this formula. First, this expression indicates our ability to move large distances D with comparably small increases in time T. Furthermore, to maintain a constant ΔN, the movement time T∝ω−3/40

. Moreover, to perform a large number of moves k for a deep circuit, we can estimate ΔN ∝ k/T4, suggesting that we can increase our number of moves from, for example, 5 to 80 by slowing each move from 200 μs to 400 μs. Move speed could be further improved with different a(t) profiles, but inevitably with finite resources such as trap depth, quantum speed limits will eventually prevent arbitrarily fast motion of qubits across the array[30](https://www.nature.com/articles/s41586-022-04592-6#ref-CR30).

We now compare equation ([2](https://www.nature.com/articles/s41586-022-04592-6#Equ2)) with our experimental observations. In Fig. [1d](https://www.nature.com/articles/s41586-022-04592-6#Fig1) we start to observe atom loss when we move 55 μm in 200 μs under a constant negative jerk. This speed limit is consistent with our above estimates: using ω0 = 2π × 40 kHz and xzpf = 38 nm, we predict ΔN ≈ 6 for this move, corresponding to the onset of tangible heating at this move speed. More quantitatively, we assume a Poisson distribution with mean N and variance N and integrate the population above some critical Nmax upon which the atom will leave the trap. From this analysis we find atom retention is given by 12(1+erf[(Nmax−N)/2N−−−√])

.

Extended Data Fig. [2a, b](https://www.nature.com/articles/s41586-022-04592-6#Fig6) measures the atom retention as a function of move time T and trap frequency ω0/2π. Using the functional form above, for both sets of measurements, we extract an Nmax of about 30, corresponding to adding about 30 excitations before exciting the atom out of the trap. Such a limit is physically reasonable as the absolute trap depth of 4 MHz implies only about 100 levels, the atom starts at finite temperature, and moreover the effective trap frequency reduces once the anharmonicity of the trap starts to play a role. We note that these estimates are only approximate (we roughly estimate ω0 for the trap depths used during the motion), but nonetheless suggests our motion limit is consistent with physical limits for our chosen a(t). Our analysis here also neglects the acoustic lensing effects associated with ramping the AOD frequency, which causes astigmatism by focusing one axis to a different plane and thus deforms the trap and reduces the peak trap intensity (and ω0) as given by the Strehl ratio.

Additional heating and loss during the circuit can also be caused by repeated short drops for performing two-qubit gates, where the tweezers are briefly turned off to avoid anti-trapping of the Rydberg state and light shifts of the ground–Rydberg transition. However, drop–recapture measurements in Extended Data Fig. [2c](https://www.nature.com/articles/s41586-022-04592-6#Fig6) suggest that the 500-ns drops we use experimentally have a negligible effect until hundreds of drops per atom (corresponding to hundreds of CZ gates). We find that atom loss and heating as a function of number of drops are well described by a diffusion model, which would then predict that reducing atom temperature by a factor of 2× (reducing thermal velocity by 2–√×

) and reducing the drop time tdrop by 2×, together would increase the number of possible CZ gates per atom to thousands.

Two-qubit CZ gates implementation

We implement our two-qubit gates and calibrations following ref. [5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5). Specifically, the two-qubit CZ gate is implemented by two global Rydberg pulses, with each pulse at detuning Δ and length τ, and with a phase jump ξ between the two pulses. The pulse parameters are chosen such that qubit pairs, adjacent and under the Rydberg blockade constraint, will return from the Rydberg state back to the hyperfine qubit manifold with a phase depending on the state of the other qubit. The numerical values for these pulse parameters are:

Δ=−0.377371Ω

ξ=−0.621089×(2π)

τ=0.683201/[Ω/(2π)]

For our experiments in Figs. [1](https://www.nature.com/articles/s41586-022-04592-6#Fig1)–[3](https://www.nature.com/articles/s41586-022-04592-6#Fig3), we operate with a two-photon Rydberg Rabi frequency of Ω/2π = 3.6 MHz, giving a theoretical τ = 190 ns and a theoretical Δ/(2π) = −1.36 MHz. We choose the negative detuning sign to help minimize excitation into the mj = +1/2 Rydberg state (mj denotes magnetic sublevel of the 70S1/2 Rydberg state), which is detuned by about 24 MHz under the field of 8.5 G (and experiences a three-times lower coupling to the Rydberg laser than the desired mj = −1/2 state owing to reduced Clebsch–Gordan coefficients). In this work, we operate with strong blockade between adjacent qubits, with Rydberg–Rydberg interactions V0/2π ranging from 200 MHz to 1 GHz. In Fig. [4](https://www.nature.com/articles/s41586-022-04592-6#Fig4), we operate with Ω/2π = 4.45 MHz for the two-qubit gates.

Managing spurious phases during CZ gates

The two-qubit gate from ref. [5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5) induces both an intrinsic single-qubit phase, as well as spurious phases that are primarily induced by the differential light shift from the 420-nm laser. Under certain configurations, the 420-nm-induced differential light shift on the hyperfine qubit can be exceedingly large (>8 MHz), yielding phase accumulations on the hyperfine qubit of about 6π. Small, percent-level variations of the 420-nm intensity can thus lead to significant qubit dephasing.

Reference [5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5) addresses this 420-induced-phase issue by performing an echo sequence: after the CZ gate, the 1,013-nm Rydberg laser is turned off, a Raman π pulse is applied and then the 420-nm laser is pulsed again to cancel the phase induced by the 420 light during the CZ gate. This method echoes out the 420-induced phase, but comes at a cost of a factor of two increase in the 420-induced scattering error, which is the dominant source of error in our two-qubit CZ gates.

Echo between CZ gates

To address these various issues, here we perform a Raman π pulse between each CZ gate to echo out spurious gate-induced phases on the hyperfine qubit (Extended Data Fig. [1](https://www.nature.com/articles/s41586-022-04592-6#Fig5)). This approach has several advantages. The 420-induced phase is now cancelled by pairs of CZ gates, without explicitly applying additional 420-nm pulses to echo each individual CZ gate, thereby reducing the scattering error of the CZ gate in this work by a factor of approximately two. We estimate that this echo technique, having reduced the scattering error incurred during each gate, roughly compensates the increased scattering rate incurred by spreading our optical power over more space in 2D, thereby giving us comparable gate fidelites to the two-qubit CZ gate fidelities of ≥97.4(2)% reported in ref. [5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5). Furthermore, the echo between CZ gates also cancels the intrinsic single-qubit phase of the CZ gate, removing errors in the calibration of this parameter, as well as cancelling any other gate-induced spurious single-qubit phases such as a roughly 0.01-rad phase induced by pulsing the traps off for 500 ns for the two-qubit gate (Extended Data Fig. [1](https://www.nature.com/articles/s41586-022-04592-6#Fig5)). In instances where the number of CZ gates we apply is odd, we perform the echo for the final CZ gate.

Sign of intermediate-state detuning

To further suppress the effect of the spurious, 420-induced phase, we operate our 420-nm laser to be red-detuned (by 2 GHz) from the 6P3/2 transition. For red detunings, the light shift on the |0⟩ state and the |1⟩ state are of the same sign, minimizing the differential light shift, whereas for blue detunings <6.8 GHz, the light shift on the |0⟩ state and the |1⟩ state have opposite signs and amplify the differential light shift.

Sensitivity to axial trap oscillations

In typical Rydberg excitation timescales with optical tweezers, the axial trap oscillation frequencies of several kilohertz are inconsequential. Here with our circuits running as long as 1.2 ms, with Rydberg pulses throughout, we find that the axial trap oscillations can have important effects. In particular, the axial oscillations cause the atoms to make oscillations in and out of the Rydberg beams: at our estimated axial temperature of about 25 μK and axial oscillation frequency of 6 kHz, we estimate an axial spread ⟨z2⟩−−−√≈1.3μm

. For our 20-μm-waist beams, the effect of this positional spread is relatively small on the pulse parameters of the CZ gate, but can be significant on the sensitive 420-induced phase we seek to cancel by echoing out the phase induced by CZ gates separated by about 200 μs (see previous section). When using 20-μm-waist beams, and a 2.5-GHz blue detuning of our 420-nm laser, we find that the dephasing due to the axial trap oscillations is significant (Extended Data Fig. [4](https://www.nature.com/articles/s41586-022-04592-6#Fig8)). To remedy this deleterious effect, we increase the beam waist of our 420-nm laser to 35 μm (while maintaining constant intensity) and change the laser frequency to be 2-GHz red-detuned, together resulting in a significant reduction in the dephasing associated with improper echoing of the 420-nm pulse.

Bell-state preparation and fidelity

In Fig. [1](https://www.nature.com/articles/s41586-022-04592-6#Fig1), we prepare the |Φ+⟩ Bell state in the same way that is done in ref. [5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5). After initializing a pair of qubits in |00⟩, we apply X(π/2) pulse–CZ gate–X(π/4) pulse. We calculate and plot the raw resulting fidelity of this |Φ+⟩ Bell state as the sum of populations in |00⟩ and |11⟩, averaged with the fitted amplitude of parity oscillations (example in Fig. [1c](https://www.nature.com/articles/s41586-022-04592-6#Fig1)), which measures the off-diagonal coherences. In Fig. [1d](https://www.nature.com/articles/s41586-022-04592-6#Fig1), upon significant loss from movement, this fidelity estimate becomes skewed because we begin measuring an artificially large population in |11⟩ (as state |1⟩ is detected as loss); accordingly, we estimate the |Φ+⟩ population as two times the population of |00⟩ once the population difference between |11⟩ and |00⟩ becomes greater than 0.1 (an arbitrary cut-off where the effects of loss start to become significant). In Fig. [1d](https://www.nature.com/articles/s41586-022-04592-6#Fig1), for moves slower than 300 μs, we achieve an average raw Bell-state fidelity after the moving of 94.8(2)%. If we do not move or attempt to preserve coherence for 500 μs (that is, if we measure immediately after preparing the Bell state), then we measure a raw Bell-state fidelity of 95.2(1)% (not plotted here), consistent with the results in ref. [5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5).

Analysis of error sources

We detail here some of our measured and estimated sources of error for an entire sequence (toric code preparation in particular, our deepest circuit). We find the total single-qubit fidelity after performing the entire sequence is roughly 96.5% for the toric code circuit, which we measure by embedding the entire experiment in a Ramsey sequence: that is, we perform a Raman π/2 pulse, do all motion and decoupling, and then do a final π/2 pulse with variable phase to measure total contrast. We are able to account for our single-qubit fidelity quantitatively as being composed of our known single-qubit errors in Extended Data Fig. [6c](https://www.nature.com/articles/s41586-022-04592-6#Fig10).

Estimated contributions to two-qubit gate error are summarized in Extended Data Fig. [6c](https://www.nature.com/articles/s41586-022-04592-6#Fig10). These estimates come from numerical simulations in QuTiP (version 4.5.0) with experimental parameters. The effects of intermediate state scattering and Rydberg decay are included via collapse operators in the Lindblad master equation solver. Other error contributions include finite-temperature random Doppler shifts and position fluctuations, as well as laser pulse-to-pulse fluctuations, all of which are simulated using classical Monte Carlo sampling of experiment parameters. Experimental parameters used for the simulations are as follows: blue and red Rabi frequencies (Ωb, Ωr) = 2π × (160, 90) MHz, 6P3/2 intermediate state detuning of 2 GHz, intermediate state lifetime of 110 ns, 70S1/2 Rydberg state lifetime of 150 μs, Rydberg blockade energy of 500 MHz, splitting to second Rydberg state of 24 MHz, radial and axial trap frequencies (ωr, ωz) = 2π × (40, 6) kHz, and temperature T = 20 μK. We can also use this modelling to project for future performance; by assuming a 10 times increase in available 1,013-nm intensity and that atoms are cooled to a temperature of 2 uK, we project a possible CZ gate fidelity of 99.7%, beyond the surface code threshold[38](https://www.nature.com/articles/s41586-022-04592-6#ref-CR38),[67](https://www.nature.com/articles/s41586-022-04592-6#ref-CR67). Alkaline-earth atoms could also offer other routes to high-fidelity operations for QEC[68](https://www.nature.com/articles/s41586-022-04592-6#ref-CR68),[69](https://www.nature.com/articles/s41586-022-04592-6#ref-CR69),[70](https://www.nature.com/articles/s41586-022-04592-6#ref-CR70).

To understand how our various single-qubit and two-qubit errors contribute to our graph-state fidelities, we perform a stochastic simulation of the quantum circuit used for graph-state preparation (Extended Data Fig. [6a, b](https://www.nature.com/articles/s41586-022-04592-6#Fig10)). We utilize the Clifford properties of our circuit, allowing for efficient numerical evaluation and random sampling of many possible error realizations. The simulation is performed under a realistic error model, where the rates of ambient depolarizing noise and atom loss are measured in the experiment (Extended Data Fig. [6c](https://www.nature.com/articles/s41586-022-04592-6#Fig10)). The resulting stabilizer and logical qubit expectation values agree well with those measured experimentally.

Rydberg beam shaping and homogeneity

We shape our Rydberg beams into tophats of variable size through wavefront control using the phase profile on an SLM[25](https://www.nature.com/articles/s41586-022-04592-6#ref-CR25). This ability allows us to match the height of our beam profile to the experiment zone size of any given experiment, thereby maximizing our 1,013-nm light intensity and CZ gate fidelities. We optimize our Rydberg beam homogeneity until peak-to-peak inhomogenities are below <1%. To this end, we correct all aberrations up to the window of our vacuum chamber, as done in ref. [25](https://www.nature.com/articles/s41586-022-04592-6#ref-CR25), which yields an inhomogeneity on the atoms of several per cent that we attribute to imperfections of the final window. To further optimize the homogeneity, we empirically tune aberration corrections on the tophat through Zernike polynomial corrections to the phase profile in the SLM plane (Fourier plane). With this procedure, we reduce peak-to-peak inhomogeneities to <1% over a range of 40–50 μm in the atom plane.

Creation and optimization of graph layouts

We outline here a description of how we optimize our graph layouts for the cluster state, Steane code, surface code and toric code preparation. Our optimization in this work is heuristic, and future work can develop appropriate algorithms for designing optimal circuits through atom spatial arrangement and AOD trajectories. Extended Data Fig. [5](https://www.nature.com/articles/s41586-022-04592-6#Fig9) shows all of the graphs we create and the process for creating them. There are several parameters we optimize for. (1) Minimize the number of parallel two-qubit gate layers. (2) Minimize the total move distance for the moving atoms. (3) Have all moving atoms in one sublattice (all graphs realized here are bipartite) to facilitate the final local rotation of one sublattice. (4) Minimize the vertical extent of the array and the number of distinct rows (to maximize 1,013-nm light intensity and minimize sensitivity to beam inhomogeneity between the rows). (5) When ordering gates, apply two-qubit gates as early as possible in the circuit. If a gate layer induces a bit-flip (X error) then that error can propagate during subsequent gates (becoming a Z error on the other qubit), so gates should be in the earliest layer possible.

Local (sublattice) hyperfine rotations

We perform local rotations in the hyperfine basis by use of our horizontally propagating 420-nm beam, which imposes a differential light of several megahertz on the hyperfine qubit and can thus be used for realizing a fast Z rotation. To realize the local Y(π/2) rotation used throughout this work, we move one sublattice of atoms out of the 420-nm beam, then apply [global Y(π/4)]–[local Z(π)]–[global Y(π/4)]. This realizes a Y(π/2) rotation on one sublattice and a Z(π) rotation on the other sublattice (which is inconsequential as it then commutes with the immediately following measurement in the Z basis). To apply a Y(π/2) on the other sublattice of atoms, we add an additional global Z(π) (implemented by jumping the Raman laser phase) between the two Y(π/4) pulses. Future experiments will benefit from an additional set of locally focused beams for performing local Raman control of hyperfine qubit states, but we find that moving atoms works so efficiently (even for moving >50 μm to move out of the 420-nm beam) that this approach is well suited for producing a high-fidelity, homogeneous rotation on roughly half the qubits.

Local Rydberg initialization

We perform local Rydberg control to initialize the |ℤ2⟩ = |rgrg...⟩ ≡ |r1r1...⟩ state for studying the dynamics of many-body scars. We achieve this local initialization by applying approximately 50-MHz light shifts between |1⟩ and |r⟩ using 810-nm tweezers generated by an SLM onto a desired subset of sites, and then apply a global Rydberg π pulse, which excites the non-light-shifted atoms. We use this approach here to prepare every other atom in each chain into |r⟩, but emphasize that as the locations of the SLM tweezers are fully programmable, this technique can be used to prepare any initial blockade-satisfying configuration of atoms in |1⟩ and |r⟩.

The 50-MHz biasing light shift is significantly larger than the Rydberg Rabi frequency Ω/2π = 4.45 MHz, leading to a Rydberg population on undesired sites of <1%. The t = 0 time point of Extended Data Fig. [10b](https://www.nature.com/articles/s41586-022-04592-6#Fig14) shows the high-fidelity preparation of the |ℤ2⟩ state using this approach. We note that with 810-nm light, even though the achieved biasing light shift is significant, the Raman-scattering-induced T1 (of the hyperfine qubit) is still about 50 ms and thus leads to a scattering error ≲4 × 10−6 during the 200-ns pulse of the light-shifting tweezers. There can also be a motional effect from the biasing tweezers, with an estimated radial trapping frequency of 150 kHz, which we also deem to be negligible during the 200-ns pulse.

Rydberg Hamiltonian

In Fig. [4](https://www.nature.com/articles/s41586-022-04592-6#Fig4), we study dynamics under the many-body Rydberg Hamiltonian

HRydℏ=Ω2∑iσxi−Δ∑ini+∑i<jVijninj,

(3)

where ℏ

is the reduced Planck constant, Ω is the Rabi frequency, Δ is the laser frequency detuning, ni = |ri⟩⟨ri| is the projector onto the Rydberg state at site i and σxi=|1i⟩⟨ri|+|ri⟩⟨1i|

flips the atomic state. For the entanglement entropy measurements in this work, we choose lattice spacings where the nearest-neighbour interaction V0 > Ω results in the Rydberg blockade, preventing adjacent atoms from simultaneously occupying |r⟩. In particular, the many-body experiments are performed on eight-atom chains, quenching to a time-independent HRyd with V0/2π = 20 MHz, Ω/2π = 3.1 MHz and Δ/2π = 0.3 MHz. Quenching to small, positive Δ = 0.0173V0 partially suppresses the always-positive long-range interactions and thereby is optimal for scar lifetime, as derived and shown experimentally in ref. [71](https://www.nature.com/articles/s41586-022-04592-6#ref-CR71).

Coherent mapping protocol

As described in the text, we implement a coherent mapping protocol to transfer a generic many-body state in the {|1⟩, |r⟩} basis to the long-lived and non-interacting {|0⟩, |1⟩} basis. To achieve this mapping, immediately following the Rydberg dynamics we apply a Raman π pulse to map |1⟩ → |0⟩, and then a subsequent Rydberg π pulse to map |r⟩ → |1⟩ (ref. [72](https://www.nature.com/articles/s41586-022-04592-6#ref-CR72)).

Even for perfect Raman and Rydberg π pulses (on isolated atoms), there are three key sources of infidelity associated with this mapping process. (1) Any population in blockade-violating states (that is, two adjacent atoms both in |r⟩) will be strongly shifted off-resonance for the final Rydberg π pulse. As such, this atomic population will be left in the Rydberg state and lost. (2) Long-range interactions, for example, from next-nearest neighbours, will detune the final Rydberg π pulse from resonance and thus reduce pulse fidelity. As the long-range interactions are not the same for all many-body microstates, this effect cannot be mitigated by a simple shift of the detuning. (3) Dephasing of the state occurs throughout the duration of the Raman π pulse, predominantly from Doppler shifts between the ground states |0⟩ and |1⟩ and the Rydberg state |r⟩. Although these random on-site detunings are also present during the many-body dynamics, turning the Rydberg drive Ω off allows the system to freely accumulate phase and makes us particularly sensitive to dephasing errors.

We now detail our mitigation of the above error mechanisms. To minimize errors from (1), we perform our many-body dynamics with Ω2/(2V20)≈0.01

. This minimizes the probability of an atom to violate blockade to be of order 1%. To help minimize errors from (2), we increase the amplitude of the 420-nm laser for the final π pulse by a factor of 2×, such that (VNNN/Ω)2 = 0.005 (where VNNN is the interactions with next-nearest neighbours), reducing pulse errors from long-range interactions to order 1%. Finally, to reduce errors from (3), we perform a fast Raman π pulse and leave only 150 ns between ending the many-body Rydberg dynamics and beginning the Rydberg π pulse. The 150-ns gap is comparably short relative to the T∗2

 ≈ 3–4 μs of the {|g⟩, |r⟩} basis, leading to a random phase accumulation of the order of about 0.02 × 2π rad per particle, but is further compounded by having entangled states of N particles in one copy accumulating a random phase relative to entangled states of N particles in the second copy. We study these various effects numerically in Extended Data Fig. [9c](https://www.nature.com/articles/s41586-022-04592-6#Fig13).

Finally, we note that the global Raman beam induces a light-shift-induced phase shift of about π on |0⟩ and |1⟩ relative to |r⟩ during the Raman π pulse. Similarly, the global 420-nm laser also induces a light-shift-induced phase shift of about π between |0⟩ and |1⟩ during the Rydberg π pulse. Although the measurements we perform here are interferometric (in other words, the singlet state we measure is invariant under global rotations) and thus not affected by these global phase shifts, in future work these phase shifts can be measured and accounted for where relevant.

Measuring entanglement entropy

The second-order Renyi entanglement entropy is given by S2(A)=

−log2Tr[ρ2A], where Tr[ρ2A] is the state purity of reduced density matrix ρA on subsystem A. The purity can be measured with two copies by noticing that Tr[ρ2A]=Tr[S^ρA⊗ρA] is the expectation value of the many-body SWAP operator S^ (ref. [10](https://www.nature.com/articles/s41586-022-04592-6#ref-CR10),[11](https://www.nature.com/articles/s41586-022-04592-6#ref-CR11)). The many-body SWAP operator is composed of individual SWAP operators s^i on each twin pair, that is, S^=Πis^i (with i ∈ A). Measuring this expectation value amounts to probing occurences of the singlet state |01⟩−|10⟩2√ (with eigenvalue −1 under s^i), as all other s^i eigenstates have eigenvalue +1. Occurences of the singlet state in each twin pair, that is, the Bell state |Ψ−⟩, is extracted by a Bell measurement circuit (with an additional local Z(π), see next paragraph), which maps |Ψ−⟩ → |00⟩ and can thereafter be measured in the computational basis. As such, after performing the Bell measurement circuit, we analyse the resulting bit-string outputs and calculate the purity of any subsystem A by calculating ⟨Πi∈As^i⟩

: that is, we measure purity as the average parity = ⟨(−1)observed |00⟩ pairs⟩ within A. In the absence of experimental imperfections, the purity will equal 1 for the whole system, and be less than 1 for subsystems that are entangled with the rest of the system.

A Bell measurement circuit can be decomposed into applying an X(π/2) rotation on one atom of the twin pair, then applying a CZ gate and then a global X(π/2) rotation. In other measurements, we realize a local X(π/2) by doing a global X(π/4) rotation, then local Z(π) rotation and then global X(π/4). However, we note that for this singlet measurement circuit, the first X(π/4) is redundant as the singlet state is invariant under global rotations, and so for the local X(π/2) we only apply the local Z(π) and then the second global X(π/4). This effectively realizes the X(π/2) on one qubit, up to a Z(π) on the other qubit (not shown in circuit diagram in Fig. [4](https://www.nature.com/articles/s41586-022-04592-6#Fig4)). Under this simplification, the Bell measurement circuit to map |Ψ−⟩ → |00⟩ can be roughly understood as the reverse of the Bell-state preparation circuit in ref. [5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5), which is precisely how we calibrate the parameters of the Bell measurement.

Calibrating and benchmarking the interferometry

To validate the interferometry measurement (and check for proper calibration), we benchmark it separately from the many-body dynamics and coherent mapping protocol. We perform this benchmarking by preparing independent qubits in identical, variable single-qubit superpositions (through a global Raman pulse of variable time) and ensuring that the interferometry rarely results in |00⟩ for all the variable initial product states (Extended Data Fig. [8a](https://www.nature.com/articles/s41586-022-04592-6#Fig12)). We find this is an important benchmarking step, because we find that small miscalibrations of the Bell measurement can lead to lower fidelity (that is, higher entropy) for different initial product states and thereby result in additional spurious signals in an entanglement entropy measurement. We note that this measurement is particularly sensitive to the single-qubit phase immediately before the final X(π/2) pulse (induced by the CZ gate and cancelled by a global Z(θ) pulse).

Additional many-body data and details

To benchmark our method of measuring entanglement entropy in a many-body system, in Extended Data Fig. [8b](https://www.nature.com/articles/s41586-022-04592-6#Fig12) we study the entanglement dynamics after initializing two proximal atoms in |1⟩ and resonantly exciting to the Rydberg state for a variable time t. Under conditions of Rydberg blockade, this excitation results in two-particle Rabi oscillations between |11⟩ and the entangled state |W⟩=12√(|1r⟩+|r1⟩)

(top panel of Extended Data Fig. [8b](https://www.nature.com/articles/s41586-022-04592-6#Fig12))[3](https://www.nature.com/articles/s41586-022-04592-6#ref-CR3),[13](https://www.nature.com/articles/s41586-022-04592-6#ref-CR13),[72](https://www.nature.com/articles/s41586-022-04592-6#ref-CR72). The state purity of this two-particle system is measured by performing Bell measurements on atom pairs from two identical copies. Locally, the measured purity of the one-particle subsystem reduces to a value of about 0.5 when the system enters the maximally entangled |W⟩ state, at which point the reduced density matrix of each individual atom is maximally mixed. In contrast, the purity of the global, two-particle state remains high. The observation that the global-state purity is higher than the local-subsystem purity is a distinct signature of quantum entanglement[11](https://www.nature.com/articles/s41586-022-04592-6#ref-CR11),[12](https://www.nature.com/articles/s41586-022-04592-6#ref-CR12).

For the data shown in Fig. [4c, e](https://www.nature.com/articles/s41586-022-04592-6#Fig4), we subtract the data by an extensive classical entropy as is done in ref. [12](https://www.nature.com/articles/s41586-022-04592-6#ref-CR12). This fixed, time-independent offset is given by the entropy per particle, that is, (global entropy at quench time t = 0) × (subsystem size)/(global system size). In Extended Data Fig. [9a](https://www.nature.com/articles/s41586-022-04592-6#Fig13), we show the raw entanglement entropy measurements alongside numerics, to indicate the size of the extensive classical entropy contribution. In plotting, we also delay the theory curves by 10 ns to account for the fact that the Raman π pulse cuts off the final 10 ns of the Rydberg evolution, which is done to keep the coherent mapping gap as short as possible and minimize Doppler dephasing. Furthermore, in Extended Data Fig. [9b](https://www.nature.com/articles/s41586-022-04592-6#Fig13) we plot the measured global purity and compare it with numerical simulations incorporating experimental errors (Extended Data Fig. [9c](https://www.nature.com/articles/s41586-022-04592-6#Fig13)).

In Extended Data Fig. [10](https://www.nature.com/articles/s41586-022-04592-6#Fig14), we show additional many-body data on the eight-atom chain system, with the same parameters as those used in the main text. We show the measured single-site entropy of each site in the eight-atom chain for the |ℤ2⟩ quench in Extended Data Fig. [10a](https://www.nature.com/articles/s41586-022-04592-6#Fig14). Furthermore, in Extended Data Fig. [10b](https://www.nature.com/articles/s41586-022-04592-6#Fig14), we plot the global Rydberg population, measured in both the {|1⟩, |r⟩} basis and the {|0⟩, |1⟩} basis.

Data availability

The data that supports the findings of this study are available from the corresponding author on reasonable request.

References

-

Cirac, J. I. & Zoller, P. A scalable quantum computer with ions in an array of microtraps. Nature 404, 579–581 (2000).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BD3cXis1GmtL8%3D) [Article](https://doi.org/10.1038%2F35007021) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=A%20scalable%20quantum%20computer%20with%20ions%20in%20an%20array%20of%20microtraps&journal=Nature&volume=404&pages=579-581&publication_year=2000&author=Cirac%2CJI&author=Zoller%2CP)

-

Preskill, J. Quantum Computing in the NISQ era and beyond. Quantum 2, 79 (2018).

[Article](https://doi.org/10.22331%2Fq-2018-08-06-79) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Quantum%20Computing%20in%20the%20NISQ%20era%20and%20beyond&journal=Quantum&volume=2&publication_year=2018&author=Preskill%2CJ)

-

Jaksch, D. et al. Fast quantum gates for neutral atoms. Phys. Rev. Lett. 85, 2208–2211 (2000).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BD3cXmt1Wku7g%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.85.2208) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Fast%20quantum%20gates%20for%20neutral%20atoms&journal=Phys.%20Rev.%20Lett.&volume=85&pages=2208-2211&publication_year=2000&author=Jaksch%2CD)

-

Urban, E. et al. Observation of Rydberg blockade between two atoms. Nat. Phys. 5, 110–114 (2009).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BD1MXhtlKktr8%3D) [Article](https://doi.org/10.1038%2Fnphys1178) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Observation%20of%20Rydberg%20blockade%20between%20two%20atoms&journal=Nat.%20Phys.&volume=5&pages=110-114&publication_year=2009&author=Urban%2CE)

-

Levine, H. et al. Parallel implementation of high-fidelity multiqubit gates with neutral atoms. Phys. Rev. Lett. 123, 170503 (2019).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC1MXit1ent7rP) [Article](https://doi.org/10.1103%2FPhysRevLett.123.170503) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Parallel%20implementation%20of%20high-fidelity%20multiqubit%20gates%20with%20neutral%20atoms&journal=Phys.%20Rev.%20Lett.&volume=123&publication_year=2019&author=Levine%2CH)

-

Nigg, D. et al. Quantum computations on a topologically encoded qubit. Science 345, 302–305 (2014).

[MathSciNet](http://www.ams.org/mathscinet-getitem?mr=3235634) [CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC2cXhtFCqs7nJ) [MATH](http://www.emis.de/MATH-item?1355.81061) [Article](https://doi.org/10.1126%2Fscience.1253742) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Quantum%20computations%20on%20a%20topologically%20encoded%20qubit&journal=Science&volume=345&pages=302-305&publication_year=2014&author=Nigg%2CD)

-

Ryan-Anderson, C. et al. Realization of real-time fault-tolerant quantum error correction. Phys. Rev. 11, 041058 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB38XhtV2qsrk%3D) [Article](https://doi.org/10.1103%2FPhysRevX.11.041058) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Realization%20of%20real-time%20fault-tolerant%20quantum%20error%20correction&journal=Phys.%20Rev.&volume=11&publication_year=2021&author=Ryan-Anderson%2CC)

-

Satzinger, K. J. et al. Realizing topologically ordered states on a quantum processor. Science 374, 1237–1241 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXislGjsrbK) [Article](https://doi.org/10.1126%2Fscience.abi8378) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Realizing%20topologically%20ordered%20states%20on%20a%20quantum%20processor&journal=Science&volume=374&pages=1237-1241&publication_year=2021&author=Satzinger%2CKJ)

-

Kitaev, A. Y. Fault-tolerant quantum computation by anyons. Ann. Phys. 303, 2–30 (2003).

[MathSciNet](http://www.ams.org/mathscinet-getitem?mr=1951039) [CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BD3sXisl2hsg%3D%3D) [MATH](http://www.emis.de/MATH-item?1012.81006) [Article](https://doi.org/10.1016%2FS0003-4916%2802%2900018-0) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Fault-tolerant%20quantum%20computation%20by%20anyons&journal=Ann.%20Phys.&volume=303&pages=2-30&publication_year=2003&author=Kitaev%2CAY)

-

Daley, A. J., Pichler, H., Schachenmayer, J. & Zoller, P. Measuring entanglement growth in quench dynamics of bosons in an optical lattice. Phys. Rev. Lett. 109, 020505 (2012).

[CAS](https://www.nature.com/articles/cas-redirect/1:STN:280:DC%2BC3s%2FhvFCqtw%3D%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.109.020505) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Measuring%20entanglement%20growth%20in%20quench%20dynamics%20of%20bosons%20in%20an%20optical%20lattice&journal=Phys.%20Rev.%20Lett.&volume=109&publication_year=2012&author=Daley%2CAJ&author=Pichler%2CH&author=Schachenmayer%2CJ&author=Zoller%2CP)

-

Islam, R. et al. Measuring entanglement entropy in a quantum many-body system. Nature 528, 77–83 (2015).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC2MXhvFKht7jI) [Article](https://doi.org/10.1038%2Fnature15750) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Measuring%20entanglement%20entropy%20in%20a%20quantum%20many-body%20system&journal=Nature&volume=528&pages=77-83&publication_year=2015&author=Islam%2CR)

-

Kaufman, A. M. et al. Quantum thermalization through entanglement in an isolated many-body system. Science 353, 794–800 (2016).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC28XhtlCju73M) [Article](https://doi.org/10.1126%2Fscience.aaf6725) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Quantum%20thermalization%20through%20entanglement%20in%20an%20isolated%20many-body%20system&journal=Science&volume=353&pages=794-800&publication_year=2016&author=Kaufman%2CAM)

-

Bernien, H. et al. Probing many-body dynamics on a 51-atom quantum simulator. Nature 551, 579–584 (2017).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC2sXhvFWht7fL) [Article](https://doi.org/10.1038%2Fnature24622) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Probing%20many-body%20dynamics%20on%20a%2051-atom%20quantum%20simulator&journal=Nature&volume=551&pages=579-584&publication_year=2017&author=Bernien%2CH)

-

Turner, C. J., Michailidis, A. A., Abanin, D. A., Serbyn, M. & Papić, Z. Weak ergodicity breaking from quantum many-body scars. Nat. Phys. 14, 745–749 (2018).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC1cXhtFGktr3J) [Article](https://doi.org/10.1038%2Fs41567-018-0137-5) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Weak%20ergodicity%20breaking%20from%20quantum%20many-body%20scars&journal=Nat.%20Phys.&volume=14&pages=745-749&publication_year=2018&author=Turner%2CCJ&author=Michailidis%2CAA&author=Abanin%2CDA&author=Serbyn%2CM&author=Papi%C4%87%2CZ)

-

Molmer, K. & Sorensen, A. Multiparticle entanglement of hot trapped ions. Phys. Rev. Lett. 82, 1835 (1999).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DyaK1MXhtlersLk%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.82.1835) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Multiparticle%20entanglement%20of%20hot%20trapped%20ions&journal=Phys.%20Rev.%20Lett.&volume=82&publication_year=1999&author=Molmer%2CK&author=Sorensen%2CA)

-

Wright, K. et al. Benchmarking an 11-qubit quantum computer. Nat. Commun. 10, 5464 (2019).

[CAS](https://www.nature.com/articles/cas-redirect/1:STN:280:DC%2BB3MfkvVylug%3D%3D) [Article](https://doi.org/10.1038%2Fs41467-019-13534-2) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Benchmarking%20an%2011-qubit%20quantum%20computer&journal=Nat.%20Commun.&volume=10&publication_year=2019&author=Wright%2CK)

-

Periwal, A. et al. Programmable interactions and emergent geometry in an array of atom clouds. Nature 600, 630–635 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXivValtrvK) [Article](https://doi.org/10.1038%2Fs41586-021-04156-0) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Programmable%20interactions%20and%20emergent%20geometry%20in%20an%20array%20of%20atom%20clouds&journal=Nature&volume=600&pages=630-635&publication_year=2021&author=Periwal%2CA)

-

Majer, J. et al. Coupling superconducting qubits via a cavity bus. Nature 449, 443–447 (2007).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BD2sXhtV2is7%2FO) [Article](https://doi.org/10.1038%2Fnature06184) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Coupling%20superconducting%20qubits%20via%20a%20cavity%20bus&journal=Nature&volume=449&pages=443-447&publication_year=2007&author=Majer%2CJ)

-

Mandel, O. et al. Coherent transport of neutral atoms in spin-dependent optical lattice potentials. Phys. Rev. Lett. 91, 010407 (2003).

[Article](https://doi.org/10.1103%2FPhysRevLett.91.010407) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Coherent%20transport%20of%20neutral%20atoms%20in%20spin-dependent%20optical%20lattice%20potentials&journal=Phys.%20Rev.%20Lett.&volume=91&publication_year=2003&author=Mandel%2CO)

-

Beugnon, J. et al. Two-dimensional transport and transfer of a single atomic qubit in optical tweezers. Nat. Phys. 3, 696–699 (2007).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BD2sXhtFSgsrfO) [Article](https://doi.org/10.1038%2Fnphys698) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Two-dimensional%20transport%20and%20transfer%20of%20a%20single%20atomic%20qubit%20in%20optical%20tweezers&journal=Nat.%20Phys.&volume=3&pages=696-699&publication_year=2007&author=Beugnon%2CJ)

-

Monroe, C. et al. Large-scale modular quantum-computer architecture with atomic memory and photonic interconnects. Phys. Rev. A 89, 022317 (2014).

[Article](https://doi.org/10.1103%2FPhysRevA.89.022317) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Large-scale%20modular%20quantum-computer%20architecture%20with%20atomic%20memory%20and%20photonic%20interconnects&journal=Phys.%20Rev.%20A&volume=89&publication_year=2014&author=Monroe%2CC)

-

Wu, T.-Y., Kumar, A., Giraldo, F. & Weiss, D. S. Stern–Gerlach detection of neutral-atom qubits in a state-dependent optical lattice. Nat. Phys. 15, 538–542 (2019).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC1MXms1Omtbw%3D) [Article](https://doi.org/10.1038%2Fs41567-019-0478-8) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Stern%E2%80%93Gerlach%20detection%20of%20neutral-atom%20qubits%20in%20a%20state-dependent%20optical%20lattice&journal=Nat.%20Phys.&volume=15&pages=538-542&publication_year=2019&author=Wu%2CT-Y&author=Kumar%2CA&author=Giraldo%2CF&author=Weiss%2CDS)

-

Pompili, M. et al. Realization of a multinode quantum network of remote solid-state qubits. Science 372, 259–264 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXovVKksrk%3D) [Article](https://doi.org/10.1126%2Fscience.abg1919) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Realization%20of%20a%20multinode%20quantum%20network%20of%20remote%20solid-state%20qubits&journal=Science&volume=372&pages=259-264&publication_year=2021&author=Pompili%2CM)

-

Pino, J. M. et al. Demonstration of the trapped-ion quantum CCD computer architecture. Nature 592, 209–213 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXosVClsLg%3D) [Article](https://doi.org/10.1038%2Fs41586-021-03318-4) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Demonstration%20of%20the%20trapped-ion%20quantum%20CCD%20computer%20architecture&journal=Nature&volume=592&pages=209-213&publication_year=2021&author=Pino%2CJM)

-

Ebadi, S. et al. Quantum phases of matter on a 256-atom programmable quantum simulator. Nature 595, 227–232 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXhsFSnsbbK) [Article](https://doi.org/10.1038%2Fs41586-021-03582-4) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Quantum%20phases%20of%20matter%20on%20a%20256-atom%20programmable%20quantum%20simulator&journal=Nature&volume=595&pages=227-232&publication_year=2021&author=Ebadi%2CS)

-

Levine, H. et al. Dispersive optics for scalable Raman driving of hyperfine qubits. Phys. Rev. A 105, 032618 (2022).

-

Vandersypen, L. M. K. & Chuang, I. L. NMR techniques for quantum control and computation. Rev. Mod. Phys. 76, 1037 (2005).

[Article](https://doi.org/10.1103%2FRevModPhys.76.1037) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=NMR%20techniques%20for%20quantum%20control%20and%20computation&journal=Rev.%20Mod.%20Phys.&volume=76&publication_year=2005&author=Vandersypen%2CLMK&author=Chuang%2CIL)

-

Gullion, T., Baker, D. B. & Conradi, M. S. New, compensated Carr–Purcell sequences. J. Magn. Reson. 89, 479–484 (1990).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DyaK3MXhs1ejug%3D%3D) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=New%2C%20compensated%20Carr%E2%80%93Purcell%20sequences&journal=J.%20Magn.%20Reson.&volume=89&pages=479-484&publication_year=1990&author=Gullion%2CT&author=Baker%2CDB&author=Conradi%2CMS)

-

Dordević, T. et al. Entanglement transport and a nanophotonic interface for atoms in optical tweezers. Science 373, 1511–1514 (2021).

[MathSciNet](http://www.ams.org/mathscinet-getitem?mr=4377593) [Article](https://doi.org/10.1126%2Fscience.abi9917) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Entanglement%20transport%20and%20a%20nanophotonic%20interface%20for%20atoms%20in%20optical%20tweezers&journal=Science&volume=373&pages=1511-1514&publication_year=2021&author=Dordevi%C4%87%2CT)

-

Lam, M. R. et al. Demonstration of quantum brachistochrones between distant states of an atom. Phys. Rev. X 11, 011035 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXovVCmsbo%3D) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Demonstration%20of%20quantum%20brachistochrones%20between%20distant%20states%20of%20an%20atom&journal=Phys.%20Rev.%20X&volume=11&publication_year=2021&author=Lam%2CMR)

-

Hein, M. et al. Entanglement in graph states and its applications. In Proc. International School of Physics “Enrico Fermi” (eds Casati, G. et al.) Vol. 162, 115–218 (IOS Press, 2006).

-

Raussendorf, R. & Briegel, H. J. A one-way quantum computer. Phys. Rev. Lett. 86, 5188 (2001).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BD3MXjvFWmsL4%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.86.5188) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=A%20one-way%20quantum%20computer&journal=Phys.%20Rev.%20Lett.&volume=86&publication_year=2001&author=Raussendorf%2CR&author=Briegel%2CHJ)

-

Schwartz, I. et al. Deterministic generation of a cluster state of entangled photons. Science 354, 434–437 (2016).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC28Xhslamt7jM) [Article](https://doi.org/10.1126%2Fscience.aah4758) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Deterministic%20generation%20of%20a%20cluster%20state%20of%20entangled%20photons&journal=Science&volume=354&pages=434-437&publication_year=2016&author=Schwartz%2CI)

-

Tóth, G. & Gühne, O. Entanglement detection in the stabilizer formalism. Phys. Rev. A 72, 022340 (2005).

[Article](https://doi.org/10.1103%2FPhysRevA.72.022340) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Entanglement%20detection%20in%20the%20stabilizer%20formalism&journal=Phys.%20Rev.%20A&volume=72&publication_year=2005&author=T%C3%B3th%2CG&author=G%C3%BChne%2CO)

-

Raussendorf, R. & Harrington, J. Fault-tolerant quantum computation with high threshold in two dimensions. Phys. Rev. Lett. 98, 190504 (2007).

[Article](https://doi.org/10.1103%2FPhysRevLett.98.190504) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Fault-tolerant%20quantum%20computation%20with%20high%20threshold%20in%20two%20dimensions&journal=Phys.%20Rev.%20Lett.&volume=98&publication_year=2007&author=Raussendorf%2CR&author=Harrington%2CJ)

-

Gottesman, D. An introduction to quantum error correction and fault-tolerant quantum computation. Preprint at https://arxiv.org/abs/0904.2557 (2009).

-

Egan, L. et al. Fault-tolerant control of an error-corrected qubit. Nature 598, 281–286 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXitFOltLvK) [Article](https://doi.org/10.1038%2Fs41586-021-03928-y) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Fault-tolerant%20control%20of%20an%20error-corrected%20qubit&journal=Nature&volume=598&pages=281-286&publication_year=2021&author=Egan%2CL)

-

Fowler, A. G., Mariantoni, M., Martinis, J. M. & Cleland, A. N. Surface codes: towards practical large-scale quantum computation. Phys. Rev. A 86, 032324 (2012).

[Article](https://doi.org/10.1103%2FPhysRevA.86.032324) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Surface%20codes%3A%20towards%20practical%20large-scale%20quantum%20computation&journal=Phys.%20Rev.%20A&volume=86&publication_year=2012&author=Fowler%2CAG&author=Mariantoni%2CM&author=Martinis%2CJM&author=Cleland%2CAN)

-

Semeghini, G. et al. Probing topological spin liquids on a programmable quantum simulator. Science 374, 1242–1247 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXislGjsrbL) [Article](https://doi.org/10.1126%2Fscience.abi8794) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Probing%20topological%20spin%20liquids%20on%20a%20programmable%20quantum%20simulator&journal=Science&volume=374&pages=1242-1247&publication_year=2021&author=Semeghini%2CG)

-

Bravyi, S., Hastings, M. B. & Verstraete, F. Lieb–Robinson bounds and the generation of correlations and topological quantum order. Phys. Rev. Lett. 97, 050401 (2006).

[CAS](https://www.nature.com/articles/cas-redirect/1:STN:280:DC%2BD28rpsVWgtQ%3D%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.97.050401) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Lieb%E2%80%93Robinson%20bounds%20and%20the%20generation%20of%20correlations%20and%20topological%20quantum%20order&journal=Phys.%20Rev.%20Lett.&volume=97&publication_year=2006&author=Bravyi%2CS&author=Hastings%2CMB&author=Verstraete%2CF)

-

Bolt, A., Duclos-Cianci, G., Poulin, D. & Stace, T. Foliated quantum error-correcting codes. Phys. Rev. Lett. 117, 070501 (2016).

[MathSciNet](http://www.ams.org/mathscinet-getitem?mr=3630926) [CAS](https://www.nature.com/articles/cas-redirect/1:STN:280:DC%2BC2szkvVaksA%3D%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.117.070501) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Foliated%20quantum%20error-correcting%20codes&journal=Phys.%20Rev.%20Lett.&volume=117&publication_year=2016&author=Bolt%2CA&author=Duclos-Cianci%2CG&author=Poulin%2CD&author=Stace%2CT)

-

Wilk, T. et al. Entanglement of two individual neutral atoms using Rydberg blockade. Phys. Rev. Lett. 104, 010502 (2010).

[CAS](https://www.nature.com/articles/cas-redirect/1:STN:280:DC%2BC3c3jslOisg%3D%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.104.010502) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Entanglement%20of%20two%20individual%20neutral%20atoms%20using%20Rydberg%20blockade&journal=Phys.%20Rev.%20Lett.&volume=104&publication_year=2010&author=Wilk%2CT)

-

Lesanovsky, I. & Katsura, H. Interacting Fibonacci anyons in a Rydberg gas. Phys. Rev. A 86, 041601 (2012).

[Article](https://doi.org/10.1103%2FPhysRevA.86.041601) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Interacting%20Fibonacci%20anyons%20in%20a%20Rydberg%20gas&journal=Phys.%20Rev.%20A&volume=86&publication_year=2012&author=Lesanovsky%2CI&author=Katsura%2CH)

-

Kaufman, A. M., Lester, B. J. & Regal, C. A. Cooling a single atom in an optical tweezer to its quantum ground state. Phys. Rev. 2, 041014 (2012).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC3sXit1Cmt7Y%3D) [Article](https://doi.org/10.1103%2FPhysRevX.2.041014) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Cooling%20a%20single%20atom%20in%20an%20optical%20tweezer%20to%20its%20quantum%20ground%20state&journal=Phys.%20Rev.&volume=2&publication_year=2012&author=Kaufman%2CAM&author=Lester%2CBJ&author=Regal%2CCA)

-

Shea, M. E., Baker, P. M., Joseph, J. A., Kim, J. & Gauthier, D. J. Submillisecond, nondestructive, time-resolved quantum-state readout of a single, trapped neutral atom. Phys. Rev. A 102, 053101 (2020).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3cXisFylsLrM) [Article](https://doi.org/10.1103%2FPhysRevA.102.053101) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Submillisecond%2C%20nondestructive%2C%20time-resolved%20quantum-state%20readout%20of%20a%20single%2C%20trapped%20neutral%20atom&journal=Phys.%20Rev.%20A&volume=102&publication_year=2020&author=Shea%2CME&author=Baker%2CPM&author=Joseph%2CJA&author=Kim%2CJ&author=Gauthier%2CDJ)

-

Cong, I., Levine, H., Keesling, A., Bluvstein, D., Wang, S.-T. & Lukin, M. D. Hardware-efficient, fault-tolerant quantum computation with Rydberg atoms. Preprint at https://arxiv.org/abs/2105.13501 (2021).

-

Wang, C., Harrington, J. & Preskill, J. Confinement-Higgs transition in a disordered gauge theory and the accuracy threshold for quantum memory. Ann. Phys. 303, 31–58 (2003).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BD3sXisl2hsw%3D%3D) [MATH](http://www.emis.de/MATH-item?1023.81024) [Article](https://doi.org/10.1016%2FS0003-4916%2802%2900019-2) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Confinement-Higgs%20transition%20in%20a%20disordered%20gauge%20theory%20and%20the%20accuracy%20threshold%20for%20quantum%20memory&journal=Ann.%20Phys.&volume=303&pages=31-58&publication_year=2003&author=Wang%2CC&author=Harrington%2CJ&author=Preskill%2CJ)

-

Bombín, H. Gauge color codes: optimal transversal gates and gauge fixing in topological stabilizer codes. New J. Phys. 17, 083002 (2015).

[MathSciNet](http://www.ams.org/mathscinet-getitem?mr=3404220) [MATH](http://www.emis.de/MATH-item?1454.81052) [Article](https://doi.org/10.1088%2F1367-2630%2F17%2F8%2F083002) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Gauge%20color%20codes%3A%20optimal%20transversal%20gates%20and%20gauge%20fixing%20in%20topological%20stabilizer%20codes&journal=New%20J.%20Phys.&volume=17&publication_year=2015&author=Bomb%C3%ADn%2CH)

-

Breuckmann, N. P. & Eberhardt, J. N. Quantum low-density parity-check codes. PRX Quantum 2, 040101 (2021).

[Article](https://doi.org/10.1103%2FPRXQuantum.2.040101) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Quantum%20low-density%20parity-check%20codes&journal=PRX%20Quantum&volume=2&publication_year=2021&author=Breuckmann%2CNP&author=Eberhardt%2CJN)

-

Pichler, H., Zhu, G., Seif, A., Zoller, P. & Hafezi, M. Measurement protocol for the entanglement spectrum of cold atoms. Phys. Rev. 6, 041033 (2016).

[Article](https://doi.org/10.1103%2FPhysRevX.6.041033) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Measurement%20protocol%20for%20the%20entanglement%20spectrum%20of%20cold%20atoms&journal=Phys.%20Rev.&volume=6&publication_year=2016&author=Pichler%2CH&author=Zhu%2CG&author=Seif%2CA&author=Zoller%2CP&author=Hafezi%2CM)

-

Brown, A. R. et al. Quantum gravity in the lab: teleportation by size and traversable wormholes. Preprint at https://arxiv.org/abs/1911.06314 (2019).

-

Koczor, B. Exponential error suppression for near-term quantum devices. Phys. Rev. 11, 031057 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXit1OrtbbL) [Article](https://doi.org/10.1103%2FPhysRevX.11.031057) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Exponential%20error%20suppression%20for%20near-term%20quantum%20devices&journal=Phys.%20Rev.&volume=11&publication_year=2021&author=Koczor%2CB)

-

Hashizume, T., Bentsen, G. S., Weber, S. & Daley, A. J. Deterministic fast scrambling with neutral atom arrays. Phys. Rev. Lett. 126, 200603 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXht1ygu7rP) [Article](https://doi.org/10.1103%2FPhysRevLett.126.200603) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Deterministic%20fast%20scrambling%20with%20neutral%20atom%20arrays&journal=Phys.%20Rev.%20Lett.&volume=126&publication_year=2021&author=Hashizume%2CT&author=Bentsen%2CGS&author=Weber%2CS&author=Daley%2CAJ)

-

Bothwell, T. et al. Resolving the gravitational redshift within a millimeter atomic sample. Nature 602, 420–424 (2022).

-

Fowler, A. G. et al. Surface code quantum communication. Phys. Rev. Lett. 104, 180503 (2010).

[Article](https://doi.org/10.1103%2FPhysRevLett.104.180503) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Surface%20code%20quantum%20communication&journal=Phys.%20Rev.%20Lett.&volume=104&publication_year=2010&author=Fowler%2CAG)

-

Zheng, X. et al. High precision differential clock comparisons with a multiplexed optical lattice clock. Nature 602, 425–430 (2022).

-

Labuhn, H. et al. Tunable two-dimensional arrays of single Rydberg atoms for realizing quantum Ising models. Nature 534, 667–670 (2016).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC28XptVajtrw%3D) [Article](https://doi.org/10.1038%2Fnature18274) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Tunable%20two-dimensional%20arrays%20of%20single%20Rydberg%20atoms%20for%20realizing%20quantum%20Ising%20models&journal=Nature&volume=534&pages=667-670&publication_year=2016&author=Labuhn%2CH)

-

Lengwenus, A., Kruse, J., Schlosser, M., Tichelmann, S. & Birkl, G. Coherent transport of atomic quantum states in a scalable shift register. Phys. Rev. Lett. 105, 170502 (2010).

[CAS](https://www.nature.com/articles/cas-redirect/1:STN:280:DC%2BC3M7htlejtg%3D%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.105.170502) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Coherent%20transport%20of%20atomic%20quantum%20states%20in%20a%20scalable%20shift%20register&journal=Phys.%20Rev.%20Lett.&volume=105&publication_year=2010&author=Lengwenus%2CA&author=Kruse%2CJ&author=Schlosser%2CM&author=Tichelmann%2CS&author=Birkl%2CG)

-

Yang, J. et al. Coherence preservation of a single neutral atom qubit transferred between magic-intensity optical traps. Phys. Rev. Lett. 117, 123201 (2016).

[Article](https://doi.org/10.1103%2FPhysRevLett.117.123201) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Coherence%20preservation%20of%20a%20single%20neutral%20atom%20qubit%20transferred%20between%20magic-intensity%20optical%20traps&journal=Phys.%20Rev.%20Lett.&volume=117&publication_year=2016&author=Yang%2CJ)

-

Wimperis, S. Broadband, narrowband, and passband composite pulses for use in advanced NMR experiments. J. Magn. Reson. A 109, 221–231 (1994).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DyaK2cXmsVertbo%3D) [Article](https://doi.org/10.1006%2Fjmra.1994.1159) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Broadband%2C%20narrowband%2C%20and%20passband%20composite%20pulses%20for%20use%20in%20advanced%20NMR%20experiments&journal=J.%20Magn.%20Reson.%20A&volume=109&pages=221-231&publication_year=1994&author=Wimperis%2CS)

-

Xia, T. et al. Randomized benchmarking of single-qubit gates in a 2D array of neutral-atom qubits. Phys. Rev. Lett. 114, 100503 (2015).

[CAS](https://www.nature.com/articles/cas-redirect/1:STN:280:DC%2BC2MnptVartQ%3D%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.114.100503) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Randomized%20benchmarking%20of%20single-qubit%20gates%20in%20a%202D%20array%20of%20neutral-atom%20qubits&journal=Phys.%20Rev.%20Lett.&volume=114&publication_year=2015&author=Xia%2CT)

-

Ozeri, R. et al. Errors in trapped-ion quantum gates due to spontaneous photon scattering. Phys. Rev. A 75, 042329 (2007).

-

Souza, A. M., Lvarez, G. A. & Suter, D. Robust dynamical decoupling. Phil. Trans. R. Soc. A 370, 4748–4769 (2012).

[Article](https://doi.org/10.1098%2Frsta.2011.0355) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Robust%20dynamical%20decoupling&journal=Phil.%20Trans.%20R.%20Soc.%20A&volume=370&pages=4748-4769&publication_year=2012&author=Souza%2CAM&author=Lvarez%2CGA&author=Suter%2CD)

-

Couvert, A., Kawalec, T., Reinaudi, G. & Guéry-Odelin, D. Optimal transport of ultracold atoms in the non-adiabatic regime. Europhys. Lett. 83, 13001 (2008).

-

Reichle, R. et al. Transport dynamics of single ions in segmented microstructured Paul trap arrays. Fortschr. Phys. 54, 666–685 (2006).

-

Carruthers, P. & Nieto, M. M. Coherent states and the forced quantum oscillator. Am. J. Phys. 33, 537 (1965).

-

Chen, Z. et al. Exponential suppression of bit or phase errors with cyclic error correction. Nature 595, 383–387 (2021).

[Article](https://doi.org/10.1038%2Fs41586-021-03721-x) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Exponential%20suppression%20of%20bit%20or%20phase%20errors%20with%20cyclic%20error%20correction&journal=Nature&volume=595&pages=383-387&publication_year=2021&author=Chen%2CZ)

-

Wilson, J. T. et al. Trapping alkaline earth rydberg atoms optical tweezer arrays. Phys. Rev. Lett. 128, 033201 (2022).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB38Xjt1yrtLo%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.128.033201) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Trapping%20alkaline%20earth%20rydberg%20atoms%20optical%20tweezer%20arrays&journal=Phys.%20Rev.%20Lett.&volume=128&publication_year=2022&author=Wilson%2CJT)

-

Madjarov, I. S. et al. High-fidelity entanglement and detection of alkaline-earth Rydberg atoms. Nat. Phys. 16, 857–861 (2020).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3cXhtVSgsLbF) [Article](https://doi.org/10.1038%2Fs41567-020-0903-z) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=High-fidelity%20entanglement%20and%20detection%20of%20alkaline-earth%20Rydberg%20atoms&journal=Nat.%20Phys.&volume=16&pages=857-861&publication_year=2020&author=Madjarov%2CIS)

-

Schine, N., Young, A. W., Eckner, W. J., Martin, M. J. & Kaufman, A. M. Long-lived Bell states in an array of optical clock qubits. Preprint at https://arxiv.org/abs/2111.14653 (2021).

-

Bluvstein, D. et al. Controlling quantum many-body dynamics in driven Rydberg atom arrays. Science 371, 1355–1359 (2021).

[MathSciNet](http://www.ams.org/mathscinet-getitem?mr=4269297) [CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXnsVyqtrs%3D) [Article](https://doi.org/10.1126%2Fscience.abg2530) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Controlling%20quantum%20many-body%20dynamics%20in%20driven%20Rydberg%20atom%20arrays&journal=Science&volume=371&pages=1355-1359&publication_year=2021&author=Bluvstein%2CD)

-

Picken, C. J., Legaie, R., McDonnell, K. & Pritchard, J. D. Entanglement of neutral-atom qubits with long ground-Rydberg coherence times. Quantum Sci. Technol. 4, 015011 (2018).

[Article](https://doi.org/10.1088%2F2058-9565%2Faaf019) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Entanglement%20of%20neutral-atom%20qubits%20with%20long%20ground-Rydberg%20coherence%20times&journal=Quantum%20Sci.%20Technol.&volume=4&publication_year=2018&author=Picken%2CCJ&author=Legaie%2CR&author=McDonnell%2CK&author=Pritchard%2CJD)

-

Steane, A. Multiple-particle interference and quantum error correction. Proc. R. Soc. Lond. A 452, 2551–2577 (1996).

[MathSciNet](http://www.ams.org/mathscinet-getitem?mr=1421749) [MATH](http://www.emis.de/MATH-item?0876.94002) [Article](https://doi.org/10.1098%2Frspa.1996.0136) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Multiple-particle%20interference%20and%20quantum%20error%20correction&journal=Proc.%20R.%20Soc.%20Lond.%20A&volume=452&pages=2551-2577&publication_year=1996&author=Steane%2CA)

-

Erhard, A. et al. Entangling logical qubits with lattice surgery. Nature 589, 220–224 (2021).

[CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BB3MXht1Okt78%3D) [Article](https://doi.org/10.1038%2Fs41586-020-03079-6) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Entangling%20logical%20qubits%20with%20lattice%20surgery&journal=Nature&volume=589&pages=220-224&publication_year=2021&author=Erhard%2CA)

-

Ho, W. W., Choi, S., Pichler, H. & Lukin, M. D. Periodic orbits, entanglement, and quantum many-body scars in constrained models: matrix product state approach. Phys. Rev. Lett. 122, 040603 (2019).

[MathSciNet](http://www.ams.org/mathscinet-getitem?mr=3925465) [CAS](https://www.nature.com/articles/cas-redirect/1:CAS:528:DC%2BC1MXnvFSqtLg%3D) [Article](https://doi.org/10.1103%2FPhysRevLett.122.040603) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Periodic%20orbits%2C%20entanglement%2C%20and%20quantum%20many-body%20scars%20in%20constrained%20models%3A%20matrix%20product%20state%20approach&journal=Phys.%20Rev.%20Lett.&volume=122&publication_year=2019&author=Ho%2CWW&author=Choi%2CS&author=Pichler%2CH&author=Lukin%2CMD)

[Download references](https://citation-needed.springer.com/v2/references/10.1038/s41586-022-04592-6?format=refman&flavour=references)

Acknowledgements

We thank I. Cong and X. Gao for their insight into QEC and for pointing us to the connections between graph states and QEC codes; M. Abobeih, H. Bernien, M. Cain, S. Choi, M. Endres, S. Evered, W. W. Ho, T. Manovitz, A. Omran, M. Serbyn and H. Zhou for reading of the manuscript and helpful discussions. We acknowledge financial support from the Center for Ultracold Atoms, the National Science Foundation, the Vannevar Bush Faculty Fellowship, the US Department of Energy (DE-SC0021013 and DOE Quantum Systems Accelerator Center, contract number 7568717), the Office of Naval Research, the Army Research Office MURI (grant number W911NF-20-1-0082) and the DARPA ONISQ programme (grant number W911NF2010021). D.B. acknowledges support from the NSF Graduate Research Fellowship Program (grant DGE1745303) and The Fannie and John Hertz Foundation. H.L. acknowledges support from the National Defense Science and Engineering Graduate (NDSEG) fellowship. G.S. acknowledges support from a fellowship from the Max Planck/Harvard Research Center for Quantum Optics. N.M. acknowledges support by the Department of Energy Computational Science Graduate Fellowship under award number DE-SC0021110. H.P. acknowledges support by the Army Research Office (grant number W911NF-21-1-0367)

Author information

Author notes

-

Harry Levine

Present address: AWS Center for Quantum Computing, Pasadena, CA, USA

Affiliations

-

Department of Physics, Harvard University, Cambridge, MA, USA

Dolev Bluvstein, Harry Levine, Giulia Semeghini, Tout T. Wang, Sepehr Ebadi, Marcin Kalinowski, Alexander Keesling, Nishad Maskara, Markus Greiner & Mikhail D. Lukin

-

QuEra Computing Inc., Boston, MA, USA

Alexander Keesling

-

Institute for Theoretical Physics, University of Innsbruck, Innsbruck, Austria

Hannes Pichler

-

Institute for Quantum Optics and Quantum Information, Austrian Academy of Sciences, Innsbruck, Austria

Hannes Pichler

-

Department of Physics and Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, MA, USA

Vladan Vuletić

Contributions

D.B., H.L., G.S., T.T.W., S.E., M.K. and A.K. contributed to building the experimental setup, performed the measurements and analysed the data. M.K., N.M. and H.P. performed theoretical analysis. All work was supervised by M.G., V.V. and M.D.L. All authors discussed the results and contributed to the manuscript.

Corresponding author

Correspondence to [Mikhail D. Lukin](mailto:lukin at physics.harvard.edu).

Ethics declarations

Competing interests

M.G., V.V. and M.D.L. are co-founders and shareholders of QuEra Computing. A.K. is an executive at and shareholder of QuEra Computing. All other authors declare no competing interests.

Peer review

Peer review information

Nature thanks Hannah Williams and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. [Peer reviewer reports](https://www.nature.com/articles/s41586-022-04592-6#MOESM1) are available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

[Extended Data Fig. 1 CZ gate echo, atomic level structure, and typical pulse sequence.](https://www.nature.com/articles/s41586-022-04592-6/figures/5)

a, The two-qubit gates we apply, in addition to applying a controlled-Z operation between the two qubits, also induce a single-qubit phase Z(ζ) to both qubits, composed of the intrinsic phase of the CZ gate[5](https://www.nature.com/articles/s41586-022-04592-6#ref-CR5) and additional spurious phases from the 420-nm Rydberg laser and pulsing the traps off. Since we apply all gates in parallel by global pulses of the Rydberg laser, if a qubit is not adjacent to another qubit, it does not perform a CZ gate but still acquires the same Z(ζ) (identical to being adjacent to another qubit in state |0⟩, which is dark to the Rydberg laser). As diagrammed, we cancel the additional, undesired Z(ζ) by applying a π pulse between pairs of CZ gates. This echo procedure removes any need to calibrate the intrinsic phase from the CZ gate, and renders us insensitive to any spurious changes in Z(ζ) slower than ~200 μs. The additional Y(π) propagates in a known way through the CZ gates and multiplies certain stabilizers by a −1 sign, which simply redefines the sign of stabilizers and logical qubits. b, Level diagram showing key 87Rb atomic levels used. Our Rydberg excitation scheme from |1⟩ to |r⟩ is composed of a two-photon transition driven by a 420-nm laser and a 1013-nm laser (see ref. [25](https://www.nature.com/articles/s41586-022-04592-6#ref-CR25) for description of laser system). A DC magnetic field of B = 8.5 G is applied throughout this work. c, A typical pulse sequence for running a quantum circuit.

[Extended Data Fig. 2 Movement characterization and multiple drop-recaptures.](https://www.nature.com/articles/s41586-022-04592-6/figures/6)

a, Atom retention as a function of average separation speed 2D/T (as is plotted in Fig. [1d](https://www.nature.com/articles/s41586-022-04592-6#Fig1) of the main text for separating Bell pairs), with subtracted background loss of 0.7%. The inset in Fig. [1d](https://www.nature.com/articles/s41586-022-04592-6#Fig1) of the main text is normalized by (Atom retention)[2](https://www.nature.com/articles/s41586-022-04592-6#ref-CR2) (without subtracting background loss). Dark curve is calculated using experimental parameters and Eq. [2](https://www.nature.com/articles/s41586-022-04592-6#Equ2), matched to the experimental data by setting Nmax = 26 and averaging over a range of ω0/2π of ±15% around an average ω0/2π = 40 kHz. b, Atom retention as a function of inverse trap frequency (2π/ω0) after the four moves of the surface code circuit. For calculating the atom loss here we set Nmax = 33 and average the trap frequencies over a range of ± 15%. We note that these quantitative estimates are sensitive to ω0 which we roughly estimate. c, Atom loss as a function of drop time and number of drop loops, with 100 μs wait between each drop. When running quantum circuits we use 500-ns drops for each CZ gate (to avoid anti-trapping of the Rydberg state and light shifts of the transition), for which we observe here that hundreds of drops can be made (corresponding to hundreds of possible CZ gates per atom) before atom loss becomes significant. d, By rescaling the x-axis of the data to tdropN−−√

, we find the data of the various tdrop collapse onto a universal curve, suggesting a diffusion model for explaining the atom loss after a certain number of drops. By modeling such a diffusion process analytically we obtain the black curve with a temperature of 10 μK and a trapping radius of 1 μm.

[Extended Data Fig. 3 Robust single-qubit control and qubit coherence.](https://www.nature.com/articles/s41586-022-04592-6/figures/7)

a, Robust BB1 single-qubit rotation in comparison to a normal single-qubit rotation, as a function of pulse area error. An arbitrary BB1(θ, ϕ) rotation on the Bloch sphere of angle θ about axis ϕ is realized with a sequence of four pulses: (π)φ + ϕ(2π)3φ + ϕ(π)φ + ϕ(θ)ϕ, where φ = cos−1 (−θ/4π)[60](https://www.nature.com/articles/s41586-022-04592-6#ref-CR60). Pulse fidelity is measured here for a π pulse, defined such that the fidelity is the probability of successful transfer from |0⟩ → |1⟩, including SPAM correction. b, Preserving hyperfine qubit coherence using dynamical decoupling (XY16 with 128 total π pulses). Qubit coherence is observed on a timescale of seconds, with a fitted coherence time T2 = 1.49(8)s. Data is measured with either a + π/2 or −π/2 pulse at the end of the sequence, and these curves are then subtracted to yield the coherence y-axis. c, Hyperfine qubit T1, measured by the difference of final F = 2 populations between measurements starting in |F = 2, mF = 0⟩ and |F = 1, mF = 0⟩. Atom loss without cooling is separately measured (predominantly arising from vacuum loss) and normalized to also measure the intrinsic spin relaxation time T′1

in the absence of atom loss. All data here is measured in 830-nm traps.

[Extended Data Fig. 4 Effect of axial trap oscillations on echo fidelity of 420-nm Rydberg pulse.](https://www.nature.com/articles/s41586-022-04592-6/figures/8)

a, Noise correlation measurement of the 420-nm Rydberg laser pulse intensity. In the blue-detuned configuration used in this figure only, the 420-nm laser induces an 8 MHz differential light shift on the hyperfine qubit, and consequently a phase accumulation of 32π during a 2-μs pulse (our CZ gates are 400-ns total). Small fluctuations of the 420-nm laser intensity lead to large fluctuations in phase accumulation of the hyperfine qubit, and thus cause significant dephasing. The echo sequence diagrammed here probes the correlation of the accumulated phase between two 420-nm pulses separated by a variable time τ, and thus informs how far-separated in time the 420-nm pulses can be while still properly echoing out fluctuations in the 420-nm intensity. b, Hyperfine coherence (a proxy for echo fidelity) versus gap time τ between the two 420-nm pulses. The echo fidelity decreases initially due to a decorrelation of the 420-nm intensity, but then increases again, showing that the correlation of the 420-nm intensity is non-monotonic. The decaying oscillations are fit to a functional form of y = y0 + Acos2(πfτ)exp[−(τ/T)2]. c, The fitted oscillation frequency f of the correlation / decorrelation of the noise follows a square-root relationship with the trap power, and is consistent with the expected axial trap oscillation frequency. These observations indicate that a significant portion of the correlation / decorrelation of the 420-nm noise arises from the several-μm axial oscillations of the atom in the trap. For this measurement, we intentionally displace the 420-nm beam by several μm in order to place the atom on a slope of the beam, increasing our sensitivity to this phenomenon. For the other experiments in our work, we minimize sensitivity to these effects by operating in the center of a larger (35-micron-waist) 420-nm beam and operating red-detuned of the intermediate-state transition.

[Extended Data Fig. 5 Movement schematics.](https://www.nature.com/articles/s41586-022-04592-6/figures/9)

Schematics showing the gate-by-gate creation of a the 1D cluster state, b the Steane code, c the surface code, and d the toric code (see also [Supplementary Video 1](https://www.nature.com/articles/s41586-022-04592-6#MOESM2)), in a side-by-side comparison. These various graph states are all generated in the same way, and encoding a desired circuit is a matter of positioning the atoms in different initial positions and applying an appropriate AOD waveform. To realize a desired circuit, atom layouts and trajectories are optimized heuristically in the way described in the Methods text. Panel c also shows the definition of surface code stabilizers as ordered in the main text.

[Extended Data Fig. 6 Error simulations and tabulated single-qubit and two-qubit error estimates.](https://www.nature.com/articles/s41586-022-04592-6/figures/10)

We compare our measured graph state fidelities to those from a stochastic Monte Carlo simulation of our system for a, the surface code and b, the toric code. We find that the simulated stabilizers agree well with the experimental data for this empirical depolarizing noise model. In addition, for the surface code (toric code) in the experiment we find 35% (20%) of measurements detect no stabilizer errors, compared to 40% (26%) in the simulation. We assume two-qubit errors are described by rates of 0.2% Y error, 0.2% X error, 0.5% Z error, and 0.5% loss per qubit per parallel layer (4 layers for surface code, 5 layers for toric code), corresponding to a 97.2% CZ-gate fidelity. We also add ambient, single-qubit errors at a rate of 0.1% Y error, 0.1% X error, 0.4% Z error, and 0.2% loss per qubit per parallel layer, as well as an initial 1% loss before the circuit begins (empirically factoring in SPAM errors). c, Tabulation of single-qubit (SQ) and two-qubit (TQ) gate errors that are measured, estimated, and extrapolated. Simulated TQ fidelities include the 0.6% scattering error from the 420-nm echo pulse. The estimated TQ fidelities are given for the experiments of the surface code and toric code, but is an underestimate of the TQ fidelities for the cluster state and Steane code measurements where we increase the 1013-nm intensity by 2× and reduce the 420-nm intensity by 2×, increasing gate fidelity. The Bell state estimate of CZ gate fidelity is similarly done with 2× higher 1013 intensity, but includes the 420-nm echo pulse, and consequently yields a similar gate fidelity as the surface and toric code estimates.

[Extended Data Fig. 7 Properties of encoded logical states.](https://www.nature.com/articles/s41586-022-04592-6/figures/11)

a, Summary of logical error probabilities for the various error correcting graphs made in this work (all in logical state |+⟩L), for raw measurements as well as implementing error correction and error detection in postprocessing. Error correction for the Steane code is implemented with the Steane code decoder[36](https://www.nature.com/articles/s41586-022-04592-6#ref-CR36),[73](https://www.nature.com/articles/s41586-022-04592-6#ref-CR73) and is implemented with the minimum-weight-perfect-matching algorithm for the surface and toric codes[38](https://www.nature.com/articles/s41586-022-04592-6#ref-CR38). For the even-distance toric code, when correction is ambiguous we do not flip the logical qubit, and accordingly the distance d = 2 logical qubit does not change under the correction procedure. We remark that the observed fidelities are comparable to similar demonstrations in state-of-the-art experiments with other platforms[8](https://www.nature.com/articles/s41586-022-04592-6#ref-CR8),[74](https://www.nature.com/articles/s41586-022-04592-6#ref-CR74). These will need to be improved to surpass the threshold for practical error correction[38](https://www.nature.com/articles/s41586-022-04592-6#ref-CR38) (see [Methods](https://www.nature.com/articles/s41586-022-04592-6#Sec7) text). b, Lifetime of the logical |+⟩L state on the surface code, with correction and detection performed in postprocessing as in a. After state preparation, the |+⟩L state is held for a variable time before projective measurement, with two π pulses applied for dynamical decoupling (lifetime can be extended significantly further by applying e.g. 128 π pulses as done in Extended Data Fig. [3b](https://www.nature.com/articles/s41586-022-04592-6#Fig7)). Some experimental parameters are slightly different here compared to those in a, hence the higher error rates here at the time 0 point. c, Logical π/2 rotation on the Steane code to prepare logical qubit state |0L⟩. The Steane code, surface code, and toric code all have transversal single-qubit Clifford operations on the logical qubit[8](https://www.nature.com/articles/s41586-022-04592-6#ref-CR8),[36](https://www.nature.com/articles/s41586-022-04592-6#ref-CR36) (including in-software rotations of the lattice), which is a high-fidelity operation in our system since the transveral rotations are implemented in parallel with our global Raman laser and the physical single-qubit fidelities are high. We show a logical π/2 rotation here for the Steane code as an example but emphasize that we can readily realize the various basis states for all of these codes.

[Extended Data Fig. 8 Benchmarking the interferometry measurement.](https://www.nature.com/articles/s41586-022-04592-6/figures/12)

a, To benchmark our gate-based interferometry technique, we prepare variable single-particle pure states (by applying a variable-length resonant Raman pulse) and then reconfigure the system and apply the interferometry circuit on twin pairs. The interferometry circuit converts the anti-symmetric singlet state |Ψ−⟩ to the computational basis state |00⟩, while converting the symmetric triplet states to other computational states. We plot the resulting twin pair output states in the left panel. We rarely observe the |00⟩ state (1.95(2)% of measurements), with a measurement fidelity independent of the initial state. This low probability P00 of observing |00⟩ corresponds to a high extracted single-particle purity of 2P00 − 1 = 0.961(3) (right panel). We find this measurement to be a useful benchmark, as interferometry miscalibrations can result in significant state-dependence of the observed purity that would then compromise the validity of the many-body entanglement entropy measurement. b, Benchmarking the entanglement entropy measurement with Bell state arrays. (Top) Microstate populations during two-particle oscillations between |11⟩ and 12√(|1r⟩+|r1⟩)

under a Rydberg pulse of variable duration. Faint lines show measurement results in the {|1⟩, |r⟩} basis, and dark lines show results in the {|0⟩, |1⟩} basis after the coherent mapping process. (Bottom) Measured local and global purities by analyzing the number parity of observed |00⟩ twin pairs in each measurement. For this two-particle data we use a gap of 230 ns in the coherent mapping sequence as opposed to the 150-ns gap used in the many-body data.

[Extended Data Fig. 9 Raw many-body data and numerical modeling of errors.](https://www.nature.com/articles/s41586-022-04592-6/figures/13)

a, Raw measured Renyi entropy without subtracting the extensive classical entropy, as a function of subsystem size for quenches from |rgrgrgrg⟩ and |gggggggg⟩. The Renyi entropy of the 4-atom subsystem is the same underlying data as the half-chain entanglement entropy plotted in Fig. [4d](https://www.nature.com/articles/s41586-022-04592-6#Fig4) of the main text. In the main text, we subtract the data by a fixed offset given by the classical entropy-per-particle, corresponding to the time = 0 offset for each subsystem size. The extensive, classical entropy offset is slightly larger for the |rgrgrgrg⟩ quench due to non-unity fidelities both of preparing |r⟩ and mapping |r⟩ → |1⟩. b, Raw global purity after the |gggggggg⟩ quench. The global purity is a sensitive proxy for the fidelity of our entire process. We find this 16-body observable, composed of three-level systems, remains > 100× the purity expected for a fully mixed state of 8 qubits (1/28) (see inset). For comparison of scale we also plot single-particle purity to the 8th power, to indicate what the global purity would be if the measurement results on each twin were uncorrelated. c, Global purity for the 8-atom quench calculated through numerical modeling of the three-level system {|0⟩, |1⟩ ≡ |g⟩, |r⟩} with a variety of simulated error sources. We model the experimentally measured purity by calculating the expectation value of the SWAP operator in the {|0⟩, |1⟩} basis between two independent chains, taking into account that residual population in |r⟩ results in atom loss and measurement associated with the +1 eigenvalue of the SWAP operator (as the twin state |00⟩ can no longer be detected). The top curve includes only errors from population left in |r⟩ following the coherent mapping step (see [methods](https://www.nature.com/articles/s41586-022-04592-6#Sec7) text). The second-from-top curve includes single-site dephasing (T∗2

) during the Rydberg dynamics and the coherent mapping gap, modeled by a random on-site detuning which is Gaussian-distributed with zero mean and standard deviation of 100 kHz. The third and fourth curves include multiplication by the experimentally observed raw global purity at quench time t = 0, and then further multiplying empirically by an exponential decay exp[−16 × t/(70 μs)] as a simple model for scattering and decay errors with an experimentally estimated rate of roughly 70 μs for each of the 16 atoms between the two chains.

[Extended Data Fig. 10 Local observables and entanglement entropy for quantum many-body scars.](https://www.nature.com/articles/s41586-022-04592-6/figures/14)

a, Experimentally measured single-site entropy for each site in the 8-atom chain when quenching from the scarred |ℤ2⟩ state, including the classical entropy subtraction. Solid curves plot exact, ideal (imperfection-free) numerics of HRyd (Eq. [3](https://www.nature.com/articles/s41586-022-04592-6#Equ3)); excellent agreement between data and numerics is found for every atom in the chain. b, (Top) Same data as Fig. [4f](https://www.nature.com/articles/s41586-022-04592-6#Fig4) of the main text, showing single-site entropy of the middle two atoms in the chain, for two different initial states[75](https://www.nature.com/articles/s41586-022-04592-6#ref-CR75). (Bottom) Measurements of the many-body state in the Z-basis with the interferometry circuit turned off. Characteristic of the scars from the |ℤ2⟩ = |rgrgrgrg⟩ state, the Rydberg excitation probability on the sublattices exhibits periodic oscillations[71](https://www.nature.com/articles/s41586-022-04592-6#ref-CR71). In the bottom row, the dark data points are measured in the {|1⟩, |r⟩} basis, and the faint data points are measured in the {|0⟩, |1⟩} basis after the coherent mapping sequence. Measurements in both bases agree well with exact numerics (solid lines), which we emphasize has no free fit parameters and does not account for any experimental imperfections, such as detection infidelity. Moreover, the data indicate the high fidelity of preparation into the |ℤ2⟩ state by use of local Rydberg π pulses. In plotting, we delay the theory curves and the {|1⟩, |r⟩} basis measurement by 10 ns to account for the fact that the Raman π pulse we apply cuts off the final 10 ns of the Rydberg evolution, when measuring in the {|0⟩, |1⟩} basis. c, Numerical simulations of the single-site Renyi entropy on two adjacent sites in the idealized ‘PXP’ model of perfect nearest-neighbor blockade[14](https://www.nature.com/articles/s41586-022-04592-6#ref-CR14). The system size is 24 atoms with periodic boundary conditions, showing the same out-of-phase oscillations in the entanglement entropy of the two sublattices. d, Diagram of the constrained Hilbert space of the system[14](https://www.nature.com/articles/s41586-022-04592-6#ref-CR14). The early-time, out-of-phase entropy oscillations[75](https://www.nature.com/articles/s41586-022-04592-6#ref-CR75) of the scars can be understood in this constrained Hilbert space picture, where the scar dynamics are known to take the state from the left end (|rgrgrgrg⟩) to the right end (|grgrgrgr⟩) (dark circles represent |r⟩ and white circles represent |g⟩)[14](https://www.nature.com/articles/s41586-022-04592-6#ref-CR14). Near these crystalline ends of this constrained Hilbert space, the Rydberg atoms can fluctuate (high entropy), but the ground state atoms are pinned (low entropy). Our analysis shows that entanglement between atoms on the same sublattice contributes to the eventual degradation of the revival fidelity of the |ℤ2⟩ state.

Supplementary information

[Peer Review File](https://static-content.springer.com/esm/art%3A10.1038%2Fs41586-022-04592-6/MediaObjects/41586_2022_4592_MOESM1_ESM.pdf)

[41586_2022_4592_MOESM2_ESM.mp4](https://static-content.springer.com/esm/art%3A10.1038%2Fs41586-022-04592-6/MediaObjects/41586_2022_4592_MOESM2_ESM.mp4)

Supplementary Video 1 Generation of the toric code using coherent transport of entangled atom arrays. The video depicts the formation of the toric code state on a torus, realized by dynamically changing the quantum processor connectivity throughout the computation time of 1 ms. Ancillas are transported across the array and entangled with stationary data qubits using a two-qubit CZ gate (red ovals), and are then moved to a separate zone for local rotation before readout of all qubits.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: text/html
Size: 426654 bytes
Desc: not available
URL: <https://lists.cpunks.org/pipermail/cypherpunks/attachments/20220422/862a26ac/attachment.txt>