Abstract
Protecting quantum information from errors is essential for large-scale quantum computation. Quantum error correction (QEC) encodes information in entangled states of many qubits and performs parity measurements to identify errors without destroying the encoded information. However, traditional QEC cannot handle leakage from the qubit computational space. Leakage affects leading experimental platforms, based on trapped ions and superconducting circuits, which use effective qubits within many-level physical systems. We investigate how two-transmon entangled states evolve under repeated parity measurements and demonstrate the use of hidden Markov models to detect leakage using only the record of parity measurement outcomes required for QEC. We show the stabilization of Bell states over up to 26 parity measurements by mitigating leakage using postselection and correcting qubit errors using Pauli-frame transformations. Our leakage identification method is computationally efficient and thus compatible with real-time leakage tracking and correction in larger quantum processors.
INTRODUCTION
Large-scale quantum information processing hinges on overcoming errors from environmental noise and imperfect quantum operations. Fortunately, the theory of quantum error correction (QEC) predicts that the coherence of single degrees of freedom (logical qubits) can be better preserved by encoding them in ever-larger quantum systems (Hilbert spaces), provided the error rate of the constituent elements lies below a fault-tolerance threshold (1). Experimental platforms based on trapped ions and superconducting circuits have achieved error rates in single-qubit (SQ) gates (2–4), two-qubit gates (2, 4, 5), and qubit measurements (3, 6–8) at or below the threshold for popular QEC schemes such as surface (9, 10) and color codes (11). They therefore seem well poised for the experimental pursuit of quantum fault tolerance. However, a central assumption of textbook QEC that error processes can be discretized into bit flips (X), phase flips (Z), or their combination (Y = iXZ) only is difficult to satisfy experimentally. This is due to the prevalent use of many-level systems as effective qubits, such as hyperfine levels in ions and weakly anharmonic transmons in superconducting circuits, making leakage from the two-dimensional computational space of effective qubits a threatening error source. In quantum dots and trapped ions, leakage events can be as frequent as qubit errors (12, 13). However, even when leakage is less frequent than qubit errors as in superconducting circuits (2, 5), if ignored, then leakage can produce the dominant damage to encoded logical information. To address this, theoretical studies propose techniques to reduce the effect of leakage by periodically moving logical information and removing leakage when qubits are free of logical information (14–17). Alternatively, more hardware-specific solutions have been proposed for trapped ions (18) and quantum dots (19). In superconducting circuits, recent experiments have demonstrated single- and multiround parity measurements to correct qubit errors with up to nine physical qubits (20–28). Parallel approaches encoding information in the Hilbert space of single resonators using cat (29) and binomial codes (30) used transmon-based photon parity checks to approach the breakeven point for a quantum memory. However, no experiment has demonstrated the ability to detect and mitigate leakage in a QEC context.
Here, we experimentally investigate leakage detection and mitigation in a minimal QEC system. Specifically, we protect an entangled state of two transmon data qubits (QDH and QDL) from qubit errors and leakage during up to 26 rounds of parity measurements via an ancilla transmon (QA). Performing these parity checks in the Z basis protects the state from X errors, while interleaving checks in the Z and X bases protects it from general qubit errors (X, Y, and Z). Leakage manifests itself as a round-dependent degradation of data-qubit correlations ideally stabilized by the parity checks: 〈Z ⊗ Z〉 in the first case and 〈X ⊗ X〉, 〈Y ⊗ Y〉, and 〈Z ⊗ Z〉 in the second. We introduce hidden Markov models (HMMs) to efficiently detect data-qubit and ancilla leakage, using only the string of parity outcomes, demonstrating restoration of the relevant correlations. Although we use postselection here, the low technical overhead of HMMs makes them ideal for real-time leakage correction in larger QEC codes.
RESULTS
A minimal QEC setup
Repetitive parity checks can produce and stabilize two-qubit entanglement. For example, performing a Z ⊗ Z parity measurement (henceforth a ZZ check) on two data qubits prepared in the unentangled state ∣++〉 = (∣0〉 + ∣1〉) ⊗ (∣0〉 + ∣1〉) /2 will ideally project them to either of the two (entangled) Bell states
Generating entanglement by measurement
Our parity check is an indirect quantum measurement involving coherent interactions of the data qubits with QA and subsequent QA measurement (Fig. 1A) (33). The coherent step maps the data-qubit parity onto QA in 120 ns using SQ and two-qubit controlled-phase (CZ) gates (5). Gate characterizations (see the Supplementary Materials) indicate state-of-the-art gate errors eSQ = {0.08 ± 0.02,0.14 ± 0.016,0.21 ± 0.06}% and eCZ = {1.4 ± 0.6,0.9 ± 0.16}% with leakage per CZ L1 = {0.27 ± 0.12,0.15 ± 0.07}%. We measure QA with a 620-ns pulse including photon depletion (7, 34), achieving an assignment error ea = 1.0 ± 0.1%. We avoid data-qubit dephasing during the QA measurement by coupling each qubit to a dedicated readout resonator (RR) and a dedicated Purcell resonator (PR) (fig. S1) (8). The parity check has a cycle time of 740 ns, corresponding to only 2.5 ± 0.2% and 5.0 ± 0.3% of the data-qubit echo dephasing times (see the Supplementary Materials).
(A) Quantum circuit for a parity measurement of the data qubits via coherent operations with ancilla QA and QA measurement. Tomography reconstructs the data-qubit output density matrix (ρ). Echo pulses (orange) are applied halfway the QA measurement when performing tomography sequential to the QA measurement. (B) Bloch sphere representation of the even-parity subspace with a marker on ∣Φ+⟩. (C to F) Plots of ρ with fidelity to the Bell states (indicated by frames) for tomography simultaneous with QA measurement (C to E) and sequetial to QA measurement (F). (C)[(D)] Conditioning on MA = +1[ −1] ideally generates ∣Φ+⟩ [∣Ψ+⟩] with equal probability P. (E)[(F)] PFU applies bit-flip correction (X on QDH) for MA = −1 and reconstructs ρ using all data for simultaneous [sequential] tomography.
The parity measurement performance can be quantified by correlating its outcome with input and output states. We first quantify the ability to distinguish even-parity (∣00〉, ∣11〉) from odd-parity (∣01〉, ∣10〉) data-qubit input states, finding an average parity assignment error ea,ZZ = 5.1 ± 0.2%. Second, we assess the ability to project onto the Bell states by performing a ZZ check on ∣++〉 and reconstructing the most likely physical data-qubit output density matrix ρ, conditioning on MA = ±1. When tomographic measurements are performed simultaneously with the QA measurement, we find Bell-state fidelities F∣Φ+〉∣MA=+1 = 〈Φ+∣ρMA=+1∣Φ+〉 = 94.7 ± 1.9% and F∣Ψ+〉∣MA=−1 = 94.5 ± 2.5% (Fig. 1, C and D). We connect ∣Ψ+〉 to ∣Φ+〉 by incorporating the PFU into the tomographic analysis, obtaining F∣Φ+〉 = 94.6 ± 0.9% without any postselection (Fig. 1E). The nondemolition character of the ZZ check is then validated by performing tomography only once the QA measurement completes. We include an echo pulse on both data qubits during the QA measurement to reduce intrinsic decoherence and negate residual coupling between data qubits and QA (fig. S3). The degradation to F∣Φ+〉 = 91.8 ± 0.5% is consistent with intrinsic data-qubit decoherence under echo and confirms that measurement-induced errors are minimal.
Protecting entanglement from bit flips and the observation of leakage
QEC stipulates repeated parity measurements on entangled states. We therefore study the evolution of F∣Φ+〉 = (1 + 〈X ⊗ X〉 − 〈Y ⊗ Y〉 + 〈Z ⊗ Z〉)/4 and its constituent correlations as a function of the number M of checks (Fig. 2A). When performing PFU using the first ZZ outcome only (ignoring subsequent outcomes), we observe that F∣Φ+〉 witnesses entanglement (>0.5) during 10 rounds and approaches randomization (0.25) by M = 25 (Fig. 2B). The constituent correlations also decay with simple exponential forms. A best fit of the form 〈Z ⊗ Z〉[M] = a · e−M/υZZ + b gives a decay time υZZ = 9.0 ± 0.9 rounds; similarly, we extract υXX = 11.7 ± 1.0 rounds (Fig. 2, C and D). By comparison, we observe that Bell states evolving under dynamical decoupling only (no ZZ checks; see fig. S4) decay similarly (υZZ = 8.6 ± 0.3, υXX = 12.8 ± 0.4 rounds). These similarities indicate that intrinsic data-qubit decoherence is also the dominant error source in this multiround protocol.
(A) The quantum circuit of Fig. 1A extended with M rounds of repeated ZZ checks. (B) Fidelity to ∣Φ+⟩ as a function of M. “No error” postselects the runs in which no bit flip is detected. “Final” applies PFU based on the last two outcomes (equivalent to MWPM). “First” uses the first parity outcome only. “Idling DD” are Bell states evolving under dynamical decoupling only (quantum circuit in fig. S4). (C) Corresponding 〈X ⊗ X〉. “Final” coincides with “First”. (D) Corresponding 〈Z ⊗ Z〉. The weak degradation observed for “Final” is the hallmark of leakage. Curves in (B) to (D) are best fits of a simple exponential decay.
To demonstrate the ability to detect X and Y but not Z errors, we condition the tomography on signaling no errors during M rounds. This boosts 〈Z ⊗ Z〉 to a constant, while the undetectability of Z errors only allows slowing the decay of 〈X ⊗ X〉 to υXX = 33.2 ± 1.7 rounds (and of 〈Y ⊗ Y〉 to υYY = 31.3 ± 1.9 rounds). Naturally, this conditioning comes at the cost of the postselected fraction fpost reducing with M (fig. S5).
Moving from error detection to correction, we consider the protection of ∣Φ+〉 by tracking X errors and applying corrections in postprocessing. The correction relies on the final two MA only, concluding even parity for equal measurement outcomes and odd parity for unequal. For this small-scale experiment, this strategy is equivalent to a decoder based on minimum-weight perfect matching (MWPM) (10, 35), justifying its use. Because our PFU strategy converts Y errors into Z errors, one expects a faster decay of 〈X ⊗ X〉 compared to the no-error conditioning; we observe υXX = 11.8 ± 1.0 rounds. Correction should lead to a constant 〈Z ⊗ Z〉. While 〈Z ⊗ Z〉 is clearly boosted, a weak decay to a steady state 〈Z ⊗ Z〉 = 0.73 ± 0.03 is also evident (Fig. 2D). As previously observed in (31), this degradation is the hallmark of leakage [see also (21, 24)]. We additionally compare the experimental results to simulations using a model that assumes ideal two-level systems (35) (no leakage) based on independently calibrated parameters of table S1 (fig. S8, A to D). At M = 1, the model and the experiment coincide for all correction strategies. At larger M, “first” and “final” correction strategies deviate substantially, consistent with a gradual buildup of leakage, which we now turn our focus to.
Leakage detection using HMMs
Both ancilla and data-qubit leakage in our experiment can be inferred from a string
(A) Simplified HMM. In each round, a hidden state (leaked or unleaked) (top) is updated probabilistically (full arrows) and produces an observable MA (bottom) with state-dependent probabilities (dashed arrows). After training, the HMM can be used to assess the likelihood of states given a produced string
Neither of the above patterns is entirely unique to leakage; each may also be produced by some combination of qubit errors. Therefore, we cannot unambiguously diagnose an individual experimental run of corruption by leakage. However, given a set of ancilla measurements MA[0], …, MA[m], the likelihood
(A) Simplified quantum circuit with preparation, repeated pairs of ZZ and XX checks, and data-qubit tomography. (B) Fidelity to ∣Φ+⟩ as a function of M, extracted from the data-qubit tomography. “No error” postselects the runs in which no error is detected (postselected fraction in fig. S5). “Final” applies PFU based on the last three outcomes (equivalent to MWPM). “Final + HMM” includes mitigation of leakage. “First” uses only the first pair of parity outcomes. (C and D) Corresponding 〈X ⊗ X〉 and 〈Z ⊗ Z〉. Curves in (B) to (D) are best fits of a simple exponential decay.
Before assessing the ability of our HMMs to improve fidelity in a leakage mitigation scheme, we first validate and benchmark them internally. A common method to validate the HMM’s ability to model the experiment is to compare statistics of the experimentally generated data to a simulated dataset generated by the model itself. As we are concerned only with the ability of the HMM to discriminate leakage,
All values are obtained by optimizing the likelihood of observing the given syndrome data except for the ancilla leakage rate (denoted as *), which is directly obtained from the experiments (as noted in the main text).
The above validation suggests that we may assume that the ratio of actual leakage events at a given Lcomp,Q is well approximated by Lcomp,Q itself (which is true for the simulated data). Under this assumption, we expose the HMM discrimination ability by plotting its receiver operating characteristic (ROC) (38). The ROC (Fig. 3F) is a parametric plot (sweeping a threshold
We now verify and externally benchmark our HMMs by their ability to improve 〈Z ⊗ Z〉 by rejecting data with a high probability of leakage. To do this, we set a threshold
Protecting entanglement from general qubit errors and mitigation of leakage
We lastly demonstrate leakage mitigation in the more interesting scenario where ∣Φ+〉 is protected from general qubit error by interleaving ZZ and XX checks (28, 31). ZZ may be converted to XX by adding π/2 y rotations on the data qubits simultaneous with those on QA. This requires that we change the definition of the syndrome to sD[m] = MA[m] · MA[m − 1] · MA[m − 2] · MA[m − 3], as we need to “undo” the interleaving of the ZZ and XX checks to detect errors. For an input state
DISCUSSION
This HMM demonstration provides exciting prospects for leakage detection and correction. In larger systems, independent HMMs can be dedicated to each qubit because leakage produces local error signals (16). An HMM for an ancilla only needs its measurement outcomes, while a data-qubit HMM only needs the outcomes of the nearest-neighbor ancillas (details in see the Supplementary Materials). Therefore, the computational power grows linearly with the number of qubits, making the HMMs a small overhead when running parallel to MWPM. HMM outputs could be used as inputs to MWPM, allowing MWPM to dynamically adjust its weights. The outputs could also be used to trigger leakage reduction units (14–17) or qubit resets (39).
In summary, we have performed the first experimental investigation of leakage detection during repetitive parity checking, successfully protecting an entangled state from qubit errors and leakage in a circuit quantum electrodynamics processor. Future work will extend this protection to logical qubits, e.g., the 17-qubit surface code (35, 40). The low technical overhead and scalability of HMMs are attractive for performing leakage detection and correction in real time using the same parity outcomes as traditionally used to correct qubit errors only.
MATERIALS AND METHODS
Device
Our quantum processor (fig. S1) follows a three-qubit–frequency extensible layout with nearest-neighbor interactions that is designed for the surface code (41). Our chip contains low- and high-frequency data qubits (QDL and QDH) and an intermediate-frequency ancilla (QA). SQ gates around axes in the equatorial plane of the Bloch sphere are performed via a dedicated microwave drive line for each qubit. Two-qubit interactions between nearest neighbors are mediated by a dedicated bus resonator (extensible to four per qubit) and controlled by individual tuning of qubit transition frequencies via dedicated flux-bias lines (42). For measurement, each qubit is dispersively coupled to a dedicated RR, which is itself connected to a common feedline via a dedicated PR. The RR-PR pairs allow frequency-multiplexed readout of selected qubits with negligible backaction on untargeted qubits (8).
Hidden Markov models
HMMs provide an efficient tool for indirect inference of the state of a system given a set of output data (36). An HMM describes a time-dependent system as evolving between a set of Nh hidden states {h} and returning one of No outputs {o} at each timestep m. The evolution is stochastic: the system state H[m] of the system at timestep m depends probabilistically on the state H[m − 1] at the previous timestep, with probabilities determined by a Nh × Nh transition matrix A
The user cannot directly observe the system state and must infer it from the outputs O[m] ∈ {o} at each timestep m. This output is also stochastic: O[m] depends on H[m] as determined by a No × Nh output matrix B
If the A and B matrices are known, along with the expected distribution
HMMs for QEC experiments
To maximize the discrimination ability of HMMs in the various settings studied in this work, we choose different quantities to use for our output vectors
One may predict the computational likelihood for data-qubit (D) leakage at timestep M in the ZZ-check experiment given
However, in the repeated ZZ-check experiment, the ancilla (A) needs to be within the computational subspace for two rounds to perform a correct parity measurement. Therefore, the computational likelihood is slightly more complicated to calculate
In the interleaved ZZ- and XX-check experiment, the situation is more complicated as we require data from the final two parity checks to fully characterize the quantum state. This implies that we need unleaked data qubits for the last two rounds and unleaked ancillas for the last three. The likelihood of the latter may be calculated by similar means to the above.
Simplest models for leakage discrimination
One need not capture the full dynamics of the quantum system in an HMM to infer whether a qubit is leaked. This is of critical importance if we wish to extend this method for the purposes of leakage mitigation in a large QEC code (as we discuss in the Supplementary Materials). The simplest possible HMM (Fig. 3A) has two hidden states: H[m] = 1 if the qubit(s) in question are within the computational subspace and H[m] = 2 if QA (or either data qubit) is leaked. The labels 1 and 2 are arbitrary here and explicitly have no correlation with the qubit states ∣1⟩ and ∣2⟩. Then, the 2 × 2 transition matrix simply captures the leakage and seepage rates of the system in question
The 2 × 2 output matrices then capture the different probabilities of seeing output O[m] = 0 or O[m] = 1 when the qubit(s) are leaked or unleaked
When studying data-qubit leakage, p0,1 simply captures the rate of errors within the computational subspace. Then, in the repeated ZZ-check experiment, p1,0 captures events such as ancilla or measurement errors that cancel the error signal of a leakage event. However, in the interleaved ZZ and XX experiments, a leaked qubit causes the syndrome to be random, so we expect p1,0 ∼ 0.5. When studying ancilla leakage, p1,0 is simply the probability of ∣2⟩ state being read out as ∣0⟩ and is also expected to be close to 0. However, p0,1 ∼ 0.5, as we do not reset QA or the logical state between rounds of measurement, and thus, any measurement in isolation is roughly equally likely to be 0 or 1. In all situations, we assume that the system begins in the computational subspace—πn(0) = δn,0. With this fixed, we may choose the parameters pleak, pseep, p0,1, and p1,0 to maximize the likelihood
Modeling additional noise
The simple model described above does not completely capture all of the details of the stabilizer measurements
Here, we choose the matrices
The error generators D(A) and D(B) may be identified as derivatives of A with respect to these error rates
This may be extended to calculate derivatives of the likelihood
HMMs used in Figs. 2 and 4
Different Markov models (with independently optimized parameters) were used to optimize ancilla and data-qubit leakage estimation for both the ZZ experiment and the experiment interleaving ZZ and XX checks. This leads to a total of four HMMs, which we label HZZ-D, HZZ-A, HZZ, XX-D, and HZZ, XX-A. A complete list of parameter values used in each HMM is given in Table 1. We now describe the features captured by each HMM. As we show in the Supplementary Materials, these additional features are not needed to increase the error mitigation performance of the HMMs but rather to ensure their closeness to the experiment and increase trust in their internal metrics.
To go beyond the simple HMM in the ZZ-check experiment when modeling data-qubit leakage (HZZ-D), we need to include additional states to account for the correlated signals of ancilla and readout errors. If we assume data-qubit errors (that remain within the logical subspace) are uncorrelated in time, then they are already well captured in the simple model. This is because any single error on a data qubit may be decomposed into a combination of Z errors (which commute with the measurement and thus are not detected) and X errors (which anticommute with the measurement and thus produce a single error signal sD[m] = 1) and is thus captured by the p0,1 parameter. When one of the data qubits is leaked, uncorrelated X errors on the other data qubit cancel the constant sD[m] = −1 signal for a single round and are thus captured by the p1,0 parameter. However, errors on the ancilla, and readout errors, give error signals that are correlated in time (separated by one or two timesteps, respectively). This may be accounted for by including extra “ancilla error states.” These may be most easily labeled by making the h labels a tuple h = (h0, h1), where h0 keeps track of whether or not the qubit is leaked, and h1 = 1,2,3 keeps track of whether or not a correlated error has occurred. In particular, we encode the future syndrome for 2 cycles in the absence of error on h1, allowing us to account for any correlations up to two rounds in the future. This extends the model to a total of 4 × 2 = 8 states. The transition and output matrices in the absence of error for the unleaked h0 = 0 states may then be written in a compact form (noting that leakage errors cancel out with correlated ancilla and readout errors to give sD[m] = +1)
Let us briefly demonstrate how the above works for an ancilla error in the system. Suppose the system was in the state h = (0,3) at time m, it would output MA[m] = −1 and then evolve to h = (0,3//2) = (0,1) at time m + 1 (in the absence of additional error). Then, it would output a second error signal [MA[m + 1] = −1] and lastly decay back to the h = (0,1//2) = (0,0) state. This gives the HMM the ability to model the ancilla error as an evolution from h = (0,0) to h(0,3). Formally, we assign the matrix
The corresponding error rate pancilla is then an additional free parameter to be optimized to maximize the likelihood. To finish the characterization of this error channel, we need to consider the effect of the ancilla error in states other than h = (0,0). Two ancilla errors in the same timestep cancel, but two ancilla errors in subsequent timesteps will cause the signature sD = …, −1, +1, −1, …. This may be captured by an evolution from h = (0,2) to h = (0,3) [instead of h = (0,1)], which implies that we should set
[Note that A0 already captures a decay from h = (0,2) → (0,1) → (0,0), which will give the desired signal.] We note that this also matches the signature of readout error, which can then be captured by a separate error channel
One can then check that ancilla errors in h = (0,3) should cause the system to remain in h = (0,3) and that ancilla or readout errors in h = (0,1) should evolve the system to h = (0,2). We note that this model cannot account for the sD = … −1, +1, +1, +1, −1 signature of readout error at time m and m + 2, but adjusting the model to include this has negligible effect.
An ancilla error in the ZZ-check experiment when the data qubits are leaked has the same correlated behavior as when they are not but may occur at a different rate. This requires that we define a new matrix
We also assume that leakage pleak and seepage pseep rates are independent of these correlated errors (i.e.,
There are a few key differences between the interleaved ZZ-XX and ZZ experiments that need to be captured in the data-qubit HMM HZZ,XX-D. First, as the syndrome is now given by sD[m] = MA[m] · MA[m − 1] · MA[m − 2] · MA[m − 3], ancilla and classical readout errors can then generate a signal stretching up to four steps in time. This implies that we require 24 possibilities for h1 to keep track of all correlations. However, as a leaked data qubit makes ancilla output random in principle, we no longer need to keep track of the ancilla output upon leakage. This implies that we can accurately model the experiment with 16 + 1 = 17 states, which we can label by h ∈ {2, (1, h1)}. The A0 and B0 matrices in the unleaked states (1, h1) follow Eq. 15, anwd we fix [A0]2,2 = 1 (as in the absence of pseep, a leaked state stays leaked). However, we allow for some bias in the leaked state error rate—B−1,2 = pdata,leaked is not fixed to 0.5. For example, this accounts for a measurement bias toward a single state, which will reduce the error rate below 0.5. The nonzero elements in the matrices
Here, a ⊕ b refers to the addition of each binary digit of a and b modulo 2. We may use this formalism to additionally keep track of Y data-qubit errors, which show up as correlated errors on subsequent XX and ZZ stabilizer checks, by introducing a new error channel
As the output used for the HZZ-A HMM is the pure measurement outcome MA, the dominant signal that must be accounted for is that of the stabilizer ZZ itself. This causes either a constant signal MA[m] = MA[m − 1] or a constant flipping signal MA[m] = −MA[m − 1]. This cannot be accounted for in the simple HMM, as it cannot contain any history in a single unleaked state. To deal with this, we extend the set of states in the HZZ-A HMM to include both an estimate of the ancilla state a ∈ {0,1,2} at the point of measurement and the stabilizer state s ∈ {0,1} and label the states by the tuple (a, s). The ancilla state then immediately defines the device output in the absence of any error
The only thing that affects the output matrices is readout error
Data-qubit errors flip the stabilizer with probability pdata
Ancilla errors flip the ancilla with probability pancilla, but these are dominated by T1 decay and so are highly asymmetric. To account for this, we used different error rates panc,a,a′ for the four possible combinations of ancilla measurement at time m − 1 and expected ancilla measurement at time m
(Note that this asymmetry could not be accounted for in the data-qubit HMM as the state of the ancilla was not contained within the output vector.) As with the data-qubit HMMs, we assume that ancilla leakage is HMM state independent, as it is dominated by CZ gates during the time that the ancilla is either in ∣+⟩ or ∣−⟩. We also assume that leakage (with rate pleak) and seepage (with rate pseep) have equal chances to flip the stabilizer state, as ancilla leakage has a good chance to cause additional error on the data qubits.
The ancilla-qubit HMMs need little adjustment between the ZZ-check experiment and the experiment interleaving ZZ and XX checks. The HZZ,XX-A HMM behaves almost identically to the HZZ-A HMM, but we include in the state information on the XX stabilizer and the ZZ stabilizer. This leaves the states indexed as (a, s1, s2). The HMM needs to also keep track of which stabilizer is being measured. This may be achieved by shuffling the stabilizer labels at each timestep: For a = 0,1, we set
Other than this, the HMM follows the same equations as above (with the additional index added as expected).
Uncertainty calculations
All quoted uncertainties are an estimation of SEM. SEMs for the independent device characterizations (see the “Generating entanglement by measurement” section and table S1) are either obtained from at least three individually fitted repeated experiments (
SEMs in the first-round Bell-state fidelities (Fig. 1 and fig. S6 and see the “Generating entanglement by measurement,” “Protecting entanglement from bit flips and the observation of leakage,” “Leakage detection using HMMs,” and “Protecting entanglement from general qubit errors and mitigation of leakage” sections) are obtained through bootstrapping. For bootstrapping, a dataset (in total 4096 runs with each 36 tomographic elements and 28 calibration points) is subdivided into four subsets, and tomography is performed on each of these subsets individually. As verification, subdivision was performed with eight subsets leading to similar SEMs. SEMs in the multiround experiment parameters (steady-state fidelities and decay constants) are also estimated from least-squares fitting by the LmFit fitting module using the covariance matrix (see the “Protecting entanglement from bit flips and the observation of leakage,” “Leakage detection using HMMs,” and “Protecting entanglement from general qubit errors and mitigation of leakage” sections) (44).
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/12/eaay3050/DC1
Supplementary Materials and Methods
Supplementary Text
Fig. S1. Quantum processor.
Fig. S2. Complete wiring diagram of electronic components inside and outside the 3He/4He dilution refrigerator (Leiden Cryogenics CF-CS81).
Fig. S3. Study of data-qubit coherence and phase accrual during ancilla measurement.
Fig. S4. Quantum circuit for Bell-state idling experiments under dynamical decoupling.
Fig. S5. Postselected fractions for the “no error” conditioning in Figs. 2 and 4.
Fig. S6. Generating entanglement by sequential ZZ and XX parity measurements and PFU.
Fig. S7. ROCs for mitigation of data-qubit and ancilla leakage during interleaved ZZ and XX checks.
Fig. S8. Comparison of experimental data and no-leakage modeling of the repeated parity check experiments of Figs. 2 and 4.
Fig. S9. Leakage mitigation for the repeated parity check experiments as a function of the chosen threshold.
Fig. S10. Leakage mitigation for the simple, two-state HMMs for repeated parity check experiments as a function of the chosen threshold.
Fig. S11. ROCs for leakage mitigation as in fig. S7 but using simple two-state HMMs.
Table S1. Measured parameters of the three-transmon device.
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).