## Optimizing quantum error correction protocols with erasure qubits

Shouzhen Gu, 1, 2, ∗ Yotam Vaknin, 2, 3, ∗ Alex Retzker, 2, 3 and Aleksander Kubica 2, 1, 4

1 California Institute of Technology, Pasadena, CA 91125, USA 2 AWS Center for Quantum Computing, Pasadena, CA 91125, USA

3 Racah Institute of Physics, The Hebrew University of Jerusalem, Jerusalem, 91904, Givat Ram, Israel 4 Department of Applied Physics, Yale University, New Haven, CT 06511, USA (Dated: August 28, 2024)

Erasure qubits offer a promising avenue toward reducing the overhead of quantum error correction (QEC) protocols. However, they require additional operations, such as erasure checks, that may add extra noise and increase runtime of QEC protocols. To assess the benefits provided by erasure qubits, we focus on the performance of the surface code as a quantum memory. In particular, we analyze various erasure check schedules, find the correctable regions in the phase space of error parameters and probe the subthreshold scaling of the logical error rate. We then consider a realization of erasure qubits in the superconducting hardware architectures via dual-rail qubits. We use the standard transmon-based implementation of the surface code as the performance benchmark. Our results indicate that QEC protocols with erasure qubits can outperform the ones with state-of-theart transmons, even in the absence of precise information about the locations of erasure errors.

## I. INTRODUCTION

Quantum error correction (QEC) with erasure qubits has been a promising avenue toward reducing the overhead associated with fault-tolerant quantum computation. In an erasure qubit, the dominant noise removes the state from the computational subspace. Information from the detection of such errors, which we refer to as erasures, can be effectively used to improve the performance of QEC protocols. While the concept of erasures via the dualrail construction had a central role in photonic quantum computation [1-5], various realizations of QEC with erasure qubits have recently been proposed for other platforms, including neutral atoms [6], trapped ions [7] and superconducting circuits [8-10], with several promising experimental demonstrations [1115]. Erasure qubits can offer improvements for arbitrary QEC codes, although recent works have focused on quantifying the benefits for topological codes, such as the surface code [16] or the honeycomb code [17].

Compared to standard QEC protocols, using erasure qubits comes at a cost. Namely, quantum circuits need to include additional elements, erasure checks (ECs) and reset operations, that are capable of diagnosing erasures and reinitializing erased qubits back to the computational subspace. These operations take additional time during which noise may accumulate. From this perspective, one would prefer to implement as few of them as possible. On the other hand, ECs provide valuable information about the location of possible erasures. Therefore,

∗ These authors contributed equally.

one expects an optimal frequency of ECs that balances these two effects. Most of previous works have assumed that ECs are performed after every entangling operation without optimizing their frequency.

In this article, we assess the benefits for QEC protocols provided by erasure qubits. Using the framework for QEC with erasure qubits introduced in Ref. [18], we analyze the performance of the surface code as a quantum memory. For various EC schedules, we obtain the correctable region in the phase space of error parameters; see Fig. 1. We also analyze the subthreshold scaling of the logical error rate. In our simulations, we use a novel decoding method that relies on an approximate conversion of QEC protocols with erasure qubits into stabilizer circuits, which may be of independent interest.

We then benchmark the performance of QEC protocols with dual-rail qubits, which are erasure qubits composed of two coupled transmons [9, 14, 19], against the standard T 1 -limited tramsmons. We observe that using dual-rail qubits results in an approximately 50% increase in QEC thresholds, provided that the measurement times of dual-rail qubits and transmons are comparable and there is a significant bias between the amplitude damping noise and Pauli errors. A key factor that indicates the advantageous performance of dual-rail qubits is the ratio of the measurement time T M to the two-qubit (2Q) gate time T 2 Q -the smaller the ratio T M /T 2 Q , the better for dual-rail qubits. We also analyze different implementations of reset operations. We compare a one-way pulse, which brings the erased state back to the computational subspace without dephasing the computational subspace, with a unitary pulse, which transforms between the erased state and the computational subspace. We find that the most effective protocol employs a one-way pulse exclusively to the

FIG. 1. The correctable region in the ( e, p, q ) phase space for the surface code with erasure qubits and EC schedules with: (a) 4 ECs, (b) 2 ECs, and (c) 1 EC, where e , p , and q denote the erasure, Pauli, and measurement error rates, respectively. Cross sections of the threshold surfaces in the ( e, p ) plane for the different EC schedules when: (d) q = 0, (e) q = 0 005, and (f) . q = 0 01. . The shaded regions indicate parameter regimes where only one schedule is below threshold. The dashed line corresponds to an optimistic scenario where ECs and reset operations do not cause additional errors. The dashed-dotted line is obtained from a scheme with noiseless reset operations but no ECs, giving an upper bound for schemes that do not use information about the locations of erasures.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-36d4b85b-dd89-4ebc-806d-1b400f236752.jpg)

qubits identified as erased, and it reduces the impact of measurement errors by approximately a factor of 2.

## A. Erasure qubits

The remainder of the paper is organized as follows. In Sec. II, we review the formalism for QEC protocols with erasure qubits and describe the surface code as a quantum memory. In Sec. III, we describe an approximate solution to the decoding problem with erasure qubits, as well as analyze the threshold surface and subthreshold performance of the different EC schedules. While the results of Sec. III are applicable to any implementation of erasure qubits, Sec. IV focuses on superconducting qubits. We compare dual-rail qubits with transmons, and also analyze the different reset operations. Possible future directions are discussed in Sec. V.

## II. PRELIMINARIES

We first review the formalism to describe QEC protocols with erasure qubits. Then, we explain how to implement the surface code with erasure qubits.

We use the formalism introduced in Ref. [18] to describe and simulate QEC protocols with erasure qubits. Compared to standard qubits, two additional operations are used when working with erasure qubits: ECs that measure if a qubit is erased, and reset operations that reinitializes a qubit in the computational subspace. For our purposes, a stabilizer circuit is one that consists of the following basic operations: single-qubit (1Q) state preparation in an eigenstate of a Pauli operator, 1Q readout in a Pauli basis, 1Q Clifford gates, and 2Q controlledPauli gates. When a stabilizer circuit is implemented with erasure qubits and enhanced with the two additional operations, we call it an erasure circuit .

To simulate erasure circuits, we associate noise sources with each operation. We do this by adding erasure locations between every operation. Then, we append Pauli noise P after every operation and include bit flip noise N on classical outcome bits (see Table I). We use e to denote erasure rates at erasure locations, p to denote Pauli error rates, and q to denote classical bit-flip rates. The noise strengths may vary depending on the operations around the error location. In a simplified model, we may consider

TABLE I. Mapping of an ideal circuit to a simulated circuit. The error channel P ( p ) is the 1Q or 2Q depolarizing channel with error rate p , and N ( q ) is the binary symmetric channel (that flips the measurement outcome) with error rate q . In general, P and N can represent arbitrary Pauli and binary channels.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-4b61ab12-5ae3-4ddd-8646-e562ad45ba09.jpg)

| operation                | ideal   | simulated            |
|--------------------------|---------|----------------------|
| state preparation        | | ψ ⟩   | | ψ ⟩ P ( p )        |
| readout                  | P       | EC P N ( q ) N ( q ) |
| erasure check with reset | EC ∗    | EC R P ( p ) N ( q ) |
| 1Q gate                  | G       | G P ( p )            |
| 2Q gate                  | G       | G P ( p )            |

Pauli channels to be depolarizing and e , p , and q to be constant within the circuit.

At an erasure location, the qubit has a probability e of being taken to an erasure subspace that is orthogonal to the computational subspace. We assume that any computational subspace measurement involving an erased qubit results in a randomized outcome. When an erased qubit interacts with another qubit through an entangling gate, we assume that a fully depolarizing channel is applied to the other qubit, unless stated otherwise; see Sec. IVD for a physical model where erasure spreads via dephasing noise. A reset operation reinitializes an erased qubit as the maximally mixed state in the computational basis while acting trivially on the computational subspace. Alternatively, we may consider a unitary reset operation that exchanges the erased state with a fixed state in the computational basis. Reset operations may be conditioned on the classical EC outcomes. We compare different reset schemes in Sec. IV C.

In our formalism, no coherences are created between the computational and erasure subspaces. 1 This allows us to efficiently sample from erasure circuits. We use one bit of information per qubit to

1 We approximate the state after a unitary reset operation as an incoherent mixture.

keep track of its erasure state and update it appropriately at erasure locations and reset operations. For qubits in the computational subspace, we use the Gottesman-Knill theorem [20, 21] to represent the state and simulate the stabilizer circuit operations.

## B. Surface code as quantum memory

The surface code is one of the most studied quantum error-correcting codes. As a CSS stabilizer code, its codespace is defined as the +1 eigenspace of a set of commuting Pauli X and Z operators associated with the faces of a square lattice; see Fig. 2(a). To measure Pauli X and Z operators, we use extra ancilla qubits and implement standard syndrome extraction circuits, as depicted in Fig. 2(b). Importantly, the order of controlled-Pauli gates is chosen to minimize the error spread and to avoid hook errors. By repeating measurements of Pauli X and Z operators, we can learn reliable information about the errors afflicting the system. Consequently, we can find an appropriate recovery operator and protect the encoded logical information from the errors.

FIG. 2. (a) The surface code with code distance d = 5. The ancilla qubits (green dots) are used to measure stabilizer checks on the adjacent data qubits (black dots). Grey and white regions are associated with X and Z stabilizers, respectively. (b) Syndrome extraction circuit for the surface code, where the locations of possible ECs with reset are labeled by A B C , , , and D .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-c1e98774-334e-4b97-89d2-81d8a2b35f37.jpg)

Implementing the surface code with erasure qubits requires including ECs and reset operations in the syndrome extraction circuit. A natural schedule, which we call the 1 EC schedule, is to perform ECs

and reset operations on data qubits at locations D in Fig. 2(b) while ancillas are being measured. Note that erasure information about the ancilla qubits will also be learned during the readout process. If ECs with reset take similar time as measurement and reinitialization of ancilla qubits, thus assuming operations are coherence limited this schedule will incur neither time nor noise penalty, while providing information about erasures.

However, performing only one EC per qubit during each syndrome extraction round may provide too little information about erasures, annulling the possible benefits of erasure qubits. We therefore consider two additional schedules with more ECs and reset operations. The 2 EC schedule has two ECs with reset per syndrome extraction round placed at locations B and D in Fig. 2(b), and the 4 EC schedule places ECs with reset at all four locations A , B , C , and D . Note that performing two or more ECs with reset per syndrome extraction round will lengthen the syndrome extraction round compared to the standard approach.

## III. ARCHITECTURE-INDEPENDENT RESULTS

In this section, we present our results on implementing the surface code with erasure qubits, without restricting to any particular architecture. The simulations are performed by initializing an eigenstate of a logical Pauli operator for the surface code and sampling using the method outlined in Sec. II A. At the end of the simulation, we decode by approximating the erasure circuit as a stabilizer circuit and applying minimum-weight perfect matching on the resulting decoding graph (see Sec. III A). An error is reported if after decoding, the value of the logical Pauli operator measured at the end of the circuit is different than the initial value. Running the experiment many times for a distanced code gives us the logical error rate per d syndrome extraction rounds, denoted p L . For more simulation details, see Appendix B.

## A. Approximate solutions to the decoding problem

Running a QEC protocol with erasure qubits requires correcting the errors that occurred using the available syndrome and EC outcomes. We propose a variation of the decoding method for erasure circuits in Ref. [18] that is suitable for schedules with infrequent ECs.

Erasure circuits can be decoded in three steps.

In the first step, we use the EC outcomes to map the erasure circuit to a stabilizer circuit with independent error mechanisms which are binary random variables. This allows us to decode using standard techniques for stabilizer circuits. In particular, we obtain a decoding hypergraph with error mechanisms as hyperedges and detectors as vertices. By definition, detectors correspond to products of measurement outcomes that are deterministic in the absence of qubit and measurement errors. Each hyperedge in the decoding hypergraph consists of all of the detectors that would be flipped if the error mechanism occurs. In the second step, we approximate the hypergraph as a graph by decomposing the hyperedges into edges, i.e., approximating error mechanism with ones that cause at most two violated detectors. In the third step, we use minimumweight perfect matching on the graph from the previous step to find a likely error causing the syndrome.

Ref. [18] described a method to do the first step exactly. Although the number of error mechanisms added in the conversion process is proportional to the size of the circuit if reset operations occur on every qubit at constant time intervals, that constant is exponential in the length of those intervals. For example, with four entangling gates between reset operations, the 1 EC schedule would introduce 1023 different error mechanisms per reset operation, making decoding impractical.

Here, we introduce a way to approximately convert erasure circuits to stabilizer circuits that results in fewer added error mechanisms. The idea is that the large number of error mechanisms is due to erasures causing correlated depolarization at many spacetime locations. Instead, we can approximate the converted stabilizer circuit by one with independent 1Q depolarizing errors, but still maintaining the marginal probability of error at each spacetime location. The number of error mechanisms added to the resulting decoding hypergraph is only linear in the number of gates between reset operations.

We proceed as follows on an erasure circuit C E . Consider a segment s of a qubit q , which is defined to be the worldline of q between two consecutive reset operations (see Fig. 3(a)). We will remove the erasure locations, ECs, and reset operations in s and add appropriate 1Q depolarizing channels to approximate the effects of erasures in s . Doing this for all segments, we convert C E to a stabilizer circuit C .

For i ∈ { 1 , . . . , r } , let G i be the i -th entangling gate in s and F i be the spacetime location placed after G i on the other qubit. For convenience, let G r +1 denote the second reset operation in s and F r +1 be the location at q after the second reset operation. (See Fig. 3(b) for an illustration.) If q was erased before G i , the operation G i will result in

FIG. 3. (a) A segment s of the qubit q in an erasure circuit C E . (b) The corresponding portion of the converted stabilizer circuit C . Depolarizing channels are placed at spacetime locations F i according to Algorithm 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-74449614-c68a-4c29-ba74-6108da271162.jpg)

full depolarization at spacetime location F i . Let ¯ a i be the probability of this event, conditioned on the EC outcomes. The distribution of errors introduced by erasures in s is approximated by independently sampling events with probabilities ¯ a , . . . a 1 ¯ r +1 and placing fully depolarizing channels at spacetime location F i if the event with probability ¯ a i is sampled. In other words, 1Q depolarizing channels of error probability 3¯ a / i 4 are placed at locations F i . Note that each of these depolarizing channels is equivalent to three independent error mechanisms which apply each of the nontrivial Pauli operators X Y , and Z with probability 1 2 ( 1 -√ 1 -¯ a i ) [22].

Algorithm 1 Approximate conversion of an erasure circuit to a stabilizer circuit with independent error mechanisms

## Input:

erasure circuit C E , erasure check outcomes ⃗ d

## Output:

stabilizer circuit C , error mechanisms { ( P i,j , p i ) }

1:

S

←{

segments in

C

E

}

- 2: for each s ∈ S do
- 3: { G i } ← entangling gates in s
- 4: {F } ← i spacetime locations associated with s

$$4 \colon \quad & \{ \mathcal { F } _ { i } \} \leftarrow \text{ spacetime locations associated with } s \quad \text{ supp} \\ 5 \colon \quad & \text{ for each } i \, \text{ do} \\ 6 \colon \quad & \bar { a } _ { i } \leftarrow \text{ Pr} \left ( q \, \text{ erased before } G _ { i } \Big | \vec { d } \right ) \quad & \text{ old } s \Big | \\ 7 \colon \quad & \text{ p} _ { i } \leftarrow \frac { 1 } { 2 } \left ( 1 - \sqrt { 1 - \bar { a } _ { i } } \right ) \quad & \text{ ment} \\ 8 \colon \quad & \text{ for each } \text{ nontrivial Pauli error } P _ { i, j } \, \text{ at } \mathcal { F } _ { i } \, \text{ do} \quad & \text{ senti} \\ 9 \colon \quad & \text{ include error mechanism } \left ( P _ { i, j }, p _ { i } \right ) \quad & \text{ frequ} \\ 10 \colon \quad & \text{ end for each } \quad & \quad \text{ we hi} \\ 11 \colon \quad & \text{ and for each } \quad & \quad \text{ pr} \quad \\ \end{cases}$$

11: end for each

12: end for each

13: C ← C E with deleted erasure checks and reset

14:

return

C

,

{

(

P

i,j

, p

i

)

}

The approximate conversion of an erasure circuit to a stabilizer circuit with independent error mechanisms is summarized in Algorithm 1. See Appendix A for example highlighting the difference between this method and the exact conversion of Ref. [18].

We remark that it is not straightforward to compare decoding using the exact conversion of Ref. [18] with using the approximate method in this section. Although the first step is more accurate using exact conversion, the resulting hyperedges of the decoding hypergraph may be decomposed less optimally (by, for instance, Stim [23]) in step two of the procedure. This hypothesis is explored in Fig. 4, where we highlight the difference in performance between the two decoding methods. In both cases, we decode using the three steps outlined in this section, the only difference being the exact conversion of Ref. [18] or the approximate conversion outlined here in step one. We see that the latter method performs just as well or even better than the former. The argument that this is due to step two is strengthened by the more pronounced performance gap in the 2 EC schedule compared to the 4 EC schedule. Because sparser reset operations results in larger correlations, it could result in larger inaccuracies from suboptimally decomposing bigger hyperedges into edges. Although the inefficiency of the first method prevents us from performing a similar comparison for the 1 EC schedule, the results for the 4 EC and 2 EC schedules demonstrate the effectiveness of the approximate circuit conversion method. In what follows, we choose to decode using the approximate conversion because it achieves similar (or lower) logical error rates while becoming computationally more efficient.

## B. Probing the correctable region

In Fig. 1(a)-(c), we plot the correctable region in the ( e, p, q ) phase space for the different EC schedules. The threshold surfaces bound the correctable regions where logical error rates are exponentially suppressed with increasing system size.

Fig. 1(d)-(f) shows cross sections of the threshold surfaces in the ( e, p ) plane for various measurement error rates q , with the different curves representing different EC schedules. As we increase the frequency of ECs, the e threshold increases because we have more information about erasures, while the p threshold decreases because there are more sources of noise. We also plot the threshold curve (dashed line) for a scenario where the ECs with reset do not induce extra noise. This optimistic situation arises if ECs are instantaneous or performed at the same time as entangling gates. The dashed-dotted line is an upper bound for standard QEC schemes that do not perform ECs. Here, we assume that leaked qubits are periodically reset but that the reset operations do not introduce noise in the system. The

FIG. 4. Comparison of decoding methods for the (a) 4 EC and (b) 2 EC schedules. The solid lines correspond to approximate conversion of the erasure circuit to a stabilizer circuit; the dashed lines correspond to the exact method. The logical error rates are plotted for different values of x along the line ( e, p, q ) = ( x, x/ 10 , x ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-658267eb-d2cb-4e46-bc45-35f77d18edda.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-89c24912-7fe1-46b0-987f-b045c86cf99c.jpg)

data indicates that the correctable region for erasure protocols almost entirely contains that of any standard scheme, and there is a significant regime of high erasure that is only correctable with erasure schemes.

It is important to note that defining EC as a separate circuit component results in different numerical values of the Pauli and erasure thresholds compared to earlier studies [8]. In those studies, the erasure rate is specified per CNOT gate without accounting for additional erasure mechanisms. Here, since each qubit can be erased independently and may be erased during both entangling gates and ECs, there is a roughly 4 × difference between the reported values caused solely by different definitions. For a more definition-independent analysis, we compare erasure and non-erasure qubits in Sec. IV B.

## C. Subthreshold performance

Next, we explore the subthreshold scaling of the logical error rate. Knowing the functional form p L would allow us to estimate the minimal distance needed to achieve any given logical error rate for a set of physical parameters ( e, p, q ). For low physical error rates, the probability of logical failure is dominated by the most likely error configurations [24]. In the phenomenological setting, a heuristic for the logical error rate under Pauli noise is p L ∝ ( p/p ∗ ) ⌈ d/ 2 ⌉ , where p ∗ is the Pauli noise threshold. Similarly, pure erasure noise would give the scaling p L ∝ ( e/e ∗ ) d , where e ∗ is the erasure threshold.

In our setting, we wish to determine the logical error rate for any point ( e, p, q ) in the correctable region bounded by the threshold surface. We propose an ansatz of the following form. Let x be a single parameter characterizing a line ℓ through the origin ( e, p, q ) = ( t e x, t p x, t q x ) for constants t e , t p , t q ; see Fig. 5(a). The line ℓ describes different physical error rates of a given noise bias. Within the correctable region, we expect the logical error rate to scale along ℓ as

$$p _ { L } \in a \left ( \frac { x } { x ^ { * } } \right ) ^ { \alpha d },$$

where x ∗ is the value of x at the intersection of ℓ with the threshold surface. The quantities a , α are fitting parameters that depend on ℓ . The parameter a is the logical error rate at the threshold, which we find to be roughly constant for all lines ℓ and EC schedules. Importantly, α characterizes the effective distance of the scheme, interpolating between 1 / 2 for pure Pauli noise and 1 for erasure noise in the phenomenological setting.

To test how well Eq. (1) describes the subthreshold logical error rate scaling, we plot p L against d for different values of x along a certain line ℓ ; see Fig. 5(b). The solid lines, indicating the expected logical error rates from the ansatz by fitting a and α , do approximate the data closely. Values of α for different noise biases are reported in the table Fig. 5(c). We see that α is close to 1 / 2 when the erasure bias, defined as the ratio e/p = t e /t p , is low, as expected for unknown Pauli noise. Its value increases with the erasure bias of ℓ . From a theoretical standpoint, in the infinite erasure bias limit and in the absence of measurement errors, α should approach 1 for the 4 EC schedule because a single erasure cannot spread to Pauli errors on more than one data qubit. In contrast, α is lower for the less frequent EC schedules which allow more spread of erasures. The plots showing how well the ansatz in Eq. (1) describes the data used to obtain the α values in the table are presented in Fig. 12 of Appendix B.

| ( t e , t p , t q )   |   4 EC |   2 EC |   1 EC |
|-----------------------|--------|--------|--------|
| (1 , 0 . 01 , 1)      |   0.84 |   0.76 |   0.64 |
| (1 , 0 . 1 , 1)       |   0.72 |   0.69 |   0.6  |
| (1 , 0 . 5 , 1)       |   0.61 |   0.59 |   0.57 |
| (1 , 5 , 5)           |   0.51 |   0.51 |   0.51 |

FIG. 5. The subthreshold scaling of the logical error rate according to the ansatz of Eq. (1). (a) We test logical error rates along the orange line ℓ in the ( e, p, q ) phase space. The orange dot is at the intersection of ℓ with the threshold surface, giving the threshold value x ∗ . (b) The logical error rates with increasing distance for various values of x along ℓ , shown in the different colors. To fit the data we use the numerical ansatz p L = a ( x/x ∗ ) αd with fitting parameters a , α . (c) Values of α for various noise biases in the different EC schedules.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-709146de-3e23-4c8e-9a64-760b1938300a.jpg)

In Fig. 6, we compare the subthreshold scaling of p L with distance for the different EC schedules. We choose three sets of error parameters, each below threshold for all three EC schedules, that illustrate when each schedule should be used. For noise with very large erasure bias, Fig. 6(a) illustrates that using four ECs per round suppresses logical error rates the best. When the noise is biased toward Pauli errors, ECs are less useful, so the 1 EC schedule with the free EC is optimal (Fig. 6(b)). For a realistic scenario with a 10 × erasure bias and nonzero measurement error, 2 EC gives the best scaling (Fig. 6(c)).

## IV. SUPERCONDUCTING ARCHITECTURES

A promising realization of an erasure qubit is via the dual-rail qubit [9, 14]. The dual-rail qubit converts energy relaxation, which is a fundamental noise mechanism in superconducting architectures, to erasures. It also exhibits robustness to the dephasing noise of the constituent transmons. We now discuss how the physical properties of the dual-rail qubit translate to its error-correcting abilities.

## A. Dual-rail qubits

A dual-rail qubit is an erasure qubit composed of two coupled transmons. The computational subspace is spanned by basis states

$$\left | \bar { b } \right \rangle = \frac { 1 } { \sqrt { 2 } } \left ( \left | 0 1 \right \rangle + ( - 1 ) ^ { b } \left | 1 0 \right \rangle \right ), \ \ b = 0, 1 \,. \quad ( 2 ) \quad \ q u b i _ { \substack { \ q u b i \\ \text{subs} } }$$

The dominant noise mechanism for superconducting qubits is amplitude damping, where the excited state | 1 ⟩ of one of the component transmons relaxes to the ground state | 0 . ⟩ This error removes the dualrail qubit from the computational subspace, and the resulting orthogonal state | 00 ⟩ can be subsequently detected.

The dual-rail qubit is made up of two coupled resonant transmons. In the rotating wave approximation, they are described by the Hamiltonian

$$\overset { \text{av us} } { \text{rates} } \quad H _ { 0 } = \sum _ { i = 1, 2 } \left ( \omega _ { \text{DR} } a _ { i } ^ { \dagger } a _ { i } + \frac { \eta } { 2 } a _ { i } ^ { \dagger } a _ { i } ^ { \dagger } a _ { i } a _ { i } a _ { i } \right ) + g _ { 1 2 } \left ( a _ { 1 } ^ { \dagger } a _ { 2 } + \text{h.c.} \right ).$$

(3)

Here, i = 1 2 labels the transmons, , a i , a † i are the ladder operators, ω DR and η &lt; 0 are the transmon frequency and nonlinearity, respectively, and g 12 is the coupling frequency. The computational basis states ∣ ∣ ¯ b 〉 are the eigenstates of the dual-rail Hamiltonian with a single excitation and energies E | ¯ b ⟩ = ω DR +( -1) b g 12 .

The amplitude damping channel describes the process by which an excitation in the transmon escapes to the cold environment. It can be modeled by the jump operators L i = √ 1 /T ( ) i 1 a i , where T ( ) i 1 denotes the energy relaxation time of transmon i = 1 2. , Currently, typical energy relaxation rates for a high-quality transmon are T 1 ≈ 100 µ s [25]. Within the computational subspace of the dual-rail qubit, the dephasing time is much slower than that of the constituent transmons and can be as high as T (DR) ϕ = 1ms since the dual-rail internal coupling behaves analogously to dynamical decoupling [9], as was demonstrated in Ref. [14]. ECs on dual-rail qubits can be implemented as described in Ref. [9].

One straightforward method to restore a dual-rail qubit from the erased state to the computational subspace is by using an XY drive, which we call a unitary pulse. One can drive the qubit from the | 00 ⟩ state to a particular state within the computational subspace, e.g., the state | ¯ . 0 ⟩ As the XY drive is a unitary operation limited to the transmon subspace, it will inevitably transform a computational state

FIG. 6. Scaling of the logical error rate with the distance d for noise parameters e , p and q in the correctable region. (a) For pure erasure noise, the 4 EC schedule performs the best. (b) For pure Pauli noise, ECs introduce unnecessary noise, so the 1 EC schedule is the best. (c) For erasure biased noise, the 2 EC schedule is the best.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-faacda3a-61f0-4cd9-b7e0-2cad26bce8e7.jpg)

to the erased state | 00 ⟩ 2 . As a result, this process should only be performed conditioned on detecting erasures, with false-positive erasure detection introducing additional erasures.

where | B + ⟩ is the closest state in the double excitation subspace to the | 11 ⟩ state (see Appendix E). The transition frequency between the two states is protected by the shift

An alternative method of implementing the reset operation is by removing the uncertainty in the erasure status of the dual-rail qubit to a different degree of freedom, such as the readout resonator, and applying a single pulse that maps the erased state back to the computational subspace without dephasing the computational subspace. We refer to such pulse as a one-way pulse. Similar pulses were used for conditional reset and leakage reduction using the transmon's readout resonator [26-29].

Labeling the total state as | dual-rail ⟩⊗| resonator , ⟩ we can implement a pulse that transforms | 00 ⟩ ⊗ | 0 ⟩ ↦→ | ¯ 0 ⟩⊗| ⟩ 1 without disturbing the computational subspace. One way to achieve this is to use a highfrequency parametric drive. If the readout resonator is coupled symmetrically to the dual-rail qubit, the combined Hamiltonian is

$$H & = H _ { 0 } + \delta \cos \left ( \omega _ { d } t \right ) \sum a _ { i } ^ { \dagger } a _ { i } & \text{em} _ { \substack { \text{the $\mathbb{C}$} \\ \text{the $\mathbb{C}$} } } \\ & + \omega _ { R O } r ^ { \dagger } r + g _ { R O } \left ( r ^ { \dagger } + r \right ) \sum _ { i = 1, 2 } \left ( a _ { i } ^ { \dagger } + a _ { i } \right ) \,, \quad \left ( 4 \right ) & \text{shift} \, \cdot \\ & \text{dual} \, \text{not} \, a \, \text{$$

where r , r † are the ladder operators of the readout resonator, ω RO is its frequency, and g RO is its coupling to the transmons. The frequency of both transmons are modulated with amplitude δ and frequency ω d .

↦

When ω d = ω DR + ω RO + g 12 , the transition | 00 ⟩ ⊗ | 0 ⟩ → | ⟩ ⊗ | ¯ 0 1 ⟩ becomes resonant with Rabi frequency g RO δ/ √ 2 ω d . The next closest resonance is the transition from | ¯ 1 ⟩⊗| ⟩ 0 to the state | B + ⟩⊗| ⟩ 1 ,

2 It is possible to design a pulse that transitions a computational state to a different non-computational state, but that would be even more detrimental.

$$\Delta = \frac { \eta } { 2 } + \sqrt { \frac { \eta ^ { 2 } } { 4 } + 4 g _ { 1 2 } ^ { 2 } } \approx \frac { 4 g _ { 1 2 } ^ { 2 } } { | \eta | } \,, \quad \ \ ( 5 )$$

with the approximation valid for g 12 ≪| | η .

Standard transmons and cavities would require a very high frequency ( ≳ 10 GHz) drive. Parametrically driving the transmons at the min/max point of the tuning curve would reduce this frequency by a factor of 2, using the quadratic response to the flux, at the expense of reducing the speed of the pulse.

It is also possible to implement the reset operation for the dual-rail qubit using quantum operations that preserve the erasure information. For example, similar to Ref. [14], the erasure detection can be implemented using an ancilla qubit. The ancilla is coupled to the dual-rail qubit, shifting its energy gap when the dual-rail qubit decays outside the computational subspace. An XY drive at the shifted frequency excites the ancilla only when the dual-rail qubit is erased, without dephasing the computational subspace. A reset pulse can be applied to the dual-rail qubit by applying a cross-resonance drive on the ancilla [30], which would only excite the dual-rail qubit when the ancilla is excited. The erasure information can later be recovered by measuring the ancilla.

## B. Comparison with conventional transmons

In Sec. III, we characterized the performance of different EC schedules in terms of error parameters e , p , and q . These general results can be used to draw conclusions about specific physical scenarios without obtaining additional data from simulations.

To do this, we express the results in terms of physical parameters. For superconducting qubits, the relevant quantities are the amplitude damping T 1 and dephasing T ϕ times, as well as the operation times T 2 Q for 2Q gates and T M for readout or ECs with reset. The goal in this section is to identify the conditions under which the performance of the dual-rail qubit exceeds that of the transmon and when it is advantageous to perform more frequent ECs.

In the dual-rail qubit, errors mainly arise from amplitude damping (causing erasures) and dephasing (causing Pauli errors). After time t ≪ T 1 , the probability of erasure is approximately e = t/T 1 . Averaging over the syndrome extraction circuit, for the l EC schedule, where l = 1 2 4, we can estimate , , the effective erasure rate per operation as

$$e = \frac { 4 T _ { 2 Q } + l T _ { M } } { ( 4 + l ) \, T _ { 1 } }.$$

The probability of Pauli error is determined by the dephasing rate T ϕ . The parameter p represents the probability of a gate failure, which is doubled when the qubit participates in a 2Q gate. Averaging over the syndrome extraction circuit again, we obtain

$$p = \frac { 4 T _ { 2 Q } + l T _ { M } } { ( 2 + l ) 2 T _ { \phi } } \,, \quad \quad ( 7 ) \quad \text{FI} \\ \text{th}$$

where the factor of 2 + l in the denominator is due to each 2Q gate contributing half an error location on the qubit.

When considering a single transmon, we consider an optimistic situation with no leakage. Hence, we approximate the scenario with e = 0 and

$$p = \frac { 4 T _ { 2 Q } + T _ { M } } { ( 2 + 1 ) 2 T _ { \phi } } + \frac { 4 T _ { 2 Q } + T _ { M } } { ( 2 + 1 ) \, 2 T _ { 1 } } \,, \quad \ ( 8 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where we take l = 1 in Eq. (7) and add an additional noise mechanism with rate 1 / (2 T 1 ) modeling the qubit energy relaxation. 3 We motivate these approximations in Appendix D.

We now express the error parameters e , p and q that describe the threshold surfaces of the correctable regions in Fig. 1 in terms of the physical parameters T 1 , T ϕ , T 2 Q and T M of the transmon. For simplicity, we assume q = 0 and use the following ansatz

$$e / e ^ { * } + p / p ^ { * } = 1 \,, \quad \quad ( 9 ) \quad \text{ex}$$

3 The factor of 2 can be thought of as the transmon occupying the ground state half the time, while the dual-rail qubit is always excited, so it decays faster.

FIG. 7. Cross sections of the threshold surfaces in the (a) ( T 2 Q , T M ) and (b) ( T , T 1 ϕ ) planes for different EC schedules. In (a), we set T 1 = 100 µ s for both the dual-rail qubit and the transmon. The dual-rail qubit dephasing time is extended to T (DR) ϕ = 1ms due to its internal coupling. For the transmon, we consider two dephasing times T ϕ = 100 µ , s ∞ . The crossover regime between the 4 EC, 2 EC and 1 EC schedules occurs around T M = T 1 × e ∗ 4EC . In (b), we set T 2 Q = T M = 150ns. The extended dephasing time of the dual-rail qubit is accounted for by setting T ( DR ) ϕ = 30 T ϕ . To compare it with the transmon, we choose two regimes that approximate the gate times in Refs. [31, 32]. Transmons with long T ϕ are outperformed by 30-40% in terms of minimal T 1 , depending on how fast they implement 2Q gates.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-443a4d2a-7719-46aa-9ece-9ca725174591.jpg)

where e ∗ and p ∗ are the erasure and Pauli thresholds, to interpolate the cross section of the ( e, p, q ) threshold surface. We present the results in Fig. 7.

When we fix T , T 1 ϕ , the dual-rail qubit can tolerate slower gates compared to the transmon. Fig. 7(a) can be interpreted as follows. When the ratio T ϕ T 1 exceeds e ∗ p ∗ , erasures become dominant. The dualrail qubit threshold curves are essentially those described by Eq. (6), with a threshold corresponding to e ∗ (which depends on the EC schedule). For the transmon curve where T ϕ = ∞ , it aligns with Eq. 7 with p ∗ = 0 007 for a standard surface code. . The transmon curve with T ϕ = T 1 incorporates

double the Pauli error. In the limit of T ϕ → ∞ (Fig. 7(a)), the maximal measurement time for the dual-rail qubit is similar to that of the transmon. We see a transition regime where less detections become preferable to 4 EC, approximately at the point T M = T 1 × e ∗ 4EC , where e ∗ 4EC is the erasure threshold value for 4 EC.

One would expect the dual-rail qubit to be advantageous compared to the transmon with the same T 1 . However, since the decay rate of the dual-rail qubit is twice the decay rate of the transmon and dualrail gates tend to be slower, this benefit will only be significant if measurements for the dual-rail qubit can be realized on a par with measurements for the transmon. At that point, the elevated threshold of the dual-rail qubit substantially reduces the minimal necessary T 1 to approximately half that of the transmon.

When we consider the the 2Q gate and measurement times T 2 Q and T M , dual-rail qubits greatly outperform transmons with bad T ϕ (Fig. 7(b)). However, even with T ϕ → ∞ , as long as the dual-rail qubit can be measured at a similar rate, e.g., in 150 ns, then the additional erasure information results in the minimal necessary T 1 reduced by 2 × .

For transmons, there is a trade-off between tunability and gate time. Tunable transmons admit fast gates but worst coherence. Fixed transmons have much better coherence at the expense of slower gates. The dual-rail qubit has the advantage of both. Namely, it has a better threshold when gates are implemented slowly, and it can be tunable with a very long lifetime. In Fig. 7(b), we plot the parameters for two devices that roughly characterize the stateof-the-art for tunable [31] and fixed [32] transmons. The low T ϕ of tunable transmons significantly inhibits their fidelity, greatly reducing the advantage of their fast gates. Given that T 1 ≳ T ϕ for fixed transmons, the dual-rail qubit's robustness to dephasing offers a significant advantage.

More generally, state-of-the-art decoherence rates compared with gate and measurement times would suggest that both transmons and dual-rail qubits surpass the QEC threshold, a result that for transmons was never achieved in experiments [31, 33]. For transmons, this discrepancy is explained by both leakage to higher excited states and the difficulty to consistently fabricate many transmons with long coherence times on a single device. Dual-rail qubits are more robust to fabrication issues since they have a wide tunability range with long coherence [14]. However, as demonstrated by the transmon scenario, coherence is not always the primary factor affecting fidelity, which could diminish the benefits of the dualrail qubit.

Dual-rail qubits are more resistant to leakage dur- ing 2Q gates. The energy gap in the computational subspace is approximately 20 × smaller than the gap to the next excited state, highly suppressing any leakage due to driven singleand 2Q gates. Nondriven 2Q gates can also leak the dual-rail qubit by moving an excitation to its partner dual-rail qubit. This type of leakage is protected only by the detuning between the two dual-rail qubits and the transmon nonlinearity, but it erases the partner dual-rail qubit, making it easier to identify. It is also reasonable to assume that leakage detection can be integrated into the erasure detection scheme. Moreover, inherent robustness of the dual-rail qubit to dephasing opens the possibility to utilize transmons versions with less coherence or other qubits.

## C. Imperfect erasure reset

Next, we investigate how false-negative and falsepositive erasure detections affects the performance of dual-rail qubits. We also compare the two implementations of reset operations and find their effect on the QEC threshold. We find that using one-way pulses, and only implementing them following erasure detection, reduces the effect of false-positive and false-negative erasure detection by roughly a half. To simplify the analysis, we simulated noise that only occurs during gates [9] and expressed the results as the averaged erasure per operation. This modification slightly shifts the value of e ∗ , but does not change the qualitative results.

The process of restoring the dual-rail qubit may need to incorporate the classical feedback and operate on a time scale similar to that of other quantum operations. Here, we show that the reset operation of the dual-rail qubit that does not depend on the erasure information amplifies the noise from erroneous erasure detections. This type of classical feedback has been demonstrated in superconducting qubits [34, 35], with an extra delay that, in principle, can go below 50 ns. Single Flux Quantum (SFQ) devices [36] have the potential to further decrease feedback latency to a minimal level.

Many current devices still suffer from a considerable feedback delay that prolongs the measurement time. The ancilla based recovery pulse that preserves the erasure information can allow some additional leeway in designing the cycle of current devices. This can be combined with a decreased frequency of ECs, reducing the effect of the measurement time on the erasure rate in Eq. (6) (albeit also lowering the QEC threshold).

In systems with a high rate of false-negative erasure detections, one might consider using the oneway pulse on all of the qubits following every EC,

independently of the EC outcome. Erasure qubits that were missed due to a false-negative detection would return to the computational subspace without raising any erasure flag. This approach would result in the loss of information regarding the qubit decay but would prevent erasures from spreading and harming additional qubits. However, we find through numerical simulations that the information about erasures is more valuable, and pulsing all the qubits effectively reduces the QEC threshold.

Assuming no Pauli errors, i.e., p = 0, we find the maximal erasure rate e for a range of false-negative erasure detection rates q fn ; see Fig. 8(a). We note that the data can be fitted with the following ansatz

$$\frac { e } { e ^ { * } } + \frac { q _ { \tt f n } } { q _ { \tt f n } ^ { * } } = 1, \quad \quad \quad ( 1 0 )$$

where q ∗ fn is a fitting parameter. When the reset pulse is applied selectively to the qubits once they are identified as erased, we obtain q ∗ fn ≈ 0 49 0 33 0 28 for the 4 EC, 2 EC and 1EC sched-. , . , . ules, respectively. When the reset pulse is applied to all qubits following every EC, we find q ∗ fn ≈ 0 3 . , 0 17 0 14. . , .

We implemented numerical simulations to assess the impact of false-positive erasure detection on the logical error rate. When a false-positive detection occurs with a one-way pulse, the qubit is not erased in the simulation, but the decoder is (incorrectly) informed about an erasure. The decoder knows the probability of false-positive detections and adjusts the probabilities on the matching graph accordingly.

To simulate the unitary pulse, following each erasure detection that failed with a false positive error, we either erase the qubit or completely depolarize it, each with probability 1 / 2. This captures our assumption that no coherences are created between the computational and erasure subspaces.

The rate of false-positive erasure detections would likely be comparable to the 2Q gate erasure rate. Since the main advantage of dual-rail qubits comes from the knowledge about the locations of erasures, the false-positive detection rate would have a noticeable effect on the QEC thresholds, as a large fraction of erasure detections would be inaccurate. This effect would be more pronounced in the 4 EC schedule compared to the 2 EC and 1 EC schedules.

Assuming no Pauli errors, i.e., p = 0, we find the maximal erasure rate e for a range of false-positive erasure detection rates q fn ; see Fig. 8(b). The data can be fitted with the following ansats

$$\frac { e } { e ^ { * } } + \frac { q _ { \text{fp} } } { q _ { \text{fp} } ^ { * } } = 1, \quad \quad \ \ ( 1 1 ) \quad \text{of} \\ \text{no}.$$

where q fp is a fitting parameter. We find that for the one-way pulse q ∗ fp ≈ 0 043 0 059 0 074 for the 4 EC, . , . , .

FIG. 8. Cross sections of the threshold surfaces in the (a) ( e, q fn ) and (b) ( e, q fp ) planes, where q fn and q fn are the rates of false-positive and false-negative erasure detections, respectively. In (a), we use either a selective reset pulse on qubits that were identified as erased, or apply a reset pulse to all qubits. The trend lines are fitted using Eq. (10). In (b), we use the one-way or unitary pulse to reset erased qubits. The trend lines fitted using Eq. (11).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-d78f17c9-98af-4a80-a23a-1993da5e5b58.jpg)

2 EC and 1 EC schedules, respectively. In contrast, for the unitary pulse q ∗ fp ≈ 0 02 0 026 0 03. . , . , .

## D. Biased erasure on dual-rail qubits

In our simulations, we incorporate the lack of knowledge about the timing of the jump operator during the 2Q gate by fully depolarizing the partner dual-rail qubit. This approach is justified since the interaction term is not oriented along the Z or X axes, being a flip-flop term X X 1 2 + Y Y 1 2 .

A possible alternative is to use the level repulsion of two dual-rail qubits connected with the same capacitive coupling, similar to Ref. [37]. If the two dual-rail qubits are not resonant, the interaction would be along the ZZ direction, evolving according

to a Hamiltonian of the form

$$H = \omega _ { 1 } Z _ { 1 } + \omega _ { 2 } Z _ { 2 } + g _ { z z } Z _ { 1 } Z _ { 2 } \,, \quad ( 1 2 ) \quad,$$

where Z i is the Pauli Z operator for the dual-rail qubit i = 1 2. , The interaction stops when one of the qubits decays to the | 00 ⟩ state. Without loss of generality, let it be the first qubit. Then, when averaged over all possible decay times, following a reset channel on the decayed qubit, the noise channel acts on the state ρ in the computational subspace of both qubits as follows

$$\mathcal { R } \left ( \rho \right ) = \mathcal { D } _ { 1 } \left [ \frac { 1 } { \pi } \int _ { 0 } ^ { \pi } e ^ { i Z _ { 1 } Z _ { 2 } \theta } \rho e ^ { - i Z _ { 1 } Z _ { 2 } \theta } \right ] \quad ( 1 3 )$$

$$= \mathcal { D } _ { 1 } \left [ \frac { 1 } { 2 } \rho + \frac { 1 } { 2 } Z _ { 1 } Z _ { 2 } \rho Z _ { 1 } Z _ { 2 } \right ] \quad \ \ ( 1 4 ) \, \stackrel { \ddots } { \text{cod} } \\ \text{for} \, \, \,$$

$$= \frac { 1 } { 2 } \mathcal { D } _ { 1 } \left [ \rho \right ] + \frac { 1 } { 2 } Z _ { 2 } \mathcal { D } _ { 1 } \left [ \rho \right ] Z _ { 2 } \,. \quad \quad ( 1 5 )$$

Here, D 1 is a channel that fully depolarizes the first qubit and acts trivially on the second qubit, modeling the decay and reset of the decayed qubit. Consequently, only the phase of the qubit that did not decay is lost, which effectively biases the noise.

Gates based on level repulsion are available in many quantum computing architectures and would always give the benefit mentioned above. In systems in which this type of gate is not available or is too slow, a biased gate is still possible through a coherent excitation of a state outside the computational subspace, but such implementation requires a diagonal term of the form Z Z 1 2 . In neutral atoms, for example, this comes from the Rydberg blockade term ∆ | rr ⟩ ⟨ rr | . We analyze this case in further details in Appendix C, similar to the case in Ref. [38].

Note the difference between this approach and Ref. [38], where the source of bias partially comes from the decay mechanism. Since only the Rydberg state decays, an erasure event effectively measures the qubit in the Rydberg state, losing only its phase. The effect on the unerased qubit depends on the specific implementation of the gate and can result in unbiased noise and leakage (see Appendix C).

We consider the XZZX surface code [39-42], which can exploit bias in the underlying noise. In Ref. [40], a syndrome extraction circuit is designed with two CZ gates and two CX gates. When implemented using CZ and Hadamards, two of the possible eight Z errors coming from the partner-qubit decaying are therefore converted to X errors. We use the gate schedule described in Ref. [31], which suppresses any hook errors. The improved threshold shows that the XZZX surface code can exploit the bias introduced by the ZZ gate; see Fig. 9.

FIG. 9. Logical error rate p L for the XZZX surface code as a function of the erasure rate e with d rounds of QEC. We find threshold values e ∗ ≈ 0 0172 0 0175 0 014 . , . , . for the 4 EC, 2 EC, and 1 EC schedules, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-19033b68-3a97-46b4-87c0-5bee9ba8a2ed.jpg)

## V. DISCUSSION

In this article, we studied the performance of the surface code with erasure qubits. We probed the threshold surfaces for different EC schedules and analyzed their subthreshold behavior. We identified the parameter regimes in which it is optimal to use each EC schedule. Notably, even with less precise erasure information from infrequent ECs, the surface code with erasure qubits can outperform the standard scheme. In terms of hardware realizations of erasure qubits, we analysed the dual-rail qubit and compared it with the standard transmon-based architecture. While the dual-rail qubit offers an advantage in suppressing dephasing errors, its effectiveness in tackling amplitude damping errors (assuming state-of-the-art parameters) requires gates and measurements which are faster than about 200 ns.

Our study is a distinctive approach to demonstrating QEC in the lab. Although fabricated superconducting devices are usually created for a particular QEC code, the QEC protocol can be modified based on the observed qubit lifetime, gate fidelities, and, in particular, measurement efficiencies. For example, the achievable measurement time, following an experimental optimization of the readout procedure, should dictate the rate of ECs within the circuit. A similar result is applicable to cavity-based dual-rail qubits, as we discuss in Appendix F 1. For erasure qubits based on neutral atoms, erasure detection is cheap enough to make the 4 EC schedule the most beneficial; see Appendix F 2.

While we focused on the surface code as a quantum memory, we expect other QEC protocols to similarly benefit from erasure qubits. Future study

could explore other QEC codes, such as quantum low-density parity-check codes [43], or protocols for implementing logical operations. It would be interesting to see if the improvement from erasure qubits is universal or if certain QEC protocols particularly benefit from erasure information.

The decoding problem with erasure qubits is another direction of future explorations. Although our work and Ref. [18] provide ways to decode general erasure circuits by converting them to stabilizer circuits, the loss in performance due to the invoked approximations is unclear. Decoders that directly work with erasure circuits, perhaps for specific QEC codes, may have more rigorous performance guarantees and use erasure information optimally.

Further optimization of QEC protocols with erasure qubits are possible. One may consider independently placing ECs and reset operations throughout the circuit. Performing multiple ECs before resetting could be useful in the presence of high falsepositive erasure detection rates in order to gain more confidence in the erasure information, or if ECs can be performed concurrently with gates so that only reset operations incur noise. Furthermore, circuits may be implemented adaptively, e.g., once an erasure is detected, we may stop extracting the syndrome of stabilizers involving that qubit to mitigate the spread of errors. All these approaches are complementary, and incorporating them would further increase the benefits of using erasure qubits.

## ACKNOWLEDGMENTS

We thank F. Brand˜ ao, A. Haim, J. Iverson, H. Levine and C. Ryan for insightful discussions. We thank A. Grimsmo for helping correctly implement the XZZX surface code. S.G. acknowledges funding from the Air Force Office of Scientific Research (FA9550-19-1-0360).

Note added.We would like to bring the reader's attention to a related, independent work by Chang et al. [44], which explores the effects of physicallymotivated imperfect erasure checks on the performance of the surface code and appeared in the same arXiv posting.

## Appendix A: Example of circuit conversion for decoding

We present an explicit example of the conversion of an erasure circuit to a stabilizer circuit, highlighting the difference between the exact conversion of Ref. [18] and that the approximate conversion of Sec. III A. Consider a segment of an erasure cir- cuit with three erasure locations and two entangling gates, as in Fig. 10(a). Assume that each erasure location has erasure rate e and the ECs have false positive and negative detection rates q . The EC outcomes have probabilities

$$\underset { \omega } { \overset { \int } { n - 1 } } \quad \Pr ( E C = 1 ) = [ 1 - ( 1 - e ) ^ { 3 } ] ( 1 - q ) + ( 1 - e ) ^ { 3 } q \,, \\ \underset { \omega } { \overset { \int } { n - 1 } } \,.$$

$$\ r a \quad \Pr ( E C = 0 ) = 1 - \Pr ( E C = 1 ) \,. \quad \quad \ ( A 2 )$$

If a positive detection occurs, the posterior probability of the erasure happening before gate G i is

$$\bar { a } _ { i } = \frac { ( 1 - ( 1 - e ) ^ { i } ) ( 1 - q ) } { \Pr ( E C = 1 ) } \quad \quad ( A 3 )$$

for i ∈ { 1 2 3 , , } (counting G 3 as the second reset operation). If EC = 0, the probability is

$$\bar { a } _ { i } = \frac { ( 1 - ( 1 - e ) ^ { i } ) q } { \Pr ( E C = 0 ) } \,. \quad \quad ( \mathbf A 4 )$$

Thus, in the approximate scheme, we place depolarizing channels with error probability 3¯ a / i 4 at each location F i in the circuit in Fig. 10(b). Equivalently, we place Pauli X , Y , and Z error channels with strength p i = 1 2 ( 1 -√ 1 -¯ a i ) at each location.

FIG. 10. (a) A segment of an erasure circuit with two entangling gates and three erasure locations. (b) A converted stabilizer circuit. Locations F 1 , F 2 , F 3 are replaced with correlated depolarizing channels to obtain an equivalent stabilizer circuit or uncorrelated depolarizing channels in the approximate conversion scheme.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-cf529f4f-6b77-409b-a51d-bfbd9805dd64.jpg)

In contrast, exact conversion requires capturing the correlated errors induced by erasures. If an erasure occurs at location E i , it would result in depolarizing at locations F j for all j ≥ i . In particular, the possibility of erasure at E 1 requires us to insert correlated depolarizing channels at locations F 1 , F 2 , F 3 with some probability. Decomposing this error into binary random variables requires introducing 63 error mechanisms, which correspond to the nontrivial Pauli strings of length three. For more details, see Ref. [18].

## Appendix B: Numerical simulation details

In this appendix, we provide some more details on our numerical simulations and present several supplementary plots.

In our simulations, we noiselessly initialize an eigenstate of a chosen logical Pauli operator of the surface code with distance d . Then, we run 3 d rounds of the noisy syndrome extraction circuit. Finally, we perform an ideal measurement of the logical operator. A logical error is reported if the final decoded value of the logical operator is different than when initialized. The stabilizer simulations are done using Stim [23]. In the decoding process, we convert the erasure circuit to a stabilizer circuit, use Stim to obtain a decoding hypergraph and further decompose the hyperedges by approximating them as independent edges, and find the minimum-weight perfect matching solution using PyMatching [45].

In the sampling process, we first sample the erasures and erasure detection events. For each erasure sample, we then sample the Pauli errors in the circuit. The reason for doing the sampling in two stages is that each erasure sample corresponds to a different Stim circuit. As the overhead of initializing a Stim circuit is significant compared to sampling from the circuit and decoding, it is efficient to sample from the same Stim circuit multiple times. The logical error rate p ′ L is the total fraction of errors obtained over at least 1000 erasure samples and 200 samples of Pauli errors per corresponding circuit. We then report the normalized logical error rate per d syndrome extraction rounds, calculated as p L = 1 2 (1 -(1 -2 p ′ L ) 1 / 3 ) ≈ p ′ L / 3.

Threshold points are calculated by sweeping a single error parameters y ∈ { e, p, q } near the suspected point in ( e, p, q ) phase space for various code distances d . We fit the universal scaling ansatz for critical points of phase transitions to the data,

$$p _ { L } = a x ^ { 2 } + b x + c \,, \quad \quad ( B 1 ) \quad \text{th}$$

where

$$x = ( y - y ^ { * } ) d ^ { \alpha } \,. \quad \quad \text{(B2)}$$

Here, a, b, c, y ∗ , α are fitting parameters, and y ∗ is the estimated threshold. For an example calculation, see Fig. 11.

Fig. 12 shows the fit of the Eq. (1) ansatz for the subthreshold logical error rate. The data presented are the ones used to obtain the α values in Fig. 5(c).

## Appendix C: Biased-noise gates

In this appendix, we consider two different gates on erasure qubits that have the biased noise mecha-

(a)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-b8b3e315-664a-422b-ac77-5c1079b99751.jpg)

(b)

FIG. 11. Sample calculation of a threshold point. (a) We test values of p near a suspected threshold point for fixed e and q . (b) Rescaled data using the universal ansatz in Eq. (B1), giving a threshold of p = 0 00294. .

nism where the erased qubit decays to an identifiable erased state, and its partner qubit for the 2Q gate is only affected by noise in the Z direction.

As mentioned in the main text, in the case of level repulsion between two transmons or dual-rail qubits, the 2Q gate stops at an unknown time, which results in an unknown phase on the partner qubit. In this case the error on the unerased qubit is biased to a single axis, the Z axis.

A second type of gate with this property is described in Ref. [46], which is implemented on neutral atoms. Assume the two atoms begin in the state

$$\text{the} \quad \left | \psi \right \rangle = a _ { 0 0 } \left | 0 0 \right \rangle + a _ { 0 1 } \left | 0 1 \right \rangle + a _ { 1 0 } \left | 1 0 \right \rangle + a _ { 1 1 } \left | 1 1 \right \rangle \,, \, \left ( C 1 \right )$$

where | 0 ⟩ , | 1 ⟩ are the computational states of the atom. We denote the excited Rydberg state as | r ⟩ . We pulse both qubits with the following Hamiltonian:

$$H = \Omega \left ( \mathbb { I } \otimes \left | 1 \right \rangle \left \langle r \right | + \left | 1 \right \rangle \left \langle r \right | \otimes \mathbb { I } \right ) + \Delta \left | r r \right \rangle \left \langle r r \right |, \text{ } \left ( C 2 \right )$$

which is written in the interaction picture with re-

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-763e1f6c-c7da-4cda-bae8-aebc8ea92ae7.jpg)

FIG. 12. Fitting the subthreshold logical error rate scaling ansatz of Eq. (1) to obtain the α values presented in Fig. 5(c). In each plot, the different colors correspond to different values of x .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-ce7dbac2-21fe-4c8d-9b0e-b496b8eb1274.jpg)

spect to the gaps:

$$H _ { 0 } = \omega _ { 0 } \left | 0 \right \rangle \left \langle 0 \right | + \omega _ { 1 } \left | 1 \right \rangle \left \langle 1 \right | + \omega _ { r } \left | r \right \rangle \left \langle r \right | \,. \quad ( C 3 ) \quad \text{$\emph{subarray}$} \,.$$

In the regime of strong Rydberg blockade (∆ ≫ Ω), the state evolves in the interaction picture with respect to H 0 according to where | W ⟩ = [ 1 | r ⟩ + | r 1 ] ⟩ / √ 2. Now assume that the first atom decayed as in Ref. [38] from its Rydberg state | r ⟩ to an identifiable state | e ⟩ . The unnormalized state of the two atoms becomes

$$| \cdot$$

$$\langle r \lhd \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \\ | \psi \left ( t \right ) \rangle & = a _ { 0 0 } \left | 0 0 \rangle + a _ { 0 1 } \left [ \cos \left ( \Omega t \right ) | 0 1 \rangle + i \sin \left ( \Omega t \right ) | 0 r \rangle \right ] \\ & + a _ { 1 0 } \left [ \cos \left ( \Omega t \right ) | 1 0 \rangle - i \sin \left ( \Omega t \right ) | r 0 \rangle \right ] \\ & + a _ { 1 1 } \left [ \cos \left ( \sqrt { 2 } \Omega t \right ) | 1 1 \rangle - i \sin \left ( \sqrt { 2 } \Omega t \right ) | W \rangle \right ] \,, \\ & \quad ( C 4 ) \quad \text{occur}$$

$$\left | 0 r \right \rangle \right ] \\ \left | \psi \left ( t \right ) \right \rangle = a _ { 1 0 } \sin \left ( \Omega t \right ) \left | e 0 \right \rangle + a _ { 1 1 } \sin \left ( \sqrt { 2 } \Omega t \right ) \left | e 1 \right \rangle / \sqrt { 2 } \,, \\ \mathcal { V } \right \rangle \right ] \,, \quad \underset { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { \dots } { \dots } \, \sum \, \frac { 1 } { \dots } \,.$$

occurring at some random time t . Continuing the interaction for some additional time t 2 , the state be-

comes

$$| \psi \left ( t + t _ { 2 } \right ) \rangle & = - i a _ { 1 0 } \sin \left ( \Omega t \right ) | e 0 \rangle \quad \begin{pmatrix} \text{the t} \\ \text{nels} \end{pmatrix} \\ & - i a _ { 1 1 } \sin \left ( \sqrt { 2 } \Omega t \right ) \cos \left ( \Omega t _ { 2 } \right ) | e 1 \rangle \, / \sqrt { 2 } \quad \begin{pmatrix} \text{scrib} \\ \text{the r} \\ \text{for} \\ \text{$-} a _ { 1 1 } \sin \left ( \sqrt { 2 } \Omega t \right ) \sin \left ( \Omega t _ { 2 } \right ) | e r \rangle \, / \sqrt { 2 }. \end{pmatrix} \\ & \quad \begin{pmatrix} \text{To} \\ \text{syste} \\ \text{trans} \\ \text{Hern} \end{pmatrix}$$

At the end of the gate, either the undecayed qubit stays in the computational manifold span {| 0 ⟩ , | 1 ⟩} or in the Rydberg | r ⟩ state. In the computational manifold, notice that its transformation to the new state commutes with the Z operator acting on the second atom. Hence, the undecayed atom only lost its phase, as any stabilizer with a Z support on the second atom would not be affected by this transformation. The amplitude of the | er ⟩ state at the end of the gate describes a leakage mechanism, which does not fit into our formalism, and might present an issue for implementing biased-noise gates on neutral atoms.

## Appendix D: Error channels on dual-rail qubits and transmons

The large benefit of detecting erasure compared with decoding Pauli noise is slightly obscured by the way we define the Pauli and erasure rates. The erasure probability is per-operation-per-qubit, while the Pauli rate is only per-operation, which for 2Q gates, would introduce an additional factor of 2. To put it differently, a 2Q gate is twice as likely to fail compared with a 1Q operation of the same length, but we give both the same probability (which is common in the QEC literature).

Therefore, we need to be careful when averaging the failure probability for Pauli and erasure noises. For Pauli noise, we take two qubits and consider the probability of either of them failing in an imaginary syndrome extraction cycle in which we apply a CNOT between them four times. This cycle will include 4 CNOTs and 2 l measurements. The probability of either of them failing during that time is 2 × (4 T 2 Q + lT M ) /T 1 , and the average per operation, i.e. the Pauli error probability, is then

$$p = \frac { 8 T _ { 2 Q } + 2 l T _ { M } } { ( 4 + 2 l ) \, T _ { 1 } } = \frac { 4 T _ { 2 Q } + l T _ { M } } { ( 2 + l ) \, T _ { 1 } } \,. \quad ( D 1 ) \quad \ U$$

For the erasure rate, which is a per-qubit-peroperation rate, we simply average the total cycle time over the number of operations in the cycle:

$$e = \frac { 4 T _ { 2 Q } + l T _ { M } } { ( 4 + l ) \, T _ { M } } \,. \quad \text{(D2)}$$

When comparing dual-rail qubits to transmons, the two systems are affected by different noise channels that are not directly comparable. We now describe the channels and how we compared them in the main text.

For the dual-rail qubit, it is easiest to analyze the system using a non-Hermitian Hamiltonian. Both transmons have some decay time T 1 , and the nonHermitian evolution of the system is given by

$$\text{ubit} \\ \text{nal} \\ \text{new} \\ \text{the} \quad \text{Woc} \ \text{discomlize to the formition not in non} \ \text{$k$} \. \ -$$

We diagonalize the Hermitian part using b ± = ( a 1 ± a 2 ) / √ 2, and then move to the rotating frame with respect to H 0 = Ω 2 [ b † + + b -b † - -b ] :

$$\underset { \text{rel} } { \text{of} } \\ \underset { \text{rel} } { \text{is-} } \quad H _ { I } & = - \frac { i } { 2 } \left ( \frac { 1 } { T _ { 1 } ^ { ( 1 ) } } + \frac { 1 } { T _ { 1 } ^ { ( 2 ) } } \right ) \left [ b _ { + } ^ { \dagger } b _ { + } + b _ { - } ^ { \dagger } b _ { - } \right ] \\ & \quad - \frac { i } { 2 } \left ( \frac { 1 } { T _ { 1 } ^ { ( 1 ) } } - \frac { 1 } { T _ { 1 } ^ { ( 2 ) } } \right ) \left [ b _ { + } ^ { \dagger } b _ { - } e ^ { i \Omega t } + \text{h.c.} \right ]. \quad ( \text{D4} ) \\ \text{ts} \quad \cdot \ \. \ \. \ \. \ \. \ \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \. \ \. \. \ \. \. \ \. \ \. \. \ \. \. \ \. \. \. \ \. \. \ \. \ \. \. \. \. \ \. \. \ \. \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.. \. \. \.. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \ \. \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \ \. \. \. \. \. \.. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.. \. \. \. \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \ \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.. \. \.. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \ \. \. \ \. \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.. \. \. \. \. \. \. \. \. \. \.. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \ \. \. \. \. \. \. \. \. \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.. \. \. \. \. \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \ \. \. \. \. \. \. \.. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.. \. \. \. \. \. \. \ \. \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.. \. \. \. \.. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \ \. \. \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \$$

In the rotating wave approximation (Ω ≫ T ( ) i 1 ), the rotating terms disappear, and both logical states decay to the ground at an average rate

$$T _ { 1 } = \left ( \frac { 1 } { 2 } \left [ \frac { 1 } { T _ { 1 } ^ { ( 1 ) } } + \frac { 1 } { T _ { 1 } ^ { ( 2 ) } } \right ] \right ) ^ { - 1 } \,, \quad \text{(D5)}$$

which is simply T 1 if both transmons have the same decay rate.

To analyze the decay of the transmon, we look at its energy relaxation Linbdladian

$$\cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ } } \, \left ( \mathbf D \, \right )$$

Applied for time t , we have

$$\underset { \begin{array} { c } \text{ur} \\ \text{alg} \\ \text{ply} \\ \text{vill} \\ \text{obs} \end{array} } \quad U _ { T _ { 1 } } \left ( \rho \right ) = \begin{pmatrix} \rho _ { 1 1 } e ^ { - t / T _ { 1 } } & \rho _ { 1 2 } e ^ { - t / 2 T _ { 1 } } \\ \rho _ { 1 2 } ^ { * } e ^ { - t / 2 T _ { 1 } } & \rho _ { 2 2 } + \rho _ { 1 1 } \left ( 1 - e ^ { - t / T _ { 1 } } \right ) \end{pmatrix} \,. \\ \text{others channel cannot have different by an unwindented by}$$

This channel cannot be directly approximated by a Pauli channel, but we can reduce it to a Pauli channel using the Pauli twirling approximation [47]. In this form, the twirled channel is

$$& \text{In this form, the twirled channel is} \\ U _ { T _ { 1 } } ^ { P } \left ( \rho \right ) & = \left ( e ^ { - t / 2 T _ { 1 } } + \frac { 1 } { 2 } e ^ { - t / T _ { 1 } } - \frac { 1 } { 2 } \right ) \rho \\ & + \frac { 1 } { 4 } \left ( 1 - e ^ { - t / T _ { 1 } } \right ) \left [ \frac { 1 } { 2 } \sigma _ { x } \rho \sigma _ { x } + \frac { 1 } { 2 } \sigma _ { y } \rho \sigma _ { y } \right ] \\ & + \left ( \frac { 1 } { 2 } - \frac { 1 - e ^ { - t / T _ { 1 } } } { 4 } - \frac { e ^ { - t / 2 T _ { 1 } } } { 2 } \right ) \sigma _ { z } \rho \sigma _ { z } \,.$$

For shot times t ≪ T 1 , we get

$$U _ { T _ { 1 } } ^ { P } \left ( \rho \right ) & \approx \left ( 1 - t / 2 T _ { 1 } \right ) \rho + \frac { t } { 2 T _ { 1 } } \left [ \frac { 1 } { 2 } \sigma _ { x } \rho \sigma _ { x } + \frac { 1 } { 2 } \sigma _ { y } \rho \sigma _ { y } \right ] \,. \quad \text{ here} \\ T _ { 1 } \shortparallel \colon \in \mathfrak { m l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l b l$$

This is a slightly biased channel along the X,Y direction. In order to approximate this channel as a depolarizing channel of the form

$$U _ { p } \left ( \rho \right ) = \left ( 1 - p \right ) \rho + \frac { p } { 3 } \sum _ { a \in \{ x, y, z \} } \left ( \sigma _ { a } \rho \sigma _ { a } \right ) \,, \quad \text{(D10)} \quad \text{rate} \\ \text{depth} \quad \text{$\dots$}$$

we assign p = t/ 2 T 1 as the probability of any Pauli error.

We model dephasing using the following Markovian Lindblad operator:

$$\mathcal { L } \left ( \rho \right ) = \frac { 1 } { 2 T _ { \varphi } } \left [ \sigma _ { z } \rho \sigma _ { z } - \rho \right ] \,. \quad \left ( D 1 1 \right ) \, \begin{array} { c } \text{bia} \\ \text{the} \\ \end{array} \,.$$

The channel after time t is

$$& U _ { T _ { \varphi }, t } \left ( \rho \right ) = \left ( \frac { 1 } { 2 } + \frac { 1 } { 2 } e ^ { - t / T _ { \varphi } } \right ) \rho + \left ( \frac { 1 } { 2 } - \frac { 1 } { 2 } e ^ { - t / T _ { \varphi } } \right ) \sigma _ { z } \rho \sigma _ { z } \,. \\ & \mathfrak { m } _ { \varphi } \, \text{s.s.s.t.s.t. } \vdots \, \text{$\times T$}$$

For short times t ≪ T φ ,

$$U _ { T _ { \varphi }, t } \left ( \rho \right ) \approx \left ( 1 - \frac { t } { 2 T _ { \varphi } } \right ) \rho + \frac { t } { 2 T _ { \varphi } } \sigma _ { z } \rho \sigma _ { z } \,. \quad ( D 1 3 ) \quad \text{$A$} \\ \text{$\text{$d$} }$$

## Appendix E: One way pulse derivation

The double excitation subspace of the dualrail qubit is spanned by a dark state | D ⟩ = √ 1 2 [ 02 / | ⟩ - | 20 ] and two additional bright states. ⟩ The internal dual-rail qubit coupling couples the state | B ⟩ = √ 1 / 2 [ 02 | ⟩ + 20 ] to the state | ⟩ | S ⟩ = 11 | ⟩ with the following interaction:

$$H = \begin{pmatrix} 2 \omega _ { D R } + \eta & 2 g _ { 1 2 } \\ 2 g _ { 1 2 } & 2 \omega _ { D R } \end{pmatrix} \,. \quad \begin{pmatrix} 0 \end{pmatrix} \,.$$

The eigenstates of this Hamiltonian have energies E ± = 2 ω DR + η 2 ± √ η 2 4 +4 g 2 12 which we label | B ± ⟩ . The transition to the lowest energy state | B + ⟩ is the most resonant during the reset pulse, giving the shift

$$\Delta = E _ { | B _ { + } \rangle \otimes | 1 \rangle } - E _ { | \mathbb { I } \rangle \otimes | 0 \rangle } - \omega _ { d } = \frac { \eta } { 2 } \pm \sqrt { \frac { \eta ^ { 2 } } { 4 } + 4 g _ { 1 2 } ^ { 2 } } \,.$$

## Appendix F: Other platforms

## 1. Cavity dual-rail qubit

When considering a cavity dual-rail qubit, a general assumption is that the limiting factor T C 1 is the cavity amplitude damping, which is dictated by the intrinsic decay and the Purcell effect, while coherence is influenced by the transmon and is fairly small. The erasure rate during a 2Q gate is e = T 2 Q T C 1 = ∆ 2 g 2 1 αT C 1 . The limit on erasure, however, will be dominated by the erasure rate of the transmon e T = T 2 Q T T 1 = ∆ 2 g 2 1 αT T 1 &gt; e C and the dephasing erasure rate on the transmon e T ϕ = T 2 Q T T ϕ = ∆ 2 g 2 1 αT T ϕ &gt; e T . As dephasing does not propagate under dispersive interaction, its main effect would be to reduce the measurement fidelity, and thus, the main erasure bottleneck would be dominated by e T . As it can be assumed that we are in the regime in which T 1 ≪ T ϕ and T ϕ is smaller than the transmon-induced T ϕ , the bias is T ϕ T T 1 , which should be large enough to allow the benefits of the doubling of the erasure distance, which is between 10 to 100; see Fig. 5(c).

## a. Numbers

As reported in Ref. [15]; the cavity amplitude damping time is T 1 ≈ 250 µ s, and the intrinsic cavity dephasing is very long. The gate time is around 2 µ s, resulting in an erasure cavity error of 0 8%. . The erasure ancilla error is 1 3% as its . T 1 is 147 µ s. The measurement time is around 12 µ s, which results in erasure rate of 4 8% during the measurement. .

The subthreshold region is shown in Fig. 13 as a function of the measurement time and the two qubit gate time. As the results are comparable to the transmon dual-rail qubit (Fig. 7), the cavity dualrail qubit operations are not fast enough currently to show subthreshold behavior.

FIG. 13. The ( T , T 1 ϕ ) threshold curves. The dephasing time is chosen to be long enough to be negligible.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-a239d3a0-88b6-423d-a708-00fad9da931e.jpg)

One major advantage of the cavity dual-rail qubit over the transmon variant is the improved measurement fidelity, attributed to the longer cavity T 1 . Additionally, switching from a single cavity to a dualrail qubit enables erasure detection and single qubit operations, while in the case of the transmon, it also significantly extends dephasing time.

Unlike the transmon dual-rail qubit, there is no single-rail analog for the cavity, as the transition | 0 ⟩ ↔ | 1 ⟩ is extremely challenging for a single cavity.

## 2. Rydberg atoms

In this section, we analyze the threshold as a function of the number of ECs in the setting of Rydberg atoms. According to Ref. [6], an erasure conversion fraction of 0.98 could theoretically be attained, and we follow this bias in our analysis. This results in a maximum bias of 50, which could potentially be suf-

- [1] I. L. Chuang and Y. Yamamoto, Simple quantum computer, Physical Review A 52 , 3489 (1995).
- [2] E. Knill, R. Laflamme, and G. J. Milburn, A scheme for efficient quantum computation with linear optics, nature 409 , 46 (2001).
- [3] E. Knill, Bounds on the probability of success of postselected nonlinear sign shifts implemented with linear optics, Physical Review A 68 , 064303 (2003).
- [4] S. Bartolucci, P. Birchall, H. Bombin, H. Cable, C. Dawson, M. Gimeno-Segovia, E. Johnston, K. Kieling, N. Nickerson, M. Pant, et al. , Fusionbased quantum computation, Nature Communications 14 , 912 (2023).
- [5] M. Silva, M. R¨ otteler, and C. Zalka, Thresholds for linear optics quantum computing with photon loss at the detectors, Physical Review A-Atomic, Molecular, and Optical Physics 72 , 032307 (2005).
- [6] Y. Wu, S. Kolkowitz, S. Puri, and J. D. Thompson, Erasure conversion for fault-tolerant quantum computing in alkaline earth rydberg atom arrays, Nature Communications 13 , 4657 (2022).
- [7] M. Kang, W. C. Campbell, and K. R. Brown, Quantum error correction with metastable states of trapped ions using erasure conversion, PRX Quantum 4 , 020358 (2023).
- [8] A. M. Kubica and A. Retzker, Heralding of amplitude damping decay noise for quantum error correction, US Patent 11,748,652 (2023).
- [9] A. Kubica, A. Haim, Y. Vaknin, H. Levine, F. Brand˜o, and A. Retzker, Erasure qubits: Overa coming the T 1 limit in superconducting circuits, Physical Review X 13 , 041022 (2023).
- [10] J. D. Teoh, P. Winkel, H. K. Babla, B. J. Chap-

ficient to achieve a significant distance; see Fig. 5. Measurement-induced dephasing could in principle be as low as 10 -5 [11], though currently this was only achieved for the decay of the metastable state and not the Rydberg one. It is reasonable to assume that the coherence time of the qubit could be prolonged enough to not pose a limitation, and thus, the bias would only be limited by the selection rules, which justifies the use of the 50 × bias used in Fig. 14. The 2Q gate would be limited by power and control limitations, i.e., the faster the gate is, the higher is the Pauli error. In Fig. 14, we assume a negligible Pauli error, which is only dictated by the erasure conversion efficiency (50) and an erasure rate which is limited by the Rydberg state lifetime (assuming 100 µ s).

Fig. 14 indicates a consistent gain due to erasure which results from the low cost of measurement in this setting. However, measurement time does make a difference as can be seen by comparing Fig. 14(a) and Fig. 14(b).

- man, J. Claes, S. J. de Graaf, J. W. O. Garmon, W. D. Kalfus, Y. Lu, A. Maiti, K. Sahay, N. Thakur, T. Tsunoda, S. H. Xue, L. Frunzio, S. M. Girvin, S. Puri, and R. J. Schoelkopf, Dual-rail encoding with superconducting cavities, Proceedings of the National Academy of Sciences 120 , e2221736120 (2023).
- [11] S. Ma, G. Liu, P. Peng, B. Zhang, S. Jandura, J. Claes, A. P. Burgers, G. Pupillo, S. Puri, and J. D. Thompson, High-fidelity gates and mid-circuit erasure conversion in an atomic qubit, Nature 622 , 279-284 (2023).
- [12] P. Scholl, A. L. Shaw, R. B.-S. Tsai, R. Finkelstein, J. Choi, and M. Endres, Erasure conversion in a high-fidelity rydberg quantum simulator, Nature 622 , 273-278 (2023).
- [13] K. S. Chou, T. Shemma, H. McCarrick, T.-C. Chien, J. D. Teoh, P. Winkel, A. Anderson, J. Chen, J. Curtis, and S. J. e. a. de Graaf, Demonstrating a superconducting dual-rail cavity qubit with erasuredetected logical measurements, arXiv:2307.03169 (2023).
- [14] H. Levine and et al. , Demonstrating a longcoherence dual-rail erasure qubit using tunable transmons, Physical Review X 14 , 011051 (2024).
- [15] A. Koottandavida, I. Tsioutsios, A. Kargioti, C. R. Smith, V. R. Joshi, W. Dai, J. D. Teoh, J. C. Curtis, L. Frunzio, R. J. Schoelkopf, and M. H. Devoret, Erasure detection of a dual-rail qubit encoded in a double-post superconducting cavity, Phys. Rev. Lett. 132 , 180601 (2024).
- [16] E. Dennis, A. Kitaev, A. Landahl, and J. Preskill, Topological quantum memory, Journal of Mathe-

FIG. 14. The minimal required T 1 of the metastable manifold as a function of the 2Q gate time to reach subthreshold behavior. We assume that the bias between decay from the optical states to the ground state (and other detectable Rydberg states) is 50 times stronger than the decay back to the metastable subspace, introducing a 50 × bias between Pauli and erasure noise during the 2Q gate. We used 3ms as the decay time from the metastable state to the ground state. In the 'No-Erasure' simulation, we considered both decay to the metastable states and to the ground states as nondetectable Pauli noise. (a) The measurement is assumed to take 20 µ s, as in Ref. [11]. (b) The measurement is assumed to take 2 µ s. The measurement time does make a difference in this regime, especially when comparing 2 and 4 ECs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-a486baf7-f207-41aa-b69d-49176ce49983.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00829v2-figure-340f33d1-851f-4aa7-a81a-8f605c3985d1.jpg)

matical Physics 43 , 4452 (2002).

- [17] M. B. Hastings and J. Haah, Dynamically generated logical qubits, Quantum 5 , 564 (2021).
- [18] S. Gu, A. Retzker, and A. Kubica, Fault-tolerant quantum architectures based on erasure qubits (2023), arXiv:2312.14060 [quant-ph].
- [19] D. L. Campbell, Y.-P. Shim, B. Kannan, R. Winik, D. K. Kim, A. Melville, B. M. Niedzielski, J. L. Yoder, C. Tahan, S. Gustavsson, and W. D. Oliver, Universal nonadiabatic control of small-gap superconducting qubits, Phys. Rev. X 10 , 041051 (2020).
- [20] D. Gottesman, The heisenberg representation of
- quantum computers, in Proc. XXII International Colloquium on Group Theoretical Methods in Physics, 1998 (1998) pp. 32-43.
- [21] S. Aaronson and D. Gottesman, Improved simulation of stabilizer circuits, Physical Review A 70 , 052328 (2004).
- [22] R. Chao, M. E. Beverland, N. Delfosse, and J. Haah, Optimization of the surface code design for Majorana-based qubits, Quantum 4 , 352 (2020).
- [23] C. Gidney, Stim: a fast stabilizer circuit simulator, Quantum 5 , 497 (2021).
- [24] F. H. E. Watson and S. D. Barrett, Logical error rate scaling of the toric code, New Journal of Physics 16 , 093045 (2014).
- [25] M. Carroll, S. Rosenblatt, P. Jurcevic, I. Lauer, and A. Kandala, Dynamics of superconducting qubit relaxation times, npj Quantum Information 8 , 132 (2022).
- [26] Y. Zhou, Z. Zhang, Z. Yin, S. Huai, X. Gu, X. Xu, J. Allcock, F. Liu, G. Xi, Q. Yu, H. Zhang, M. Zhang, H. Li, X. Song, Z. Wang, D. Zheng, S. An, Y. Zheng, and S. Zhang, Rapid and unconditional parametric reset protocol for tunable superconducting qubits, Nature Communications 12 , 5924 (2021).
- [27] P. Magnard, P. Kurpiers, B. Royer, T. Walter, J.C. Besse, S. Gasparinetti, M. Pechal, J. Heinsoo, S. Storz, A. Blais, and A. Wallraff, Fast and unconditional all-microwave reset of a superconducting qubit, Physical Review Letters 121 , 060502 (2018).
- [28] K. Geerlings, Z. Leghtas, I. M. Pop, S. Shankar, L. Frunzio, R. J. Schoelkopf, M. Mirrahimi, and M. H. Devoret, Demonstrating a driven reset protocol for a superconducting qubit, Physical Review Letters 110 , 120501 (2013).
- [29] J. F. Marques, H. Ali, B. M. Varbanov, M. Finkel, H. M. Veen, S. L. M. van der Meer, S. VallesSanclemente, N. Muthusubramanian, M. Beekman, N. Haider, B. M. Terhal, and L. DiCarlo, Allmicrowave leakage reduction units for quantum error correction with superconducting transmon qubits, Physical Review Letters 130 , 250602 (2023).
- [30] E. Magesan and J. M. Gambetta, Effective hamiltonian models of the cross-resonance gate, Phys. Rev. A 101 , 052308 (2020).
- [31] R. Acharya, I. Aleiner, R. Allen, T. I. Andersen, M. Ansmann, et al. , Suppressing quantum errors by scaling a surface code logical qubit, Nature 614 , 676-681 (2023).
- [32] M. Steffen, J. Chow, S. Sheldon, and D. McClure, IBM Quantum's most performant system yet (2022).
- [33] C. K. Andersen, A. Remm, S. Lazar, S. Krinner, N. Lacroix, G. J. Norris, M. Gabureac, C. Eichler, and A. Wallraff, Repeated quantum error detection in a surface code, Nature Physics 16 , 875 (2020).
- [34] A. D. C´ orcoles, M. Takita, K. Inoue, S. Lekuch, Z. K. Minev, J. M. Chow, and J. M. Gambetta, Exploiting dynamic quantum circuits in a quantum algorithm with superconducting qubits, Physical Review Letters 127 , 100501 (2021).

- [35] L. Stefanazzi, K. Treptow, N. Wilcer, C. Stoughton, C. Bradford, S. Uemura, S. Zorzetti, S. Montella, G. Cancelo, S. Sussman, et al. , The QICK (Quantum Instrumentation Control Kit): Readout and control for qubits and detectors, Review of Scientific Instruments 93 , 044709 (2022).
- [36] O. Mukhanov, A. Kirichenko, C. Howington, J. Walter, M. Hutchings, I. Vernik, D. Yohannes, K. Dodge, A. Ballard, B. Plourde, et al. , Scalable quantum computing infrastructure based on superconducting electronics, in 2019 IEEE International Electron Devices Meeting (IEDM) (2019) pp. 31.2.1-31.2.4.
- [37] I. Zuk, D. Cohen, A. V. Gorshkov, and A. Retzker, Robust gates with spin-locked superconducting qubits, arXiv preprint arXiv:2306.09149 (2023).
- [38] K. Sahay, J. Jin, J. Claes, J. D. Thompson, and S. Puri, High-threshold codes for neutral-atom qubits with biased erasure errors, Physical Review X 13 , 041013 (2023).
- [39] J. P. Bonilla Ataides, D. K. Tuckett, S. D. Bartlett, S. T. Flammia, and B. J. Brown, The xzzx surface code, Nature Communications 12 , 1 (2021).
- [40] A. S. Darmawan, B. J. Brown, A. L. Grimsmo, D. K. Tuckett, and S. Puri, Practical quantum error cor-
- rection with the xzzx code and kerr-cat qubits, PRX Quantum 2 , 030345 (2021).
- [41] Q. Xu, N. Mannucci, A. Seif, A. Kubica, S. T. Flammia, and L. Jiang, Tailored xzzx codes for biased noise, Physical Review Research 5 , 013035 (2023).
- [42] A. Dua, A. Kubica, L. Jiang, S. T. Flammia, and M. J. Gullans, Clifford-deformed surface codes, PRX Quantum 5 , 010347 (2024).
- [43] N. P. Breuckmann and J. N. Eberhardt, Quantum low-density parity-check codes, PRX Quantum 2 , 040101 (2021).
- [44] K. Chang, S. Singh, J. Claes, K. Sahay, J. Teoh, and S. Puri, Surface code with imperfect erasure checks, arXiv:2408.XXXXX (2024).
- [45] O. Higgott and C. Gidney, Sparse blossom: correcting a million errors per core second with minimumweight matching, arXiv:2303.15933 (2023).
- [46] H. Levine, A. Keesling, G. Semeghini, A. Omran, T. T. Wang, S. Ebadi, H. Bernien, M. Greiner, V. Vuleti´, H. Pichler, and M. D. Lukin, Parallel imc plementation of high-fidelity multiqubit gates with neutral atoms, Phys. Rev. Lett. 123 , 170503 (2019).
- [47] M. R. Geller and Z. Zhou, Efficient error models for fault-tolerant architectures and the pauli twirling approximation, Phys. Rev. A 88 , 012314 (2013).
Shouzhen Gu Yotam Vaknin Alex Retzker Aleksander Kubica
Optimizing quantum error correction protocols with erasure qubits
Erasure qubits offer a promising avenue toward reducing the overhead of
quantum error correction (QEC) protocols. However, they require additional
operations, such as erasure checks, that may add extra noise and increase
runtime of QEC protocols. To assess the benefits provided by erasure qubits, we
focus on the performance of the surface code as a quantum memory. In
particular, we analyze various erasure check schedules, find the correctable
regions in the phase space of error parameters and probe the subthreshold
scaling of the logical error rate. We then consider a realization of erasure
qubits in the superconducting hardware architectures via dual-rail qubits. We
use the standard transmon-based implementation of the surface code as the
performance benchmark. Our results indicate that QEC protocols with erasure
qubits can outperform the ones with state-of-the-art transmons, even in the
absence of precise information about the locations of erasure errors.