Erasure conversion in a high-fidelity Rydberg quantum simulator

Minimizing and understanding errors is critical for quantum science, both in noisy intermediate scale quantum (NISQ) devices1 and for the quest towards fault-tolerant quantum computation2,3. Rydberg arrays have emerged as a prominent platform in this context4 with impressive system sizes5,6 and proposals suggesting how error-correction thresholds could be significantly improved by detecting leakage errors with single-atom resolution7,8, a form of erasure error conversion9–12. However, two-qubit entanglement fidelities in Rydberg atom arrays13,14 have lagged behind competitors15,16 and this type of erasure conversion is yet to be realized for matter-based qubits in general. Here we demonstrate both erasure conversion and high-fidelity Bell state generation using a Rydberg quantum simulator5,6,17,18. When excising data with erasure errors observed via fast imaging of alkaline-earth atoms19–22, we achieve a Bell state fidelity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ge 0.997{1}_{-13}^{+10}$$\end{document}≥0.9971−13+10, which improves to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ge 0.998{5}_{-12}^{+7}$$\end{document}≥0.9985−12+7 when correcting for remaining state-preparation errors. We further apply erasure conversion in a quantum simulation experiment for quasi-adiabatic preparation of long-range order across a quantum phase transition, and reveal the otherwise hidden impact of these errors on the simulation outcome. Our work demonstrates the capability for Rydberg-based entanglement to reach fidelities in the 0.999 regime, with higher fidelities a question of technical improvements, and shows how erasure conversion can be utilized in NISQ devices. These techniques could be translated directly to quantum-error-correction codes with the addition of long-lived qubits7,22–24.

Minimizing and understanding errors is critical for quantum science, both in noisy intermediate scale quantum (NISQ) devices 1 and for the quest towards fault-tolerant quantum computation 2,3 .Rydberg arrays have emerged as a prominent platform in this context 4 with impressive system sizes 5,6 and proposals suggesting how error-correction thresholds could be significantly improved by detecting leakage errors with single-atom resolution 7,8 , a form of erasure error conversion [9][10][11][12] .However, two-qubit entanglement fidelities in Rydberg atom arrays 13,14 have lagged behind competitors 15,16 and this type of erasure conversion is yet to be realized for matter-based qubits in general.Here we demonstrate both erasure conversion and high-fidelity Bell state generation using a Rydberg quantum simulator 5,6,17,18 .When excising data with erasure errors observed via fast imaging of alkalineearth atoms [19][20][21][22] , we achieve a Bell state fidelity of ≥0.9971 +10 −13 , which improves to ≥0.9985 +7 −12 when correcting for remaining state preparation errors.We further apply erasure conversion in a quantum simulation experiment for quasi-adiabatic preparation of long-range order across a quantum phase transition, and unveil the otherwise hidden impact of these errors on the simulation outcome.Our work demonstrates the capability for Rydberg-based entanglement to reach fidelities in the ∼0.999 regime, with higher fidelities a question of technical improvements, and shows how erasure conversion can be utilized in NISQ devices.The shown techniques could be translated directly to quantum errorcorrection codes with the addition of long-lived qubits 7,[22][23][24] .
We begin by detailing our erasure conversion scheme and how it is employed in conjunction to Bell state generation, resulting in fidelities competitive with other state-of-the-art platforms 15,16,25,26 .Our experimental apparatus has been described in detail before 13 , and is based on trapping individual strontium atoms in arrays of optical tweezers 19,20 (Methods).Strontium features a rich energy structure, allowing us to utilize certain energy levels as a qubit subspace to perform entangling operations and separate levels for detection of leakage errors (Fig. 1a).
To controllably generate entanglement between atoms, we employ Rydberg interactions [27][28][29] .When two atoms in close proximity are simultaneously excited to highlying electronic energy levels, called Rydberg states, they experience a distance-dependent van der Waals interaction V = C 6 /r 6 , where r is the interatomic spacing, and C 6 is an interaction coefficient.If the Rabi frequency, Ω, which couples the ground, |g⟩, and Rydberg, |r⟩, states is much smaller than the interaction shift, Ω/V ≪ 1, the two atoms cannot be simultaneously excited to the Rydberg state (Fig. 1b, inset), a phenomena known as Rydberg blockade.In this regime, the laser drives a unitary operation, Û (t), that naturally results in the two atoms forming a Bell state, |Ψ + ⟩ = 1 √ 2 (|gr⟩ + |rg⟩), between the ground and Rydberg states (Fig. 1b).
This Bell state generation has several major practical limitations.Of particular interest here are leakage errors to the absolute ground state, 1 S 0 , which are converted to erasure errors in our work as described below (and in Ext.Data Fig. 1).The first error of this type is imperfect preparation of atoms in |g⟩ prior to applying Û (t).The second arises from decay out of the Rydberg state along multiple channels.We distinguish decay into 'bright' states, which we can image, and 'dark' states, which are undetected (Ext.Data Fig. 2).The former primarily refers to low-lying energy states which are repumped to 1 S 0 as part of the imaging process or decay to 1 S 0 via intermediate states, while the latter mainly consists of nearby Rydberg states accessed via blackbody radiation.
Here we employ a scheme, theoretically proposed 7 but not yet demonstrated, which allows us to detect the location of such leakage errors (Fig. 1b), converting them into so-called erasure errors, i.e., errors with a known location 9 .To this end, we demonstrate fast, 24 µs imaging of atoms in 1 S 0 (Ext.Data.Fig. 1) with single-site resolution and 0.980 +1 −1 fidelity.Such fast imaging had previously been performed for a few, freely-propagating, alkali atoms 30 , but not for many trapped atoms in tweezer arrays or alkaline-earth atoms (Methods).
Our general procedure is shown in Fig. 1b (further detailed in Ext.Data Fig. 3).We first rearrange 31,32 atoms into pairs, coherently transfer them to |g⟩, and then perform the entangling Û operation.Immediately after, we auto-ionize the atoms to project the populations of the resultant state.  1 | Erasure conversion for high-fidelity entanglement.a, Level structure used in this work.We distinguish two subspaces: a qubit subspace in which the atoms interact via their Rydberg states, and a measurement subspace used to detect leakage errors from the qubit subspace with single-site resolution, realizing erasure conversion.b, Sketch of the erasure conversion scheme, as applied to Bell pair generation.After arranging atoms into pairs (top) we prepare them in |g⟩, and entangle them via the Rydberg blockade mechanism (right), denoted by a unitary operation Û (t).Immediately afterwards, we auto-ionize atoms in |r⟩, effectively projecting the populations of the Bell states, and follow with a fast erasure conversion image to detect leakage out of the qubit subspace during the preparation or evolution periods.This is followed by the final detection of atoms in |g⟩, yielding two separate, independent images.We can discard data from pairs where atoms are detected in the erasure-error image, termed erasure-excision in the following.Atom fluorescence images are single-shot, with post-processing applied to improve detection fidelity 30 (Methods).c, Lower bounds for Bell state fidelities with (blue) and without (red) the erasure-excision, and using incoherent pumping to reduce preparation errors instead of erasure-excision (green, Methods).We present the results for the raw data, corrected for measurement errors, and corrected for state preparation and measurement errors (SPAM).All data is averaged over eight pairs of atoms which are excited in parallel.Error bars represent a 68% confidence interval (Ext.Data Fig. 5, Methods).
We then perform the fast erasure image; any atoms which are detected are concluded to be the result of some leakage error process.Importantly, the erasure image does not affect atoms remaining in |g⟩, and is extremely short compared to its lifetime, resulting in a survival probability in |g⟩ of 0.9999954 +12 −12 (Ext.Data Fig. 1, Methods).Hence, the erasure image does not perturb the subsequent final readout.Thus, we obtain two separate images characterizing a single experimental repetition, the final image showing the ostensible result of Û , and the erasure image revealing leakage errors with single-site resolution.
We note that this work is not a form of mid-circuit detection as no superposition states of |g⟩ and |r⟩ exist at the time of the erasure image.Instead, our approach is a noise mitigation strategy via erasure-excision, where experimental realizations are discarded if erasures are detected.In contrast to other leakage mitigation schemes previously demonstrated in matter-based qubit platforms [33][34][35] , we directly spatially resolve leakage errors in a way which is decoupled from the performed experiment, is not post-selected on the final qubit readout, and does not require any extra qubits to execute.
However, the coherence between |g⟩ and |r⟩ can in principle be preserved during erasure detection for future applications; in particular, we see no significant difference in Bell state lifetime with and without the imaging light for erasure detection on (Ext.Data Fig. 4, Methods).We also expect long-lived nuclear qubits encoded in |g⟩ to be unperturbed by our implementation of erasure conversion 7,[22][23][24] .

Bell state generation results
With a procedure for performing erasure conversion in hand, we now describe its impact on Bell state genera-tion.Experimentally, we only obtain a lower-bound for the Bell state generation fidelity 13 (Methods, Ext.Data Fig. 5); the difference of this lower bound to the true fidelity is discussed further below.
We first coherently transfer atoms to |g⟩ as described before, and then consider three scenarios (Fig. 1c, Ext.Data Table 1).In the first, as a baseline we perform the entangling unitary Û without considering any erasure detection results (red bars).In the second, we excise data from any pairs of atoms with an observed erasure error (blue bars).Finally, we compare against another strategy for mitigating preparation errors through incoherent repumping 13 , but without erasure detection (green bars).Notably, the raw value for the Bell state lowerbound with erasure-excision is ≥0.9962 +10 −13 , significantly higher than with the other methods.This difference mainly comes from erasure excision of preparation errors and, to a much lower degree, Rydberg decay.These contribute at the level of ∼5 × 10 −2 and 1.2 +3 −3 × 10 −4 , respectively (Methods).
Correcting for final measurement errors, we find a lower bound of ≥0.9971 +10 −13 , which quantifies our ability to generate Bell pairs conditioned on finding no erasure events.To quantify the quality of the Rydberg entangling operation Û (t) itself, we further correct for remaining preparation errors that are not detected in the erasure image (Methods), and find a state preparation and measurement (SPAM) corrected lower bound of ≥0.9985 +7 −12 .To our knowledge, these bare, measurement-corrected, and SPAM-corrected values are respectively the highest two-qubit entanglement fidelities measured for neutral atoms to date, independent of the means of entanglement generation.While Bell state generation as demon- Fig. 2 | Predicting infidelities at the 10 −3 level.a, SPAM-corrected Bell state fidelity as a function of the ratio of interaction energy and Rabi frequency, V /Ω.Error-free simulations (blue line) show fidelities continually increase with increasing V /Ω, in agreement with results from perturbation theory (dashed line).For large enough interaction strength (V /Ω > 50), other error sources become dominant, and we employ a noisy open system dynamics simulation from which we obtain an estimate of the true fidelity (light grey fill) and for the lower bounding procedure used in experiment (dark grey fill).We find good agreement between simulation and experimental results with erasure-excision (red markers).b, Predicted Bell state fidelity for V /Ω = 140 from simulations turning on a single noise term at a time.Dominant limitations come from laser frequency and intensity noise, as well as decay of the Rydberg state into dark states.We also show the results when taking into account all errors (Methods), for both the true fidelity and the lower bound estimation (right).The lower-bound significantly underestimates the true fidelity.Shaded areas in a and error bars in b represent the standard deviation of the mean over 5000 trajectories.strated here is not a computational two-qubit quantum gate -which requires additional operations -our results are indicative of the fidelities achievable in Rydberg based gate operations.
Error modelling Importantly, we understand remaining errors in the entangling operation as well the nature of detected erasure errors from a detailed ab-initio error model simulation for SPAM-corrected fidelities (Methods, Fig. 2).We identify limited interaction strength as a dominant effect that restricted SPAM-corrected entanglement fidelities in our previous work 13 (Fig. 2a); in particular, one major difference here is that we operate at smaller distance and hence larger V /Ω.In line with experimental data (red markers), fidelities at large distances are limited to F Bell ≤ 1 − 5 8 (Ω/V ) 2 obtained from perturbation theory (black dashed line, Methods).
For strong enough interaction, V /Ω > 50, corresponding to distances r < 3 µm, other error sources become limiting.In this short-distance regime, the experimental SPAM-corrected fidelity lower-bound is in good agreement with the error model prediction of ≥0.99881 +3 −3 (dark grey fill).
Our error model results show that the lower bound procedure significantly underestimates the true fidelity (light grey fill), found to be 0.99931 +6 −6 .This effect arises because the lower bound essentially evaluates the fidelity of Û by a measurement after performing Û twice (Methods), meaning particular errors can be exaggerated.Given the good match of the error model and experimental fidelity lower bounds, we expect this effect to be present in experiment as well, and to underestimate the true SPAM-corrected fidelity by about 5 × 10 −4 .
The remaining infidelity is a combination of multiple errors.In Fig. 2b, we report an error budget for the most relevant noise source contributions to the Bell state infidelity (Methods) at the experimentally chosen V /Ω = 140.Frequency and intensity laser noise are dominant limitations, but could be alleviated by improving the stability of laser power, and reducing its linewidth, for instance via cavity filtering 36 .Eliminating laser noise completely would lead to fidelities of ∼0.9997 in our model.The other major limit is Rydberg state decay into dark states, which cannot be converted into an erasure detection with our scheme.This decay is mostly blackbody induced 7,37 , and thus could be greatly reduced by working in a cryogenic environment 38 , leaving mostly spontaneous decay that is bright to our erasure detection.Accounting for these improvements, it is realistic that Rydberg-based Bell state generation in optical tweezers arrays could reach >0.9999 fidelity in the coming years.

Quantum simulation with erasure conversion
Having demonstrated the benefits of erasure-excision for the case of improving two-qubit entanglement fidelities, we now show it can be similarly applied to the case of many-body quantum simulation, demonstrating the utility of erasure detection for NISQ applications.As part of this investigation, we also distinguish erasure errors from preparation and Rydberg spontaneous decay, the latter of which becomes more visible in a many-body setting and for longer evolution times.
As a prototypical example, we explore a quasiadiabatic sweep into a Z 2 -ordered phase (Fig. 3a) through the use of a varying global detuning 39 (Fig. 3b).In this ordered phase, ground and Rydberg states form an antiferromagnetic (AFM) pattern, with long-range order appearing at a quantum phase transition.Unlike previous examples 17,40 , we operate in the effectively attractive interacting regime of the Rydberg blockaded space 39 , which features a true two-fold degenerate ground state for systems with an even number of atoms, even for open boundary conditions (Methods), and without explicitly modifying the boundary 40 .The ground state in the deeply ordered limit consists of two oppositely ordered AFM states, |grgr...gr⟩ and |rgrg...rg⟩.Staying adiabatic during ground state preparation requires evolution over microseconds, orders of magnitude longer than the two-qubit entanglement operation shown before, which magnifies the effect of Rydberg decay.In order to differentiate between leakage out of the qubit manifold due to either preparation errors or Rydberg decay, we perform two erasure images, one before the adiabatic sweep which captures preparation errors, and one after (Fig. 3c).The second image allows us to measure Rydberg decay into the detection subspace throughout the sweep.For a system size of N = 26 atoms (Fig. 3d), we see the number of detected preparation erasures (orange markers) stays constant over the course of a 3 µs sweep; conversely, the number of detected decay erasures (green markers) grows over time, in good agreement with the measured Rydberg lifetime and erasure image infidelities (green solid line, Methods).
With the ability to distinguish these effects, we plot the total probability to form either of the AFM states, P AFM = P (|grgr...gr⟩) + P (|rgrg...rg⟩) (Fig. 3e).At the conclusion of the sweep, we find P AFM = 0.33 +2 −2 without any erasure-excision (pink markers).By excising instances with preparation erasures, this fidelity is improved to 0.44 +2 −2 (orange markers), and is then further improved to 0.49 +2 −2 by additionally excising Rydberg decay erasures.The sharpness of the signal, exemplified by the derivative of P AFM with respect to the detuning, is similarly improved near the phase boundary (Fig. 3e, inset).We also observe that the gain in P AFM from erasure-excision increases with system size (Ext.Data Fig. 6).
We further explore how errors affect quantities reflecting higher-order statistics.To this end, we explore the probability distribution to find magnetic order of different magnitude by studying the AFM magnetization operator, defined as where ẐS = j∈S Ẑj is the total magnetization operator in sub-lattice S=A (odd sites) or S=B (even sites) respectively, N S is the number of atoms in each sublattice, and Ẑj = |r⟩⟨r|−|g⟩⟨g| is the local magnetization at site j.We plot the probability to find a specific eigenvalue, M , of M as a function of detuning (Fig. 3f).While the values of M are initially tightly grouped around M = 0 in the disordered phase, as the sweep progresses the probability distribution bifurcates, forming two separate dominant peaks in the Z 2 phase, consistent with aforementioned two-fold spontaneous symmetry breaking across the quantum phase transition.We find that erasure-excision improves the sharpness of the distribution in both the disordered and Z 2 phases (Fig. 3g).Near the phase transition, the distribution is close-toflat, consistent with order appearing at all length scales.These results demonstrate improvements in fidelity for preparation of long-range-ordered ground states with erasure-excision in quantum simulation experiments, a first proof-of-principle for utilizing erasure conversion in NISQ-type applications.
Learning from erasure errors Finally, we turn to studying a new tool enabled by our implementation of erasure conversion: exploring the Post-processed 30 single-shot atom fluorescence images (Methods).After arranging the array (top panel) we use the fast erasure images to learn how detected errors (middle panels) affect detection outcomes (lower panel).b, Conditional probability to detect no atom in the final image as a function of sweep progress and distance from a hypothetical erasure event.In the case of no erasures (left), the resulting profile is uniform.However, when conditioning on detecting a preparation erasure (middle), the error breaks the Z2 symmetry by establishing a single AFM order.In the case of conditioning on decay errors (right) the situation is more complex.c, AFM magnetization (see Eq. 1) as a function of time.Preparation erasures (orange markers) lead to a growth of a single AFM order with Rydberg excitations predominantly on sub-lattice A (defined as sites an odd distance from the erasure position).Decay erasures (green markers) follow a similar trend at early times by acting as effective preparation errors, but past the critical point (vertical dashed line) their behavior reverses: a decay erasure spontaneously breaks the two-fold symmetry, where the neighboring sites must have been in the ground state due to the Z2 structure, yielding Rydberg excitations on sub-lattice B. d, For the maximum sweep duration, the probability of detecting a Rydberg decay erasure (bars) is significantly increased (decreased) from the baseline level (red fill) at a distance of one (two) away from the site of a detected preparation erasure, induced by the altered Rydberg population on these sites due to the crystal formation in b, indicating direct detection of correlations between errors through erasure imaging.
effect of errors on experimental outcomes at a microscopic level and studying correlations between different error sources, which is enabled by having three separate images for a given experimental run (Fig. 4a).In particular, we consider the joint probability distribution, P(e 3 ), that atoms at sites i, j, and k are detected respectively in the preparation erasure im-age (e 1 ), the decay erasure image (e 2 ) and the final state detection image (e 3 ).
We again consider adiabatic sweeps into the Z 2 phase as in Fig. 3, but now with a total duration of 8 µs.We first study P(e (j+d) 3

= 0|e
(j) 1 = 0), equivalent to finding a Rydberg excitation on site j + d, conditioned on finding no preparation erasure on site j.We plot this quantity (Fig. 4b, left) as a function of both d and the sweep duration.We explicitly average over choices of j and find a signal essentially uniform in d.
However, if we instead consider P(e = 0|e (j) 1 = 1), the probability to find a Rydberg excitation on site j + d conditioned on detecting a preparation erasure on site j, markedly different behavior emerges (Fig. 4b, middle).For simplicity, we further post-select on instances where only a single erasure is detected across the entire array.At intermediate sweep times, we observe an AFM order forms around the preparation erasure error position.We interpret the error as breaking the atom chain into two shorter chains; excitations will naturally form at the system edges of these shorter chains in order to maximize the Rydberg density in the attractive regime in which we operate (Methods).This effectively pins the Rydberg density around the error, which then establishes a preferred AFM order further out into the array.Interestingly, the equivalent quantity for decay erasures, P(e ), shows a more complex behavior.To quantify this behavior more explicitly, we consider a variant of the AFM magnetization (Eq. 1) conditioned on the erasure location, where sublattice A (B) is now defined as being sites an odd (even) distance away from an erasure.In Fig. 4c we plot the mean AFM magnetization for both the preparation (orange circles) and decay erasure (green circles) cases.Preparation erasures develop a negative, single AFM order as they pin Rydberg excitations at odd distances away from the erasure.
Decay erasures behave similarly before the critical point, as Rydberg decay acts effectively as a preparation error.However, past the critical point, this behavior changes: decay now acts as a measurement on the AFM superposition ground state, selecting one of these orders.In this case, assuming perfect Z 2 states, the neighboring sites must have been in the ground state to detect a decay, meaning the AFM order is reversed from the preparation case.This leads to data that first dips to negative values and then grows to positive values past the phase transition (green markers in Fig. 4c).
We also study correlations between preparation errors and Rydberg decay.In particular, a preparation error forces atoms at odd intervals from the preparation erasure to have a higher probability to be in Rydberg states, meaning they should also be more likely to decay.As shown in Fig. 4d, we directly observe this effect at the end of the sweep by considering P(e (j+d) 2

= 1|e
(j) 1 = 1), the probability to detect a decay erasure at a distance d away from a preparation erasure.For d = 1 (d = 2), this probability is significantly increased (decreased) from the unconditional decay erasure probability, in line with the increased (decreased) Rydberg population on these sites, which shows that errors are correlated.
Before concluding, we note that erasure-excision for preparation errors using the first erasure image can be considered heralding the subsequent quantum simulation on the presence of atoms in tweezers in the correct initial state.For erasure-excision of Rydberg decay using the second erasure image, we interpret the post-selected results as coming from a non-jump trajectory in a Monte Carlo wavefunction approach 41 .

Discussion and Outlook
Our results could have broad implications for quantum science and technology.First, our two-qubit entanglement fidelity values and associated error modelling imply that Rydberg arrays, which have already demonstrated scalability to hundreds of atoms 5,6 , can be simultaneously equipped with high-fidelity two-qubit operations, a unique combination across all platforms.Besides our current demonstration of ∼0.999 SPAM corrected two-qubit fidelity, modeling implies that values of ∼0.9997 could be possible with laser noise improvements alone.Further, utilizing a cryogenic environment could freeze out blackbody decay to a large degree 38 , with remaining decay detected as an erasure, leaving almost no intrinsic decoherence.In this context, we note very recent results for improved computational gate fidelities 42 .
Second, the demonstrated erasure conversion techniques could find wide-spread applications for both classical and quantum error correction.For classical correction, our techniques could be modified to correct for state-preparation errors via subsequent atom rearrangement 31,32 , instead of just excising such events.Further, thermal excitations could be converted to erasures and subsequently removed by driving a blue sideband transition between 1 S 0 and 3 P 0 (Fig. 1a) prior to the fast image and subsequent atom-rearrangement 31,32 , effectively realizing erasure-based atomic cooling.
For quantum-error correction, our techniques could be combined with a long-lived qubit which is dark to the fast image, e.g., realized with the 3 P 0 nuclear qubit in neutral Sr 43 and Yb 22,24 , or S 1/2 in Ca + and Ba + ions 10 .Similarly, schemes for implementing erasure conversion in superconducting circuits have been put forward 11,12 .Such techniques could lead to drastically reduced quantum error-correction thresholds 7,8 for faulttolerant quantum computing.
Third, our results also show clearly how NISQ applications 1 can benefit from erasure conversion.Our demonstrated improvements for analog quantum simulation of ground-state physics could be extended to nonequilibrium dynamics, for example targeting regimes generating large entanglement entropies 18 , with the potential to reach a quantum advantage over classical simulations 44 .We note that while our implementation of erasure-excision slows down the effective sampling rate of the quantum device (Ext.Data Fig. 7), the classical cost can increase highly non-linearly with the resulting fidelity increase, and we hence expect a gain for such tasks.Further, we envision erasure-excision improving other tasks such as quantum optimization 45 and potentially quantum metrology 46 .
Finally, insights into erasure-error correlations, as in Fig. 4, could be used to understand error processes in NISQ devices in unprecedented detail, in particular if erasure detection could be made time-resolved with respect to the many-body dynamics.This could also be used to realize post-measurement physics with erasure detection, such as measurement-induced phase transitions 47,48 and measurement-altered quantum criticality 49 .
Note-During completion of this work we became aware of work performing erasure detection with ytterbium atoms 50 .P.S., A.L.S., R.B.T., R.F., and J.C. contributed to the experimental set-up.P.S., A.L.S. and M.E.wrote the manuscript with input from all authors.M.E.supervised this project.

Fast imaging on the erasure detection subspace
Here we describe how we perform the erasure imaging which allows us to detect site-resolved leakage errors 30 .In order to both avoid any extra heating coming from the imaging beams and optimize imaging fidelity, we shine two identical counter-propagating beams with crossed π-polarization and Rabi frequencies Ω/2π ≃ 40 MHz on the 1 S 0 → 1 P 1 transition (see Ext. Data Fig. 1a).This minimizes the net force on an atom, and the crossed polarization avoids intensity interference patterns.
We highlight the characteristic features of this imaging scheme experimentally.We show in Ext.Data Fig. 1b the survival probability of atoms in 1 S 0 as a function of imaging time.After 4 µs, more than 80% of the atoms are lost.However, the number of detected photons continues to increase: even though the kinetic energy of the atoms is too large to keep them trapped, their mean position remains centered on the tweezers.Importantly, for our implementation of erasure-excision, atom loss during the erasure image is inconsequential for our purposes as long as the initial presence of the atom is correctly identified, but in any case, other fast imaging schemes may alleviate this effect 51 .After ∼24 µs, the atomic spread becomes too large and the number of detected photons plateaus.The obtained detection histogram is shown in Ext.Data Fig. 1c.We present the results both for empty (blue) and filled (red) tweezers, which we achieve by first imaging the atoms using usual, high survival imaging for initial detection in a 50% loaded array, then perform the fast image.We obtain a typical detection fidelity of 0.980 +1 −1 of true positives and true negatives, limited by the finite probability for atoms in 1 P 1 to decay into 1 D 2 (see Ext. Data Fig. 1a).
This imaging scheme is sufficiently fast to avoid perturbing atoms in 3 P 0 , as measured by losses from 3 P 0 as a function of imaging time (Ext.Data Fig. 1d).We fit the data (circles) using a linear function (solid line), and obtain a loss of 0.0000046 +12 −12 per image, consistent with the lifetime of the 3 P 0 state 52 of ∼5 s for the trap depth of 45 µK used during fast imaging.
As to the nature of the detected erasure errors for the Bell state generation, we find that preparation errors contribute the vast majority of erasure events compared to bright Rydberg decay, and excising them has a more significant impact on reducing infidelities.In particular, application of Û lasts for only ∼59 ns, which is significantly shorter than the independently measured bright state decay lifetime of 168 +14 −14 µs (Ext.Data Fig. 2).The error model described in Fig. 2 suggests that excising such errors results in an infidelity reduction of only 1.2 +3 −3 × 10 −4 (Methods).Conversely, preparation errors account for ∼5 × 10 −2 infidelity per pair due to the long time between preparation in |g⟩ and Rydberg excitation (Ext.Data Fig. 3).Hence, the gains in fi-delity from erasure-conversion mainly come from nearly eliminating all preparation errors, which has the added benefit of significantly reducing error bars on the SPAM corrected values.Still, SPAM corrected values might also benefit from the small gain in eliminating the effect of bright state decay, and from avoiding potential deleterious effects arising from higher atomic temperature in the repumper case.
For erasure detection employed in the context of manybody quantum simulation, we adjust the binarization threshold for atom detection to raise the false positive imaging fidelity to 0.9975, while the false negative imaging fidelity is lowered to ∼0.6 (Fig. 3d); this is done as a conservative measure to prioritize maximizing the number of usable shots while potentially forgoing some fidelity gains (Ext.Data Fig. 7).
We note that the scheme we show here is not yet fundamentally limited, and there are a number of technical improvements which could be made.First, the camera we use (Andor iXon Ultra 888) has a quantum efficiency of ∼ 80%, which has been improved in some recent models, such as qCMOS devices.Further, we presently image atoms only from one direction, when in principle photons could be collected from both objectives 53 .This would improve our estimated total collection efficiency of ∼ 4% by a factor of 2, leading to faster imaging times with higher fidelity (as more photons could be collected before that atoms were ejected from the trap).Furthermore, the fidelity may be substantially improved by actively repumping the 1 D 2 state back into the imaging manifold so as to not effectively lose any atoms via this pathway.

Details of Rydberg excitation
Our Rydberg excitation scheme has been described in depth previously 13 .Prior to the Rydberg excitation, atoms are initialized from the absolute ground state 5s 2 1 S 0 to the metastable state 5s5p 3 P 0 (698.4nm) through coherent drive.Subsequently, tweezer trap depths are reduced by a factor of 10 to extend the metastable state lifetime.
For Rydberg excitation and detection, we extinguish the traps, drive to the Rydberg state (5s61s 3 S 1 , m J = 0, 317 nm), and finally perform auto-ionization of the Rydberg atoms 13 .Auto-ionization has a characteristic timescale of ∼5 ns, but we perform the operation for 500 ns to ensure total ionization.We report a more accurate measurement of the auto-ionization wavelength as ∼407.89nm.In the final detection step, atoms in 3 P 0 are readout via our normal imaging scheme 13,54 .
Atoms can decay from 3 P 0 between state preparation and Rydberg excitation, which is 60 ms to allow time for magnetic fields to settle.In previous work 13 , we supplemented coherent preparation with incoherent pumping to 3 P 0 immediately prior to Rydberg operations.However, during the repumping process, atoms can be lost due to repeated recoil events at low trap depth, which is not detected by the erasure image, and thus can lower the bare fidelity.Even with SPAM-correction of this effect, we expect the fidelity with repumping to be slightly inferior due to an increased atomic temperature for pumped atoms.

Rydberg Hamiltonian
The Hamiltonian describing an array of Rydberg atoms is well approximated by which describes a set of interacting two-level systems, labeled by site indices i and j, driven by a laser with Rabi frequency Ω and detuning ∆.The interaction strength is determined by the C 6 coefficient and the lattice spacing a. Operators are Xi = |r⟩ i ⟨g| i +|g⟩ i ⟨r| i and ni = |r⟩ i ⟨r| i , where |g⟩ i and |r⟩ i denote the metastable ground and Rydberg states at site i, respectively.For the case of measuring two-qubit Bell state fidelities we set Ω/2π = 6.2 MHz.Interaction strengths in Fig. 2a are directly measured at inter-atomic separations of 4 and 5 µm, and extrapolated via the predicted 1/r 6 scaling to the level at 2.5 µm.Mean atomic distances are calibrated via a laser-derived ruler based on shifting atoms in coherent superposition states 55 .We calibrate C 6 /2π = 230(25) GHz µm 6 using maximum-likelihoodestimation (and associated uncertainty) from resonant quench dynamics 18 , which additionally calibrates a systematic offset in our global detuning.
For performing many-body quasi-adiabatic sweeps, the detuning is swept symmetrically in a tangent profile from +30 MHz to -30 MHz, while the Rabi frequency is smoothly turned on and off with a maximum value of Ω/2π = 5.6 MHz.For an initially positive detuning, the |r⟩ state is energetically favorable, making the all ground initial state, |gg...gg⟩, the highest energy eigenstate of the blockaded energy sector, where no neighboring Rydberg excitations are allowed.For negative detunings, where |g⟩ is energetically favorable, the highest energy state uniquely becomes the symmetric AFM state (|grgr...gr⟩+|rgrg...rg⟩)/ √ 2 in the deeply ordered limit.Thus, considering only the blockaded energy sector, sweeping the detuning from positive to negative detuning (thus remaining in the highest energy eigenstate) is equivalent to the ground state physics of an effective Hamiltonian with attractive Rydberg interaction and inverted sign of the detuning.This equivalence allows us to operate in the effectively attractive regime of the blockaded phase diagram of Ref. 39 .For our Hamiltonian parameters, we use exact diagonalization numerics to identify the infinite-size critical detuning using a scaling collapse near the finite-system size minimum energy gap 56 .

Error modeling
Our error model has been described previously 13,18 .We perform Monte Carlo wavefunction based simulations 57 , accounting for a variety of noise sources including: time-dependent laser intensity noise, timedependent laser frequency noise, sampling of the beam intensity from the atomic thermal spread, Doppler noise, variations of the interaction strength from thermal spread, beam pointing stability, and others.All of the parameters which enter the error model are independently calibrated via selective measurements directly on an atomic signal if possible, as shown in Ext.Data Table 2. Parameters are not fine-tuned to match the measured Bell state fidelity, and the model equally well describes results from many-body quench experiments 18 .

Extraction of the Bell state fidelity
In order to extract the Bell state fidelities quoted in the main text, we use a lower bound method 13 which relies on measuring the populations in the four possible states P gr , P rg , P gg and P rr during a Rabi oscillation between |gg⟩ and |Ψ + ⟩.The lower bound on Bell state fidelity is given by: where P 2π i are the measured probabilities for the four states at 2π, and P π gr+rg is the probability P gr + P rg measured at π.In order to measure these probabilities with high accuracy, we concentrate our data-taking around the π and 2π times (Ext.Data Fig. 5a), and fit the obtained values using quadratic functions f (t) = p 0 + p 1 (t − p 2 ) 2 .We first detail the fitting method, then how we obtain the four probabilities, and finally the extraction of the Bell state fidelity from these.

Fitting method
We perform a fit which takes into account the underlying Beta distribution of the data, and prevents systematic errors arising from assuming a Gaussian distribution of the data.The aim of the fit is to obtain the three-dimensional probability density function Q(p 0 , p 1 , p 2 ) of f , using each experimental data point i defined by its probability density function P i (x), where x is a probability.In order to obtain a particular value of Q(p 0 , p1 , p2 ), we look at the corresponding probability density function value P i (f (t i )) for each data point i, where f (t i ) = p0 + p1 (t i − p2 ) 2 , and assign the product of each P i (f (t i )) to the fit likelihood function: We repeat this for various [p 0 , p1 , p2 ].
The result of such fitting method is shown in Ext.Data Fig 5b (black line), where we present f (t) = p 0 + p 1 (t − p 2 ) 2 for [p 0 , p 1 , p 2 ] corresponding to the maximum value of Q(p 0 , p 1 , p 2 ).We emphasize that this results in a lower peak value than a standard fitting procedure which assumes underlying Gaussian distributions of experimentally measured probabilities (red line).Choosing this lower peak value eventually will provide a more conservative, but more accurate value for the Bell state fidelity lower bound than the naïve Gaussian approach.

Obtaining the four probability distributions
Our method to obtain the probability density functions of the four probabilities at π and 2π times ensures both that the sum of the four probabilities always equals one, and that their mutual correlations are preserved.We first extract the Beta distribution of P rr by gathering all the data around the π and 2π times (Ext.Data Fig. 5c).In particular, the mode of the obtained Beta distribution at π is P rr ≃ 0.0005.The distribution of P gr+rg and P gg are obtained by fitting the data in the following way.We perform a joint fit on P gr+rg using a fit function f 1 (t), and on P gg using a fit function f 2 (t).The fit functions are expressed as: which ensures that the sum of the four probabilities is always equal to 1.We then calculate the joint probability density function Q 1,2 (p 0 , p 1 , p 2 ) of both f 1 and f 2 using the method described above.In particular: where P gr+rg i (P gg i ) is the probability density function associated with P gr+rg (P gg ) for the i-th experimental data point.In particular, we impose that p 0 ≤ 1 − P rr to avoid negative probabilities.We show the resulting Q 1,2 (p 0 , p 1 , p 2 ) in Ext.Data Fig. 5d as two-dimensional maps along (p 0 , p 1 ) and (p 0 , p 2 ).
We then obtain the one dimensional probability density function for p 0 by integrating over p 1 and p 2 (see Ext. Data Fig. 5d).This provides the fitted probability density function of P gr+rg , and hence P gg = 1 − P rr − P gr − P rg at π time.We repeat this process for various values of P rr , for both π and 2π times.
At the end of this process, we obtain different probability density functions for each P rr value.The asymmetry between P gr and P rg is obtained by taking the mean of P gr − P rg at π and 2π times.We assume the underlying distribution to be Gaussian, as P gr − P rg is centered on 0, and can be positive or negative with equal probability.

Bell state fidelity
Now that we have the probability density function for all four probabilities at π and 2π times, we move on to the Bell state fidelity extraction.For both π and 2π, we perform a Monte-Carlo sampling of the Beta distribution of P rr which then leads to a joint probability density function for P gr+rg and P gg .We then sample from this, and use Eq. 3 to obtain a value for the Bell state fidelity lower bound.We repeat this process 1 million times, and fit the obtained results using a Beta distribution (see Ext. Data Fig. 5e).We observe an excellent agreement between the fit and the data, from which we obtain F Bell ≥ 0.9962 +10 −13 , where the quoted value is the mode of the distribution, and the error bars represent 68% confidence interval.
We use the same method to obtain the measurementcorrected Bell fidelity, and the state preparation and measurement (SPAM) corrected one.After drawing the probabilities from the probability density functions, we infer the SPAM-corrected probabilities from our known errors, described in detail previously 13 .We use the values reported in Ext.Data Table 2.During this process, there is a finite chance that the sum of probabilities does not sum up to one.This comes from the fact that the probability density functions and the SPAM correction are uncorrelated, an issue which is avoided for raw Bell fidelity extraction thanks to the correlated fit procedure described above.We employ a form of rejection sampling to alleviate this issue by restarting the whole process in the case of such event.We perform this 1 million times, and fit the obtained results using a Beta distribution (see Ext. Data Fig. 5f).We observe an excellent agreement between the fit and the data, from which we obtain a SPAM-corrected fidelity F Bell ≥ 0.9985 +7 −12 , where the quoted value is the mode of the distribution, and the error bars represent 68% confidence interval.

Interaction limitation for Bell fidelity
We estimate the theoretically expected Bell state fidelity using perturbation analysis.Specifically, the resonant blockaded Rabi oscillation for an interacting atom pair is described by the following Hamiltonian where V = C 6 /r 6 is the distance-dependent, interaction strength between two atoms separated at distance r (see Eq. 2).Since the two-atom initial ground state, |ψ(0)⟩ = |gg⟩, has even parity under the left-right reflection symmetry, the Rabi oscillation dynamics can be effectively solved in an even-parity subspace with three basis states of |gg⟩, |rr⟩, and In the Rydberg-blockaded regime where V ≫ Ω, we can perform perturbation analysis with the perturbation pa-rameter η = Ω/ √ 2V and find that the energy eigenvectors of the subspace are approximated as , and E 3 ≈ V (1 + η 2 /2), respectively.Rewriting the initial state using the perturbed eigenbasis, we solve to obtain the analytical expression of the maximum achievable Bell state fidelity, F Bell , at a given perturbation strength η.Keeping the solution up to the second order of η, we find obtained at t = π/ √ 2Ω.

Statistics reduction due to erasure-excision
Our demonstration of erasure-excision explicitly discards some experimental realizations (Ext.Data Fig. 6), which can be seen as a downside of the method.However, this is a controllable trade-off: by adjusting the threshold for detecting an erasure error, we can balance gains in fidelity versus losses in experimental statistics (as shown in Ext.Data Fig. 7) for whatever particular task is of interest.In general, the optimum likely always includes some amount of erasure-excision, as it is usually better to remove erroneous data than keeping them.
Ext. Data Fig. 1 | Fast imaging on the erasure detection subspace.a, Sketch of the involved energy levels.We detect atoms in 1 S0 by strongly driving the 1 S0 ↔ 1 P1 transition.b, Survival of atoms in 1 S0 (green) and number of detected photons (red) as a function of the imaging time.We observe an increase of detected photons whereas the atoms are already lost: even though the kinetic energy of the atoms is too large to keep them trapped, their mean position remains centered on the tweezers thanks to the use of two counter-propagating beams with equal power.After ∼24 µs, the atomic spread becomes too large to measure a significant increase in detected photons.c, Typical histograms of the number of detected photons for 24 µs imaging.Using a slow, high-fidelity image prior to the fast image, we can detect if a tweezer is empty (blue) or filled (red).The typical detection fidelity which corresponds to equal error probability in detecting absence or presence of an atom is 98.0(1)%.d, Losses from 3 P0 as a function of time, expressed in number of fast images.The survival probability of an atom in 3 P0 is 99.99954(12)% for one image, consistent with its 5 second lifetime.We present the results with (blue) and without (red) performing the erasure imaging during the holding time.We observe no significant difference between the two conditions, which suggests that the erasure detection imaging light, in principle, does not destroy the coherence of the Bell state.
Ext. Data Fig. 7 | Controlling fidelity gains from erasure-excision.Erasure-excision explicitly trades improved experimental fidelity for a reduced number of experimental repetitions.However, this process is controllable by adjusting the threshold used for detecting atoms during erasure images.Changing the threshold essentially changes the false positive and false negative rate for detecting erasures correctly (inset).We plot the total AFM probability (green markers) after a sweep as in Fig. 3 of the main text, and vary the detection threshold used for identifying erasures.For too high a threshold, many erasure events go unnoticed, and so erroneous outcomes become relatively more prevalent, reducing the overall fidelity.As the detection threshold is lowered, more true erasures, where an atom is actually present in the erasure image, are correctly detected, which improves the fidelity.However, lowering the threshold too far (in our case past ∼5 photons) increases the likelihood of seeing false positive erasures; excising data based on these events discards experimental statistics with relatively little gain in fidelity.In the main text, we select a detection threshold of 5 photons.

Fig.
Fig.1| Erasure conversion for high-fidelity entanglement.a, Level structure used in this work.We distinguish two subspaces: a qubit subspace in which the atoms interact via their Rydberg states, and a measurement subspace used to detect leakage errors from the qubit subspace with single-site resolution, realizing erasure conversion.b, Sketch of the erasure conversion scheme, as applied to Bell pair generation.After arranging atoms into pairs (top) we prepare them in |g⟩, and entangle them via the Rydberg blockade mechanism (right), denoted by a unitary operation Û (t).Immediately afterwards, we auto-ionize atoms in |r⟩, effectively projecting the populations of the Bell states, and follow with a fast erasure conversion image to detect leakage out of the qubit subspace during the preparation or evolution periods.This is followed by the final detection of atoms in |g⟩, yielding two separate, independent images.We can discard data from pairs where atoms are detected in the erasure-error image, termed erasure-excision in the following.Atom fluorescence images are single-shot, with post-processing applied to improve detection fidelity 30 (Methods).c, Lower bounds for Bell state fidelities with (blue) and without (red) the erasure-excision, and using incoherent pumping to reduce preparation errors instead of erasure-excision (green, Methods).We present the results for the raw data, corrected for measurement errors, and corrected for state preparation and measurement errors (SPAM).All data is averaged over eight pairs of atoms which are excited in parallel.Error bars represent a 68% confidence interval (Ext.Data Fig.5, Methods).

Fig. 3 |
Fig.3| Erasure conversion in quantum simulation.a, b, We perform quasi-adiabatic sweeps with N = 26 spins in the effective ground state manifold of an attractive Rydberg Hamiltonian (Methods), starting from the initially disordered phase, and ending in the Z2 ordered phase.c, We perform two erasure images, one preceding the quantum simulation (to detect preparation errors), and one following (to partially detect Rydberg decay errors).d, The probability for detecting a preparation error (orange markers) remains constant while the probability for detecting a decay error (green markers) grows in agreement with the Rydberg lifetime and detection infidelities (solid lines, Methods).e, The total probability, PAFM, for forming either of the AFM Z2 states is improved by performing erasure-excision on all errors (blue markers), as compared to only on preparation errors (orange markers) or performing no excision (pink markers).The sensitivity of PAFM with respect to a change in ∆ also increases with erasure-excision (inset).f, The probability distribution for measuring a given AFM magnetization is initially peaked at 0 in the disordered phase, before bifurcating when entering the Z2 phase, consistent with spontaneous symmetry breaking.g, Deep in either phase, erasure-excision leads to a sharpening of the probability distribution (left, right).Around the phase transition, we observe a close-to-flat distribution (middle).

Fig. 4 |
Fig.4| Learning from erasure errors.a, Post-processed30 single-shot atom fluorescence images (Methods).After arranging the array (top panel) we use the fast erasure images to learn how detected errors (middle panels) affect detection outcomes (lower panel).b, Conditional probability to detect no atom in the final image as a function of sweep progress and distance from a hypothetical erasure event.In the case of no erasures (left), the resulting profile is uniform.However, when conditioning on detecting a preparation erasure (middle), the error breaks the Z2 symmetry by establishing a single AFM order.In the case of conditioning on decay errors (right) the situation is more complex.c, AFM magnetization (see Eq. 1) as a function of time.Preparation erasures (orange markers) lead to a growth of a single AFM order with Rydberg excitations predominantly on sub-lattice A (defined as sites an odd distance from the erasure position).Decay erasures (green markers) follow a similar trend at early times by acting as effective preparation errors, but past the critical point (vertical dashed line) their behavior reverses: a decay erasure spontaneously breaks the two-fold symmetry, where the neighboring sites must have been in the ground state due to the Z2 structure, yielding Rydberg excitations on sub-lattice B. d, For the maximum sweep duration, the probability of detecting a Rydberg decay erasure (bars) is significantly increased (decreased) from the baseline level (red fill) at a distance of one (two) away from the site of a detected preparation erasure, induced by the altered Rydberg population on these sites due to the crystal formation in b, indicating direct detection of correlations between errors through erasure imaging.

Ext. Data Fig. 2 |
Rydberg lifetime.a, Pseudo-level diagram of the Rydberg dynamics and associated decay channels.During Rydberg evolution with Rabi frequency Ω, the Rydberg atom can decay into either a set of states which is 'bright' to the imaging process (including both the erasure images and the final detection images), e.g states like 5s5p 3 P2, or into states which are 'dark' to the imaging process, e.g.nearby Rydberg states.A small percentage of decays into 'bright' states can go directly into 3 P0 where they can be re-excited by the Rydberg driving; note that such decays are dark in the erasure image, but bright in the final detection image.b, Measurement of the dark state decay lifetime, measured by performing a π-pulse on the Rydberg transition, waiting a variable amount of time, and then returning atoms to the ground state (inset).c, Measurement of the bright state decay lifetime, measured by performing a Rydberg π-pulse, waiting, and then performing an auto-ionization pulse to destroy any remaining Rydberg or dark state excitations.d, We prepare atoms into 3 P2 (a bright state), and then perform continuous Rydberg driving.Atoms are lost from the trap at a rate which increases with increasing Rabi frequency.e, The lifetime of atoms in 3 P2 scales inversely with the square of the Rabi frequency (i.e.scales inversely with the intensity of the Rydberg beam).We attribute this to a photo-ionization process which can convert bright state decay into dark state decay through prolonged Rydberg excitation, as shown in a. Markers in b-d are experimental data where error bars are often smaller than the marker sizes, and solid lines represent exponential fits.

Ext. Data Fig. 3 |
Experimental sequence.Sketch of the experimental sequence including the erasure detection for a the Bell state generation experiment, and b the many-body experiment.Both experiments have the same global architecture: we start by loading the atoms into the desired geometry, then initialize the atoms in |g⟩, perform the Bell state generation or quantum simulation, and finally read out by auto-ionizing atoms in |r⟩ and imaging atoms in |g⟩.The main difference between both experiments concerns the erasure detection.In a, we utilize a single erasure detection, placed after auto-ionizing atoms in |r⟩.In b, we perform two erasure images: one before applying Û (t), and one after auto-ionization.Ext.Data Fig. 4 | Coherence preservation during erasure detection.a,b, We prepare the Bell state |Ψ + ⟩, and measure the relevant populations for Bell fidelity extraction Pgr + Prg and Pgg (after an extra π pulse) as a function of holding time.We perform a 2π pulse in the middle of the holding time to get rid of dephasing effect due to e.g.Doppler effect.