Demonstration of measurement-free universal logical quantum computation
Friederike Butt, Ivan Pogorelov, Robert Freund, Alex Steiner, Marcel Meyer, Thomas Monz, Markus Müller
TL;DR
This paper demonstrates a new approach to quantum computing that avoids mid-circuit measurements, enabling fault-tolerant operations and algorithms.
Contribution
A measurement-free universal logical quantum computation toolbox is proposed and demonstrated experimentally on a trapped-ion processor.
Findings
Modular logical state teleportation is achieved without mid-circuit measurements.
A fault-tolerant universal gate set is realized using coherent operations only.
Grover’s algorithm is implemented fault-tolerantly on three logical qubits encoded in eight physical qubits.
Abstract
The ability to perform quantum error correction (QEC) and robust gate operations on encoded qubits opens the door to demonstrations of quantum algorithms. Contemporary QEC schemes typically require mid-circuit measurements with feed-forward control, which are challenging for qubit control, often slow, and susceptible to relatively high error rates. In this work, we propose and experimentally demonstrate a universal toolbox of fault-tolerant logical operations on error-detecting codes without mid-circuit measurements on a trapped-ion quantum processor. We present modular logical state teleportation between two four-qubit error-detecting codes without measurements during algorithm execution. Moreover, we realize a fault-tolerant universal gate set on an eight-qubit error-detecting code hosting three logical qubits, based on state injection, which can be executed by coherent gate…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 10
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9- —100011039ODNI | Intelligence Advanced Research Projects Activity (IARPA)
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsQuantum Computing Algorithms and Architecture · Quantum Information and Cryptography · Quantum-Dot Cellular Automata
Introduction
The practical implementation of quantum algorithms depends on their resilience to errors, alongside the ability to perform arbitrary quantum operations. Quantum error correction (QEC) enables the detection and correction of errors arising during computation by encoding information across multiple physical qubits^1–3^. Computations on these encoded qubits can be realized through a discrete, universal set of gates^3^. These operations have to be implemented in a robust, fault-tolerant (FT) fashion, meaning that local faults in the underlying gate operations do not proliferate uncontrollably across the logical qubits^4^. However, no QEC code intrinsically supports a full, inherently FT universal gate set^5^. Completing this FT universal gate set is a key challenge for realizing a potential advantage beyond the reach of algorithms that can be efficiently simulated classically. Recent experiments have demonstrated QEC cycles on trapped-ion quantum processors^6–10^, superconducting architectures^11–14^, as well as neutral-atom platforms^15,16^. FT universal gate sets have been realized on these platforms by means of code switching^17,18^, where information is transferred between two codes with complementary sets of inherently FT gates, as well as magic state injection^12,19–21^, which requires high-fidelity magic states as a resource^22,23^. These advancements in the practical and scalable implementations of logical qubits enabled the execution of the first, small quantum algorithms run on encoded qubits, such as the Bernstein–Vazirani algorithm^24,25^, one-bit addition^26,27^, Grover search on two logical qubits^16,28^ or the quantum Fourier transform on three logical qubits^29^.
Many practical protocols rely on measurements during algorithm execution and feed-forward operations conditioned on these measurement outcomes, which is experimentally demanding on many hardware platforms and limits their success probability: In both atomic and superconducting quantum processors, measurements remain orders of magnitude slower than typical gate times, which poses speed limitations and results in decoherence of idling qubits during measurements. Moreover, fluorescence read-out in atomic setups requires additional cooling during and after measurements, as atoms are heated during this process^17,20,30–32^.
Following early works^33,34^, recent theoretical works have proposed practical measurement-free protocols for logical state preparation^35^, rounds of QEC^36–38^ and the implementation of an FT universal gate set^39–41^. In these protocols, stabilizer information is transferred onto auxiliary qubits, allowing decoding and coherent feedback to be carried out within the quantum algorithm itself. This approach avoids the need for mid-circuit measurements or feed-forward operations entirely. At the end, auxiliary qubits are replaced or reset to be reused, effectively removing the entropy introduced by the noise.
In this work, we develop and experimentally demonstrate a complete toolbox of logical operations needed for FT universal quantum computing on an ion-trap quantum processor, without mid-circuit measurements or feed-forward operations. First, we construct protocols for modular logical quantum state teleportation, such that different encoded blocks are never directly coupled to one another, which is a key desideratum for scaling up quantum computations to large numbers of logical qubits. We analyze the performance of these measurement-free protocols for different logical input states, accompanied by numerical simulations. We then complete an FT, measurement-free universal gate set for an eight-qubit error-detecting code by constructing and implementing circuits for a logical Hadamard gate on an encoded qubit. Finally, we use this implementation as a building block for Grover’s algorithm to search for two elements out of eight, for the first time demonstrating a small-scale FT and measurement-free universal quantum algorithm.
Results
Experimental setup
The experimental data were obtained with a 16-qubit quantum computing device based on trapped ions^42^. The chain of 16 ^40^Ca^+^ ions is confined in a linear Paul trap. The physical qubits are encoded in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 0\rangle=| {4}^{2}{{{{\rm{S}}}}}_{1/2},{m}_{J}=-1/2 \rangle$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 1\rangle=| {3}^{2}{{{{\rm{D}}}}}_{5/2},{m}_{J}=-1/2\rangle$$\end{document} Zeeman sub-levels. The state of each qubit can be manipulated individually by optically addressing the ions with 729 nm laser light. Two-qubit gates are realized as a Mølmer–Sørensen (MS) interaction^43^, providing all-to-all two-qubit-gate connectivity. Overall, the native gate set of the device includes arbitrary-angle rotation gates \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R(\theta,\phi )=\exp (-i\frac{\theta }{2}[X\cos \phi+Y\sin \phi ])$$\end{document} , ‘virtual’ Z-gates \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}_{Z}(\theta )=\exp (-i\frac{\theta }{2}Z)$$\end{document} , and maximally-entangling two-qubit gates \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$XX(\pi /2)=\exp (-i\frac{\pi }{4}X\otimes X)$$\end{document} . A description of the experimental setup can be found in refs. ^21,42,44^.
Our trapped-ion platform is capable of performing mid-circuit measurement operations, as was shown in ref. ^7^. However, such an operation, together with an additional feed-forward, represents a substantial experimental overhead in both sequence duration and infidelity. In our protocols, we do not need to perform these operations, but require only resets of the quantum state of certain qubits. The reset is faster than our current implementation of the mid-circuit measurement (1.7 ms vs. ≈30 ms), and the preservation of the data qubit’s state is higher, which is discussed further in the ‘Methods’ Subsection ‘Qubit reset’. Instead of re-initializing physical qubits, one can also replace them with fresh physical qubits. In our experiments, we make use of the full 16-ion register and use fresh physical qubits whenever possible.
Logical state teleportation without mid-circuit measurements
In this section, we discuss how to teleport a logical state between two four-qubit registers without mid-circuit measurements or feed-forward operations, and demonstrate this concept experimentally. We consider a [[4, 1, 2]]-code instance that encodes k = 1 logical qubit in n = 4 physical qubits and has distance d = 2, meaning that any single error can be detected. The stabilizers and logical operators defining the code are shown in Fig. 1a.Fig. 1. Measurement-free logical state teleportation with the [[4, 1, 2]]-code.a Stabilizers SZ, SX and logical operators of the [[4, 1, 2]]-code. b Experimental logical quantum state tomography for FT logical state initialization. The black dashed boxes correspond to ideal values in a fault-free case. Any deviation from the black dashed boxes, e.g., unboxed gray, blue or brown areas, indicates that noise reduces the overall fidelity. c High-level circuit for measurement-based modular logical teleportation. The source (S) and target (T) code blocks are merged by measuring the joint logical \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}_{{{{\rm{L}}}}}^{{{{\rm{S}}}}}{X}_{{{{\rm{L}}}}}^{{{{\rm{T}}}}}$$\end{document} -operator via an auxiliary register (Aux.) and applying a Z-type feedback operation based on the measurement outcome (first green box). The two blocks are then split again by measuring \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Z}_{{{{\rm{L}}}}}^{{{{\rm{S}}}}}$$\end{document} and applying an X-type operation to the target register conditioned on the measurement outcome (second green box). d Schematic illustration of modular state teleportation, where the source and target registers are never directly coupled to one another, but only interact via an auxiliary quantum register (Aux). e High-level circuits for measurement-free logical teleportation and experimental logical quantum state tomography. We replace the measurements and feed-forward operations with coherent feedback operations to teleport a state without mid-circuit measurements (blue). An additional HL is applied to the target state using the circuit shown in orange. The reset operation can either be carried out explicitly by physically resetting the auxiliary qubits and reusing them afterwards, or implemented by replacing them with fresh qubits.
A standard approach for teleporting a state between registers is based on lattice surgery^45,46^, which is illustrated in Fig. 1c. First, two code blocks are merged by measuring the joint logical X-operator. Based on this measurement outcome, one applies a logical Z-operation to the target register. In a second step, the two blocks are split again by measuring the logical Z-operator of the source register and applying a conditional logical X-operation to the target register. The measurement-based approach has been realized experimentally on various platforms^12,20,47–49^.
Instead of performing measurements and conditional operations based on the measurement outcomes, we now map the respective operators to an auxiliary register and apply coherent feedback operations, as illustrated in Fig. 1e in blue. In the first step, we couple both logical qubit registers to the auxiliary register by applying pairs of CNOT-gates to map the information about the joint logical operator \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}_{{{{\rm{L}}}}}^{{{{\rm{S}}}}}{X}_{{{{\rm{L}}}}}^{{{{\rm{T}}}}}$$\end{document} of the source (S) and the target (T) register to the auxiliary qubits. The conditional logical Z-operation can then be implemented coherently with a combination of CZ-gates that act on the auxiliary and target register, as shown in the green dashed box in Fig. 1e. In the second step, we map the logical \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Z}_{{{{\rm{L}}}}}^{{{{\rm{S}}}}}$$\end{document} to the auxiliary register and apply a coherent feedback with a combination of CNOT-gates. The scheme is made FT by repeating subroutines, i.e., by mapping multiple stabilizer-equivalent logical operators onto auxiliary qubits, as discussed further in the ‘Methods’ Section. The explicit circuits can be found in the ‘Methods’ Subsection ‘Circuits’.
We construct a similar protocol that enables the implementation of a logical HL-gate (shown in Fig. 1e in orange) with the same resources as the bare teleportation protocol. We find this circuit by inserting a physical H-gate to the target qubit and propagating it back through the circuit, such that no H-gate has to be performed explicitly. This means that no HL-gate has to be applied to a [[4, 1, 2]] instance when shifting this circuit to the logical level.
We experimentally perform logical state tomography for three protocols: state initialization, logical state teleportation and the application of a logical HL-gate for logical input states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 0\rangle }_{{{{\rm{L}}}}}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+\rangle }_{{{{\rm{L}}}}}$$\end{document} , which is shown in Fig. 1b, e. Further details on the measurement bases and number of shots can be found in the ‘Methods’ Section. We achieve fidelities of up to 93(2)% for state teleportation and 95(3)% for an HL-gate, which is lower than the respective physical gate operation with a fidelity of ≈0.996. The difference in fidelities for the two logical input states can be traced back to two sources. First, dephasing on idling qubits due to fluctuations in the magnetic fields introduces a strong bias towards Z-type errors. Furthermore, we measure the qubits in the Z-basis in the end and determine the logical value from this measurement, if the target state is a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 0\rangle }_{{{{\rm{L}}}}}$$\end{document} -state. Based on these outcomes, we perform a classical round of error detection and postselect on the two Z-stabilizers of the [[4, 1, 2]]-code. When the target state is the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+\rangle }_{{{{\rm{L}}}}}$$\end{document} -state, and we determine the logical X-value, we can only postselect on one X-stabilizer. By accepting fewer runs, we effectively also discard a fraction of runs where higher-weight errors lead to a failure, and fidelities increase the more we postselect.
Our FT logical state teleportation scheme can, in principle, be scaled to higher-distance surface codes, which is discussed further in the ‘Methods’ Section. Here, the key idea is to use d representations of logical operators on a distance-d code to ensure that no weight-d fault leads to a logical failure.
FT toolbox for universal operations on the [[8, 3, 2]]-code
In this section, we discuss circuit constructions for a FT universal gate set on an eight-qubit error-detecting code, which we use to implement Grover’s search algorithm on three logical qubits experimentally in the following section. The [[8, 3, 2]]-code is the smallest instance of a three-dimensional color code^26,27,50–52^ that encodes k = 3 logical qubit in n = 8 physical qubits and has distance d = 2, meaning that any single error can be detected. The X-stabilizer and a Z-stabilizer of this code have support on all eight qubits, while three additional weight-4 Z-stabilizers are defined on three faces of a cube, intersecting on edges, as shown in Fig. 2a. The three logical Pauli X-operators of this code have support on the weight-4 faces of the cube, while the logical Z-operators are defined on edges of weight 2.Fig. 2FT logical operations on an [[8, 3, 2]]-code.a Definition of stabilizers and logical operators on the [[8, 3, 2]]-code^26,27,52^. b The upper left circuit implements measurement-based H-gate injection, where an auxiliary qubit is prepared in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|+\rangle$$\end{document} and entangled with the data qubit in state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| \psi \rangle$$\end{document} . In this protocol, one would measure the auxiliary qubit and apply a Pauli operation that depends on the measurement outcome m. The measurement and conditional operation (green dashed box) can be replaced with a combination of CNOT-gates (lower left), such that no mid-circuit measurements or feed-forward operations are required. We shift this scheme to the logical level by replacing the data qubit with one logical qubit of the [[8, 3, 2]]-code and the auxiliary qubit with one logical qubit of a [[4, 2, 2]]-code, which supports a natively transversal HL-gate, up to a simple relabeling.
The [[8, 3, 2]]-code supports a transversal non-Clifford gate^26,27,52^: the CCZ-gate can be implemented by applying single-qubit T- and T^†^-gates to individual qubits as illustrated in Fig. 3a, such that errors do not propagate within one code block. A logical CNOT-gate between qubits that are encoded within the same encoded block can be implemented by swapping pairs of qubits. In the following, we implement these CNOT-gates within one block by relabeling pairs of qubits, which does not require any physical gate operations. The [[8, 3, 2]]-code has in the past been used for multiple experimental demonstrations^15,26,27^. Recent theoretical works have proposed constructions for measurement-free, FT universal quantum computing^39,40^, but require a substantial overhead in gate operations and qubit count. Here, we introduce an implementation of a HL-gate for the [[8, 3, 2]]-code that does not rely on mid-circuit measurements or feed-forward operations, and, together with the CCZ-gate, completes a FT universal gate set. Our construction for the FT single-qubit logical HL-gate is based on state injection^53^. State injection makes use of a suitable resource state^54^, which is injected onto the data qubit by, first, entangling the two qubits, then measuring the resource qubit and, finally, applying a Clifford operation to the data qubit conditioned on the measurement outcome in the second step. Figure 2b1 shows the circuit that may be used to apply an HL-gate to a state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| \psi \rangle }_{{{{\rm{L}}}}}$$\end{document} by means of state injection. Here, an auxiliary qubit is prepared in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|+\rangle$$\end{document} as a resource state, then entangled with the data qubit with a combination of a CNOT- and a CZ-gate. Finally, the auxiliary qubit is measured in the X-basis, and either a Pauli X- or Z-flip is applied to the data qubit, depending on the measurement outcome m.Fig. 3. Two-solution Grover search on a database of size N = 2^3^.a We recompile the Grover search algorithm that includes a phase oracle^55,62^ into the available FT gate operations within one [[8, 3, 2]]-code {HL, CNOT_L_, CCZL}. b Experimentally obtained (dark) and numerically simulated (light) probabilities for each possible solution. The two solutions 011 and 101 can be clearly distinguished, and the total probability to find one of the two solutions (green) is psuccess = p011 + p101 = 0.40(4).
We now replace the measurement of the auxiliary qubit with a coherent feedback operation comprising two CNOT-gates, as shown in Fig. 2b2. A measurement can always be replaced with a quantum circuit^3,34^, but does not automatically obey fault-tolerant circuit design principles. In our circuit construction, we have to apply H-gates to the auxiliary qubit in order to achieve the desired H-gate injection to the data qubit. If both qubits corresponded to logical qubits of the same code, there would be no benefit in using this approach, because it would require an HL-gate in order to inject one. We therefore use different types of codes to inject the desired gate operation. Specifically, we consider the three encoded qubits of the [[8, 3, 2]]-code and inject an HL-gate onto one of the logical qubits by means of an auxiliary [[4, 2, 2]]-code prepared in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+0\rangle }_{{{{\rm{L}}}}}$$\end{document} as a resource state. In this circuit, we require a CNOT-gate that acts on two logical qubits that are encoded in two different code blocks. This gate can be implemented with a non-transversal, yet FT, gate implementation, implying that any single fault that may propagate through the full circuit remains detectable afterwards. This logical inter-block CNOT-gadget and the full circuit for the injection of an HL-gate are depicted in Suppl. Fig. 4. Our implementation of the FT logical HL-gate requires 4 auxiliary qubits and 26 two-qubit gates.
We perform experimental logical state tomography for each logical qubit, considering FT logical state initialization, the single-logical HL-gate and the transversal CCZL-gate on the [[8, 3, 2]]-code. We achieve fidelities of up to 81(3)% for HL on logical qubit 0, accepting 10% of the runs after postselection. Moreover, we find fidelties between 65(6)% and 99.89(14)% for the two idling logical qubits, depending on the logical input state and its sensitivity to dephasing. All results are shown and further analyzed in the ‘Methods’ Section. Notably, we identify dephasing of idling qubits as a major error source, which we estimate to account for almost two-thirds of the overall logical error rate, as discussed further in ‘Methods’ Subsection ‘Error budget’.
The presented FT universal gate set on the [[8, 3, 2]]-code unlocks the capability to run minimal logical algorithms without relying on explicit mid-circuit measurement or feed-forward operations. In the next step, we use it to implement an FT Grover search on three logical qubits encoded in the [[8, 3, 2]]-code.
Grover search on logical qubits
Grover’s search algorithm^55,56^ enables quantum computers to search through unsorted databases significantly more efficiently than classical methods. It achieves a quadratic speedup by reducing the number of queries required to find a desired item, and can be used as a subroutine for other quantum algorithms^57–60^. Grover’s algorithm consists of three steps^3,61^:
- Initialization: Prepare all qubits in an equal-weight superposition of the computational basis states with the Hadamard transform, i.e., apply single-qubit H-gates to all qubits.
- Grover iteration: Perform (a) and (b) j times to amplify the amplitude of the solution-states s:
- Apply an oracle operator O that marks the solutions by flipping the sign of these states.
- Apply a diffusion operator D that reflects the state about the initial state.
- Measurement: Measure the qubit register in the computational basis.
We implement Grover’s algorithm on three logical qubits, thus searching a database of size N = 2^n^ = 8 bits. As an example, we consider the phase oracle that marks s = 2 solution-states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 011\rangle$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 101\rangle$$\end{document} . In this setting, the probability of finding a solution after one Grover iteration in a noise-free setting is 1, which is discussed further in ‘Methods’ Subsection ‘Grover’s search algorithm’. The optimal classical search corresponds to performing a single query, followed by a random guess, and the probability to find a solution in this case is s/N + (N − s)/N ⋅ s/(N − 1) ≈ 0.46 in our case. Grover’s search algorithm has been implemented on physical qubits on trapped ions^62,63^, superconducting architectures^64–66^, on spin qubits in silicon^67^, and on molecules using NMR techniques^68^. It has also been realized on two logical qubits encoded in a [[4, 2, 2]]-code with better-than-physical logical Clifford operations^16,28^, searching a database of N = 4. This algorithm with N = 4 does not require a universal set of gates, but can be realized with Clifford-gate operation only.
We implement the three-qubit Grover’s algorithm on logical qubits encoded in the [[8, 3, 2]]-code by utilizing the universal FT gate set {HL, CNOT_L_, CCZL}. We recompile the initial circuit^3,55,62^ into the available FT gates introduced in the previous section, as shown in Fig. 3a. We then implement this circuit on our experimental trapped-ion quantum processor, accompanied by numerical simulations according to a multi-parameter noise model specified in the ‘Methods’ Section. Figure 3b shows the determined probabilities for each of the eight possible final states, two of which correspond to the correct solution states as marked in green. The total probability to find a solution using the experimental data is psuccess = p011 + p101 = 0.40(4). This overall probability to find a solution in a single shot is slightly lower than the optimal classical probability of 0.46, as determined above. However, as discussed further in the ‘Methods’ Section, only slight enhancements to the current setup are sufficient to outperform the optimal classical algorithm. Numerical simulations show that reducing, e.g., the two-qubit-gate error rate by 1% to p2 ≈ 0.015, which has been demonstrated on experimental trapped-ion platforms^31,69,70^, leads to an overall success rate of ≈ 0.52, which clearly outperforms the optimal classical strategy. Instead of reducing p2, also extending the coherence time to T2 = 100 ms, which has been shown in independent technical demonstrations^71–75^, leads to a success probability of psuccess ≈ 0.67. This demonstrates that for only slightly smaller error rates on idling qubits and two-qubit gate operations, a regime where the measurement-free quantum algorithm outperforms its classical counterpart is reachable today.
Our scheme for Grover’s algorithm can be scaled to a larger search space, provided enough qubits and sufficiently reliable gate operations are available. One can implement the FT gate set {HL, CNOT_L_, CCZL} on logical qubits encoded within one [[8, 3, 2]] block. In addition, one can apply an inter-block CNOT-gate between two logical qubits of two distinct [[8, 3, 2]]-codes^76^. These operations enable the implementation of an oracle and the amplification on more than three qubits by decomposing the required gates into the available gate sets^3^.
Discussion
In this work, we introduce and experimentally implement a complete toolbox of operations for fault-tolerant (FT) universal quantum computing without mid-circuit measurements. Our work presents the first experimental realization of a FT universal gate set that operates without mid-circuit measurements and marks the FT implementation of Grover’s algorithm on a search space of up to N = 8 on encoded logical qubits, demonstrating for the first time a FT logical algorithm without mid-circuit measurements.
Our schemes are tailored toward trapped-ion architectures that provide all-to-all connectivity^9,31,75,77^, but they can be analogously implemented on other architectures. For example, neutral atom platforms have demonstrated the capabilities required for implementing the presented code constructions^15,78^. These architectures offer long-range connectivity and high-fidelity single- and two-qubit gates, while mid-circuit measurements and real-time feedback are still experimentally demanding due to relatively long measurement times^15,78–80^. These features make our measurement-free implementations ideal candidates for neutral atom platforms, potentially enhancing performance by avoiding costly circuit components, provided that mid-circuit measurements are no noisier than a two-qubit gate. Future developments may depend on improvements in the measurement fidelities.
Future work will include the analysis of our protocols for higher-distance codes, as outlined above, and the investigation of thresholds and required overheads in terms of qubit count and gate operations, including extensions to fault-tolerant realizations under restricted qubit-connectivity^81^. Moreover, we have identified dephasing on idling qubits during two-qubit gates as a major logical error source in our experimental demonstration. Further adjustment of our schemes to a biased noise setting^41,82^, which is often given in experimental architectures^15,17,79^, could therefore potentially boost the performance while reducing overheads.
Our work presents the first demonstration of measurement-free fault-tolerant quantum computation and lays the ground for further exploring the full potential of this new paradigm of fault-tolerant quantum information processing without mid-circuit measurements.
Methods
Qubit reset
The qubit reset procedure allows for a selective reinitialization of some of the qubits to state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 0\rangle$$\end{document} . Physically, the reset is performed by quenching the lifetime of the 3^2^D_5/2_ manifold by illuminating the ion chain with the 854 nm laser with subsequent optical pumping to reinitialize the qubits in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 0\rangle$$\end{document} state. The data qubits are hidden in the 4^2^S_1/2_ manifold \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| {4}^{2}{{{{\rm{S}}}}}_{1/2},{m}_{J}=-1/2\rangle,| {4}^{2}{{{{\rm{S}}}}}_{1/2},{m}_{J}=+ 1/2\rangle$$\end{document} by means of the electron shelving technique^83^ during the lifetime quenching to preserve their state. This re-encoding of the data qubits in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| {4}^{2}{{{{\rm{S}}}}}_{1/2},{m}_{J}=-1/2\rangle,| {4}^{2}{{{{\rm{S}}}}}_{1/2},{m}_{J}=+ 1/2\rangle$$\end{document} levels instead of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| {4}^{2}{{{{\rm{S}}}}}_{1/2},{m}_{J}=-1/2\rangle,| {3}^{2}{{{{\rm{D}}}}}_{5/2},{m}_{J}=-1/2\rangle$$\end{document} levels results in a higher sensitivity to magnetic field noise and, consequently, lower coherence time. Therefore, we perform two dynamical decoupling pulses (DD) with the radio-frequency (RF) antenna driving the transition between the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| {4}^{2}{{{{\rm{S}}}}}_{1/2},{m}_{J}=-1/2\rangle$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| {4}^{2}{{{{\rm{S}}}}}_{1/2},{m}_{J}=+ 1/2\rangle$$\end{document} levels. The sketch of the procedure is shown in Fig. 4a, while additional details can be found in ref. ^7^. The reset procedure does not require recooling of the ion chain, unlike the full mid-circuit measurement, since the reset ions emit only a few photons during the procedure. Consequently, the reset is faster than our current implementation of the mid-circuit measurement (1.7 ms vs. ≈ 30 ms) and the preservation of the data qubit’s state is higher (process fidelity 0.955(9) vs. 0.908(12)). The χ-matrix for the reset procedure obtained via quantum process tomography is depicted in Fig. 4b. The main fidelity limitation for both processes at the moment is the coherence of the ground state qubit ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| {4}^{2}{{{{\rm{S}}}}}_{1/2},{m}_{J}=-1/2\rangle$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| {4}^{2}{{{{\rm{S}}}}}_{1/2},{m}_{J}=+ 1/2\rangle$$\end{document} ) of 5.3(6) ms. The recent hardware modifications, reducing the magnetic field noise, allowed us to extend this coherence time to 113(17) ms, which should significantly improve the fidelity of mid-circuit operations. Moreover, other trapped-ion setups have demonstrated even higher coherence times of the order of seconds, e.g. ^84^. The fidelities can be further boosted by improving the hiding pulses fidelities by the use of composite pulses^85^. The results will be presented in future works.Fig. 4. Experimental details of the qubit reset procedure.a The pulse sequence implementing the qubit reset procedure. Firstly, the data qubits are hidden in the 4^2^S_1/2_ manifold by applying a global RF pulse and individually addressing each data qubit with a π-pulse. After that, the auxiliary qubits are reinitialized to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 0\rangle$$\end{document} . Finally, the data qubits' encoding is restored to the original. The dynamical decoupling (DD) pulses are inserted to mitigate decoherence of the data qubits in the 4^2^S_1/2_ manifold. b The χ-matrix representation of the reset procedure as the process acting on data qubits, averaged over all data qubits. The area and the color coding of the squares correspond to the absolute value and the phase of an element of the χ-matrix, respectively. The dashed square represents an ideal outcome, specifically the identity process.
We make use of the full 16-ion register and use fresh auxiliary qubits as long as possible. The reset procedure is only used in our implementation of Grover’s algorithm (see Fig. 3) to reset one auxiliary qubit that is used for the FT preparation of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+00\rangle }_{{{{\rm{L}}}}}$$\end{document} for the [[8, 3, 2]]-code (see Suppl. Fig. 3e). This auxiliary qubit is later used for the mapping of one Z-stabilizer of the [[8, 3, 2]]-code, as discussed in ‘Methods’ Subsection ‘Tomography’.
Anticipated performance of measurement-free state teleportation
The measurement-free logical teleportation schemes are made FT as illustrated in Fig. 5a. First, we prepare auxiliary two-qubit GHZ-states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| {\psi }_{{{{\rm{aux}}}}}\rangle=(| 00\rangle+| 11\rangle )/\sqrt{2}$$\end{document} stabilized by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{X}^{{{{\rm{aux}}}}}={X}_{8}{X}_{9}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{Z}^{{{{\rm{aux}}}}}={Z}_{8}{Z}_{9}$$\end{document} , which ensures that no single fault on an auxiliary qubit propagates to a logical error when the two registers are coupled. In addition, we map two representations of the joint logical operator \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}_{{{{\rm{L}}}}}^{{{{\rm{S}}}}}{X}_{{{{\rm{L}}}}}^{{{{\rm{T}}}}}$$\end{document} , that have fully disjoint support, onto the auxiliary register such that no single fault on a data qubit leads to a logical error on the output state (panels 1 and 2 in Fig. 5a). Here, the information about each representation of the joint logical operator is stored in one physical auxiliary qubit, which then acts as a control qubit in the coherent feedback operation consisting of CZ-gates. The same strategy is used in step 3, where two representations of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Z}_{{{{\rm{L}}}}}^{{{{\rm{S}}}}}$$\end{document} with disjoint support are mapped onto two physical qubits. The explicit circuits can be found in the ‘Methods’ Subsection ‘Circuits’. The non-FT (nFT) counterparts of these protocols make use of a bare physical auxiliary qubit and only map a single representation of the respective operators onto this auxiliary qubit, which is then used to control the coherent feedback operation.Fig. 5. Fault tolerance and scaling to higher code distances for measurement-free logical state teleportation.a We prepare an auxiliary two-qubit GHZ-state to prevent single faults on auxiliary qubits from causing a logical failure. In addition, we map two stabilizer-equivalent representations of the joint logical operator with fully disjoint qubit support onto the auxiliary registers (1, 2). The same strategy is used in step 3, where two equivalent but fully disjoint logical Z-operators of the source register are mapped onto two physical qubits. b Numerically determined scaling of the logical infidelity for FT and non-FT logical teleportation and the HL-gate, averaged over initial states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 0\rangle }_{{\mathrm{L}}}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+\rangle }_{{\mathrm{L}}}$$\end{document} . We fix the error parameters \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\vec{p}=({p}_{{{{\rm{1}}}}},{p}_{{{{\rm{2}}}}},{p}_{{{{\rm{m}}}}},{p}_{{{{\rm{i}}}}},{p}_{{{{\rm{idle}}}}})$$\end{document} to experimental error rates^7,17^ and scale these with a common improvement factor λ. We identify a quadratic scaling of the infidelity for the FT protocols with λ, which indicates—as expected—that no single fault leads to a logical failure. The FT teleportation protocol outperforms its non-FT counterpart already for the current experimental noise parameters (λ = 1). The inset shows the logical infidelities at λ = 1 obtained from the experiment (darker color) and numerical simulations (lighter color). c Scaling measurement-free state teleportation to surface codes with higher distances d > 2. Each string of qubits connecting opposing boundaries supports a representation of a logical \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}_{{{{\rm{L}}}}}^{{{{\rm{S}}}}}$$\end{document} (upper lattice) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}_{{{{\rm{L}}}}}^{{{{\rm{T}}}}}$$\end{document} (lower lattice); one exemplary representation is shown in red. There are d equivalent representations that have fully disjoint support. Each one can be mapped onto an auxiliary d-qubit GHZ-state, and coherent feedback steps can be applied, which are controlled by the state on d physical auxiliary qubits.
In Fig. 5b, we simulate the scaling of the logical infidelity for FT and non-FT measurement-free logical state teleportation and the application of the HL-gate operation by means of teleportation. Here, we consider a multi-parameter noise model, attributing different error rates to each type of component in the circuits. Specifically, we consider depolarizing noise on single-qubit gates with a probability p1 = 3.6 × 10^−3^, two-qubit depolarizing noise on two-qubit gates with a probability p2 = 2.5 × 10^−2^, flipped physical qubit initializations with a probability pi = 3 × 10^−3^, and flips before the final projective measurements with a probability pm = 3 × 10^−3^. The values of the error rates correspond to the ones in our experimental setup^7,17,21,42^. We implement dephasing on all idling qubits, where a Z-fault is applied to each idling physical qubit with a probability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${p}_{{{{\rm{idle}}}}}=(1-{e}^{-{t}_{{{{\rm{gate}}}}}/{T}_{2}})$$\end{document} given the gate time of the respective operation and the coherence time T2 = 50 ms. Note that gates in the experimental setup can only be executed sequentially, which increases the total dephasing time. The noise channels and numerical methods are explicitly given in ‘Methods’ Subsection ‘Numerical methods’. We scale the error parameters \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\vec{p}(\lambda )=\lambda \cdot ({p}_{1},{p}_{2},{p}_{i},{p}_{m},{p}_{{{{\rm{idle}}}}})$$\end{document} with a common factor λ, such that λ = 1 corresponds to the set of parameters as given in the current experimental setup. As expected, the FT protocols scale quadratically with λ, indicating that the required fault-tolerance properties are fulfilled.
Our approach for logical state teleportation without mid-circuit measurements is, in principle, straightforwardly scalable to higher distance d > 2 surface codes, as illustrated in Fig. 5c. These codes have d equivalent representations of the logical Pauli operators which do not share support. For the FT mapping of the weight-d logical operators, we then have to prepare the auxiliary register fault-tolerantly in a d-qubit GHZ-state, and apply coherent feedback steps controlled on the state of d physical auxiliary qubits. However, this approach requires the FT preparation of potentially high-weight GHZ-states. An alternative route to scalability could entail concatenating the presented schemes to avoid the need for high-weight physical GHZ-states. Supplementary Fig. 1 shows the averaged simulated logical error rates for the [[4, 1, 2]]-code. For simplicity, here we implement a single-parameter noise model where each two-qubit gate introduces an error with probability p. Additionally, we count the lowest-weight fault-locations that lead to a logical failure in the [[4, 1, 2]]-code. Specifically, 1792 weight-two fault-locations lead to a failure for the initial state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+\rangle }_{L}$$\end{document} and 2784 fault-locations for the initial state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 0\rangle }_{L}$$\end{document} . Note that for each two-qubit gate, we consider 15 possible locations corresponding to the 15 distinct two-qubit Pauli errors. These counts provide an estimate of the logical failure rate, which scales as pL ≈ c × p^2^. Given this polynomial, we can infer the failure rates for higher levels of concatenation because \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${p}_{L}^{(n)}={c}^{{2}^{n}-1}\cdot {p}^{{2}^{n}}$$\end{document} , i.e., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${p}_{L}^{(1)}=c\cdot {p}^{2}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${p}_{L}^{(2)}=c\cdot {p}_{L}^{(1)2}={c}^{3}{p}^{4}$$\end{document} and so on. As shown in Fig. 2, the calculated polynomials (dashed lines) for the first two concatenation levels agree with the numerically simulated ones (solid curves). These results demonstrate that our approach for measurement-free logical state teleportation can be scaled to higher distances without the need for large GHZ-states. However, within the error-detecting framework, we always have to discard a fraction of runs, and this fraction increases with more layers of concatenation, as can be seen in Suppl. Fig. 1b, because more and higher-weight errors can be detected. Another approach to achieving scalability is to rely on different auxiliary states that still maintain fault tolerance in the presented protocols. Preliminary results indicate that MF FT operation of distance-3 surface code by means of carefully crafted quantum circuits without the need of FT prepared GHZ ancillary states is also viable, but a detailed analysis and generalization to higher code distances is left for future follow-up work.
Numerical methods
We use Monte-Carlo simulations to estimate the logical infidelities of our protocols^86^. Each circuit component is modeled by first applying the respective ideal operation, followed by an error E occurring with probability p. We simulate a depolarizing noise channel after every single- and two-qubit gate. With probabilities p1 and p2, an error from the respective sets is applied. These probabilities define the corresponding error channels
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{rcl}{{{{\mathcal{E}}}}}_{1}(\rho ) &=& (1-{p}_{1})\rho+\frac{{p}_{1}}{3}{\sum }_{i=1}^{3}{E}_{1}^{i}\rho {E}_{1}^{i}\\ {{{{\mathcal{E}}}}}_{2}(\rho ) &=& (1-{p}_{2})\rho+\frac{{p}_{2}}{15}{\sum }_{i=1}^{15}{E}_{2}^{i}\,\rho \,{E}_{2}^{i}.\end{array}$$\end{document}with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E}_{1}^{k}\in \{X,Y,Z\}$$\end{document} for k = 1, 2, 3 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E}_{2}^{k}\in \{IX,XI,XX,IY,YI,YY,IZ,ZI,ZZ,XY,YX,XZ,ZX,YZ,ZY\}$$\end{document} for k = 1, …, 15. All qubits are initialized and measured in the Z-basis, at the very end of the respective protocols. To simulate faults in these operations, we apply X-flips after initialization and before measurement, each occurring with probabilities pinit and pmeas, respectively. Moreover, qubits that remain idle during gate operations may experience dephasing, which we model with the error channel
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{{\mathcal{E}}}}}_{{{{\rm{idle}}}}}(\rho )=(1-{p}_{{{{\rm{idle}}}}})\rho+{p}_{{{{\rm{idle}}}}}Z\rho Z.$$\end{document}The probability pidle depends on the execution time t of the performed gate and the qubit coherence time T2 = 50 ms
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${p}_{{{{\rm{idle}}}}}=\frac{1}{2}\left[1-\exp \left(-\frac{t}{{T}_{2}}\right)\right].$$\end{document}In our simulations, we use t1 = 70 μs as the gate time of single-qubit gates and t2 = 350 μs as the gate time of two-qubit gates, as summarized in Suppl. Tab. 1.
We measure the final state in the logical X-, Y-, and Z-basis for each protocol, as described in ‘Methods’ Subsection ‘Tomography’, and calculate the state fidelity between the ideal logical state ρ1 and the reconstructed density operator ρ2 obtained after postselecting, as
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F({\rho }_{1},{\rho }_{2})={{{\rm{Tr}}}}{\left[\sqrt{\sqrt{{\rho }_{1}}{\rho }_{2}\sqrt{{\rho }_{1}}}\right]}^{2}.$$\end{document}We use Qiskit’s Quantum Information package to calculate fidelities^87^.
Circuits
Supplementary Fig. 2 shows the explicit circuits that were implemented for FT logical state teleportation discussed in Section ‘Logical state teleportation without mid-circuit measurement’. Supplementary Fig. 3 shows the circuit constructions for the FT logical state initializations on the [[4, 1, 2]]- and the [[8, 3, 2]]-code without mid-circuit measurements as implemented in the demonstrated protocols. In these circuits, fault tolerance is maintained even without measurements by means of a flag-qubit-controlled reduction of potentially dangerous weight-2 errors to weight-1 configurations. Supplementary Fig. 4 depicts the circuit construction for the application of a single-logical HL-gate on the [[8, 3, 2]]-code.
Tomography
[[4, 1, 2]]-code
We perform logical state tomography for two logical input states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 0\rangle }_{{{{\rm{L}}}}}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+\rangle }_{{{{\rm{L}}}}}$$\end{document} considering logical state preparation, state teleportation and the application of a HL-gate on the [[4, 1, 2]]-code, as shown in Fig. 1. To this end, we measure in the X-, Y-, and Z-basis to extract the respective logical expectation values. For measurements in the X-basis, we measure all physical qubits in the X-basis in the end and infer the logical X-operator and the X-type stabilizer from this measurement. Analogously, we can extract the logical Z-operator and the Z-type stabilizers for measurements of all physical qubits in the Z-basis. However, we cannot simply determine the required stabilizers and the logical value at the same time for measurements in the Y-basis, because they share support but are of different Pauli-type, as for example YL = Y0X1Z2 and the Y-type stabilizer SY = Y0Y1Y2Y3. We therefore map out the Y-stabilizer onto a physical auxiliary qubit with the circuit shown in Suppl. Fig. 6. Here, the gate ordering ensures that no hook error, i.e., a fault on an auxiliary qubit that may propagate onto multiple data qubits, leads to a logical flip in the subsequent measurement, as any single propagated fault is still detected in the end by a Z-stabilizer. We finally measure the 4 qubits in the Y-, X-, and Z-basis, allowing us to determine YL and one additional Z-stabilizer \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{Z}^{2}={Z}_{2}{Z}_{3}$$\end{document} .
[[8, 3, 2]]-code
We perform logical state tomography for each state prepared with the specified protocol, i.e., logical state preparation and the logical operations HL and CCZL on [[8, 3, 2]] for different input states, as shown in Fig. 6. We consider each individual logical qubit and perform tomography for each one independently. For measurements in the X-basis, we measure all physical qubits in the X-basis and determine the three logical Pauli-operators \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}_{{{{\rm{L}}}}}^{0}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}_{{{{\rm{L}}}}}^{1}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}_{{{{\rm{L}}}}}^{2}$$\end{document} and the X-stabilizer SX, as defined in Fig. 2a. Analogously, we extract the logical Z-operators simultaneously, along with the Z-stabilizers, when measuring all physical qubits in the Z-basis. In this case, we additionally map out the X-stabilizer SX onto an auxiliary qubit when performing the HL-gate to achieve fault tolerance. The circuit that is used for measurements in the Z-basis is shown in Suppl. Fig. 5a. In this circuit, a single fault may propagate as illustrated in red, but is detected by the Z-stabilizers afterwards.Fig. 6. Fidelities for different logical input states and non-trivial logical operations.Experimentally obtained logical state fidelities for state initialization, a single-logical HL-gate and the transversal CCZL-gate for the [[8, 3, 2]] code.
For measurements in the Y-basis, we have to take into account that the different logical operators may share support but are of different Pauli-type, as for example \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Y}_{{{{\rm{L}}}}}^{0}={Y}_{0}{X}_{1}{Z}_{2}{X}_{4}{X}_{5}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Y}_{{{{\rm{L}}}}}^{1}={Y}_{0}{X}_{1}{X}_{2}{X}_{3}{Z}_{4}$$\end{document} , so they cannot be extracted simultaneously in a single measurement. We therefore perform three sets of independent experiments and determine \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Y}_{{{{\rm{L}}}}}^{0}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Y}_{{{{\rm{L}}}}}^{1}$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Y}_{{{{\rm{L}}}}}^{2}$$\end{document} individually. For measurements in the Y-basis, we also map the Y-stabilizer SY = Y0Y1Y2Y4Y3Y5Y6Y7 onto an auxiliary qubit when performing the HL-gate. We then measure the physical qubits in different bases to extract the respective logical Y-operator and one additional Z-stabilizer. An exemplary circuit that is used for measurement in the Y-basis for the extraction of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Y}_{{{{\rm{L}}}}}^{0}$$\end{document} is shown in Suppl. Fig. 5b. When extracting \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Y}_{{{{\rm{L}}}}}^{1}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Y}_{{{{\rm{L}}}}}^{2}$$\end{document} , we measure the physical qubits in the bases Y0X1X2X3Z4Z5Z6Z7 and Y0Z1X2Z3X4Z5X6Z7, respectively.
Moreover, the logical auxiliary qubit is still intact after performing the single-logical HL-gate. We also projectively measure the logical auxiliary qubit, extract the stabilizers of the [[4, 2, 2]]-instance, and postselect for a trivial syndrome to increase the fidelities in our protocols. Here, we map the Z-stabilizer of the logical auxiliary qubit onto another physical auxiliary with the circuit shown in Suppl. Fig. 5c.
When we run the full logical Grover search algorithm on the three qubits of the [[8, 3, 2]]-code, we additionally map two Z-stabilizers \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{Z}^{1}={Z}_{0}{Z}_{1}{Z}_{2}{Z}_{3}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{Z}^{2}{S}_{Z}^{3}={Z}_{1}{Z}_{2}{Z}_{5}{Z}_{6}$$\end{document} onto physical auxiliary qubits in the end in order to maintain fault tolerance.
Number of measurements
In the tomography experiments described in ‘Methods’ Subsection ‘Tomography’, each logical state was measured in three measurement bases {X, Y, Z} with the same number of measurements for each basis. The teleportation experiment with the [[4, 1, 2]]-code (see Fig. 1) took 40000 shots for each logical state per measurement basis. The initialization and logical operations with the [[8, 3, 2]]-code (see Fig. 6) took 7500 shots for each logical state and logical qubit per measurement basis. The Grover’s algorithm demonstration (see Fig. 3) took 37,500 shots per measurement basis. All data sets were split into 12 equal subsets; the tomography was performed for every subset, yielding 12 values for the fidelity for every experiment. The final fidelity numbers are the mean and the standard deviation of these 12 values.
Performance of the FT universal gate set on the [[8, 3, 2]]-code
Figure 6 shows the logical state fidelities that were obtained experimentally for FT logical state initialization, the single-logical HL-gate and the transversal CCZL-gate on the [[8, 3, 2]]-code. We find that fidelities are higher if the final target state is a Z-eigenstate, as opposed to an X-eigenstate, due to dephasing, which does not affect the fidelity for Z-eigenstates. Additionally, postselection based on the four Z-stabilizers is more selective than only a single X-stabilizer, which boosts the fidelities in these cases. The degree of postselection is reflected in the acceptance rates: the average acceptance rates in the experiment[simulation] after the state initialization are 0.6[0.64] in the X-basis, 0.48[0.53] in the Y-basis, and 0.3[0.48] for measurements in the Z-basis; the numbers in brackets indicate the acceptance rate obtained in the simulation. After the injection of an HL-gate, these are 0.3[0.2] for measurements in the X-basis, 0.2[0.13] in the Y-basis and 0.1[0.07] in the Z-basis. The fidelities for the state initialization of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+00\rangle }_{{{{\rm{L}}}}}$$\end{document} and CC \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Z}_{{{{\rm{L}}}}}{|+00\rangle }_{{{{\rm{L}}}}}$$\end{document} agree with each other within the given uncertainty interval, since the CCZL consists entirely of virtual Z-rotations, thus no additional operations are physically applied to the qubits.
Figure 7 shows the simulated scaling of the logical infidelity for the logical HL-gate on the [[8, 3, 2]]-code for each logical qubit. We scale the noise parameters \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\vec{p}(\lambda )=\lambda \cdot ({p}_{1},{p}_{2},{p}_{i},{p}_{m},{p}_{{{{\rm{idle}}}}})$$\end{document} given the same values as specified above, such that λ = 1 corresponds to the set of parameters as given for the current experimental setup. The inset shows the state fidelities for the different logical input states obtained from experiment and simulation.Fig. 7. Scaling of the logical infidelity for different logical qubits during the single logical HL-gate.We fix the error parameters \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\vec{p}=({p}_{{{{\rm{1}}}}},{p}_{{{{\rm{2}}}}},{p}_{{{{\rm{m}}}}},{p}_{{{{\rm{i}}}}},{p}_{{{{\rm{idle}}}}})$$\end{document} and scale these with λ. λ = 1 corresponds to the current values in the experimental setup. We determine numerically the scaling of the logical infidelity for the logical qubit 0 (a), on which the HL is applied, and idling logical qubits 1 (b) and 2 (c). The inset shows the logical infidelities at λ = 1 obtained from the experiment (darker color) and from numerical simulations (lighter color). We find that experimentally and numerically obtained state fidelities agree within their uncertainty intervals for qubits 1 and 2, while these values differ on qubit 0 for certain states, as discussed further in ‘Methods’ Subsection ‘Performance of the FT universal gate set on the [[8, 3, 2]]-code’.
We find that the numerically and experimentally obtained fidelities agree for the two idling logical qubits, while the fidelities of the first logical qubit obtained from simulation, shown in blue in Fig. 7a, differ from the experimental result by more than 14% for logical states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+00\rangle }_{{{{\rm{L}}}}}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 000\rangle }_{{{{\rm{L}}}}}$$\end{document} . We attribute this deviation to global dephasing effects due to random fluctuations in the effective magnetic field that act on all physical qubits simultaneously^88^, instead of locally and uncorrelated on each individual qubit. The effect of this global dephasing on the eight-qubit state can be estimated by considering the explicit basis states, for example,
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 000\rangle }_{{{{\rm{L}}}}}=\frac{1}{\sqrt{2}}(| 00000000\rangle+| 11111111\rangle ).$$\end{document}Local dephasing on this state leads to decay of the off-diagonal elements of the density matrix with a factor of e^−Δn/2⋅γ**t^, where γ is a decay constant, and t is time. Δn is the number of positions in the basis states, where the entries of two basis states differ, and corresponds to the Hamming distance. For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 000\rangle }_{{{{\rm{L}}}}}$$\end{document} , Δn = 8 and the decay factor is given by e^−4γ**t^. Global dephasing, on the other hand, will cause the off-diagonal elements to decay with a factor of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${e}^{-{(\Delta m/2)}^{2}/2\cdot \gamma t}$$\end{document} ^88^. Δm is the difference in magnetization of the basis states, where the magnetization of a state is given by the difference between the number of qubits in the ground state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 0\rangle$$\end{document} and the remaining number of bits in the excited state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 1\rangle$$\end{document} ^88^.
For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 000\rangle }_{{{{\rm{L}}}}}$$\end{document} , Δm = 16 and, thus, this prefactor is given by e^−32γ**t^. This means that, for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 000\rangle }_{{{{\rm{L}}}}}$$\end{document} , the off-diagonal elements decay eight times faster for global dephasing than for local. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 000\rangle }_{{{{\rm{L}}}}}$$\end{document} is most sensitive to this global effect, since it is an eight-qubit GHZ-state with maximal difference in the magnetization between its basis states. This effect is expected and found to be less pronounced for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+00\rangle }_{{{{\rm{L}}}}}=\frac{1}{2}(| 00000000\rangle+| 11111111\rangle+| 11001100\rangle+| 00110011\rangle )$$\end{document} , where only some coherences decay according to Δm = 8 and some with Δm = 4, so twice and eight times faster than for local dephasing. We only account for local dephasing in our simulations, as characterized in ‘Methods’ Subsection ‘Numerical methods’, which may partly explain the observed differences between the numerically determined and the experimental fidelities. The deviation between experiment and simulation is less pronounced for the two idling logical qubits, because, considering the states that are most sensitive to this global dephasing effect, i.e., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 000\rangle }_{L}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+00\rangle }_{L}$$\end{document} , these qubits stay in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 0\rangle$$\end{document} -state throughout the whole protocol and are therefore not affected by Z-type errors.
Error budget
Logical errors can be correlated in quantum codes that encode multiple logical qubits, such as block codes^89^ or quantum low-density parity-check codes^90,91^. We investigate these correlated errors by determining the probabilities for each logical error configuration, including single and correlated errors, for the non-transversal single-logical HL-gate on the [[8, 3, 2]]-code. To this end, we prepare logical state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|+00\rangle }_{{{{\rm{L}}}}}$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| 0++ \rangle }_{{{{\rm{L}}}}}$$\end{document} ), then apply the HL-gate to the first qubit and measure destructively in the Z(X)-basis. From this, we infer if one, two or all three logical qubits have been flipped, which corresponds to the probability of logical X(Z)-errors. Figure 8a, b shows the probabilities for logical X- and Z-error configurations on the experimental setup. Notably, logical errors do not occur independently as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p({X}_{{{{\rm{L}}}}}^{i}{X}_{{{{\rm{L}}}}}^{j})\ne p({X}_{{{{\rm{L}}}}}^{i})p({X}_{{{{\rm{L}}}}}^{j})$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p({Z}_{{{{\rm{L}}}}}^{i}{Z}_{{{{\rm{L}}}}}^{j})\ne p({Z}_{{{{\rm{L}}}}}^{i})p({Z}_{{{{\rm{L}}}}}^{j})$$\end{document} , as theoretically predicted in previous works on quantum LDPC codes^92^. Figure 8c, d shows numerical data for a setting without dephasing on idling qubits and for perfectly initialized logical states, to isolate the contribution of the HL-gate protocol. Logical error probabilities decrease substantially, while the overall distribution is maintained. Notably, we find that dephasing attributes for a large part of the overall logical error rate: without dephasing, the logical Z-error rate on qubit 0 drops from almost 0.3 (left-most light blue column in Fig. 8b) to less than 0.1 (left-most orange column in Fig. 8d).Fig. 8. Error budget for HL-gate injection.a, b Simulated and experimentally obtained probabilities for logical X- and Z-error configurations. The numbers on top of the bars correspond to the ratio between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p({X}_{{{{\rm{L}}}}}^{i}{X}_{{{{\rm{L}}}}}^{j})$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p({X}_{{{{\rm{L}}}}}^{i})p({X}_{{{{\rm{L}}}}}^{j})$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p({Z}_{{{{\rm{L}}}}}^{i}{Z}_{{{{\rm{L}}}}}^{j})$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p({Z}_{{{{\rm{L}}}}}^{i})p({Z}_{{{{\rm{L}}}}}^{j})$$\end{document} . For example, the determined probability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p({X}_{{{{\rm{L}}}}}^{0}{X}_{{{{\rm{L}}}}}^{1})$$\end{document} is 2 times larger than the probability one would expect from independent errors on logical qubits 0 and 1, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p({X}_{{{{\rm{L}}}}}^{0})p({X}_{{{{\rm{L}}}}}^{1})$$\end{document} . c, d Probabilities for each logical X- and Z-error configuration without dephasing, and for a perfectly initialized input state. Here, the probabilities for X-error configurations are shown on a logarithmic scale for visibility.
Grover’s search algorithm
The number of required Grover iterations n providing the highest amplification of the solution-states depends on the size of the search space N (N = 2^3^ = 8 in our case) and the number of solutions s. In this work, we use a phase oracle^61,62^ with two solutions (s = 2) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$w\in \{| 011\rangle,| 101\rangle \}$$\end{document} : states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 011\rangle$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| 101\rangle$$\end{document} are marked by the oracle of the form
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$O={{{{\rm{C}}}}}_{{{{\rm{1}}}}}{Z}_{2}\cdot {{{{\rm{C}}}}}_{{{{\rm{0}}}}}{Z}_{2}.$$\end{document}The initial equal-superposition state can be represented as a superposition of solutions and non-solution states^3^
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{rcl}| \psi \rangle &=& \frac{1}{\sqrt{8}}{\sum }_{{\psi }^{{\prime\prime} }\notin \{011,101\}}| {\psi }^{{\prime\prime} }\rangle+\frac{1}{\sqrt{8}}(| 011\rangle+| 101\rangle )\\ &=& \sqrt{\frac{N-s}{N}}| {\psi }^{{\prime} }\rangle+\sqrt{\frac{s}{N}}| w\rangle=\cos \theta | {\psi }^{{\prime} }\rangle+\sin \theta | w\rangle \end{array}$$\end{document}with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s/N}=\sin \theta$$\end{document} , i.e., θ = π/6 in our case. The probability of obtaining a valid solution w when measuring in the computational basis is s/N = 1/4, and the probability of obtaining an orthogonal non-solution state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\psi }^{{\prime} }$$\end{document} equals (N − s)/N = 3/4. One Grover step is a product of two reflections, first about the solution states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| w\rangle$$\end{document} with the oracle O and then about the initial state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$| \psi \rangle$$\end{document} with the diffusion operator D. This corresponds to an overall rotation of the initial state, and the rotation angle can be identified to be 2θ = π/3 in our case. A single application of the Grover iteration, including the oracle O and the diffusion operator D, amplifies the probability of success to 1^3,61^, since
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D\cdot O| \psi \rangle=\cos ((2+1)\theta )\cdot | {\psi }^{{\prime} }\rangle+\sin ((2+1)\theta )\cdot | w\rangle=| w\rangle,$$\end{document}meaning that a solution in the fault-free case is found with certainty. Analogously, the probability of finding a solution after k Grover iterations in a noise-free setting is given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\sin }^{2}((2k+1)\theta )$$\end{document} .
A quantum circuit implementing this algorithm is shown in Fig. 3a. This original circuit can be simplified to allow for a simpler implementation with logical qubits compiled into available logical gates of the [[8, 3, 2]]-code, reducing the number of required HL-gates to one. All operations in the resulting circuit can be fault-tolerantly implemented within the [[8, 3, 2]]-code as described in Section ‘FT toolbox for universal operations on the[[8,3,2]]-code’.
Projected performance of Grover’s algorithm
We simulate Grover’s algorithm on logical qubits for different sets of noise parameters in order to estimate how much physical error rates have to improve to gain an advantage over the classically optimal success probability of 0.46. Figure 9 shows the simulated probabilities to find each possible solution state for the initial set of noise parameters (blue), for a two-qubit-gate error rate reduced by 1% and for an increased coherence time \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{2}^{{\prime} }=2{T}_{2}=100\,\,{{{\rm{ms}}}}$$\end{document} . Both projected scenarios outperform the classical counterpart, indicating that even minor enhancements to the current setup could push performance beyond this break-even point.Fig. 9. Anticipated performance of the two-solution Grover search on logical qubits.We simulate Grover’s search algorithm for the set of noise parameters as characterized by the current experimental setup (blue), indicating a success probability of psuccess = 0.40(4). For a slightly lower two-qubit-gate error rate of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${p}_{2}^{{\prime} }={p}_{2}-0.01=0.015$$\end{document} , we already obtain a total success probability of 0.52(1), which is above the classical optimal success probability of 0.46, as discussed in Section ‘Grover search on logical qubits’. If instead of lowering p2, we increase T2 by a factor of 2–100 ms (orange), we find even higher success rates of psuccess = 0.67(1).
Grover search on physical qubits
We implement Grover’s search algorithm on physical qubits, as compiled in Fig. 3a on our experimental trapped-ion setup, accompanied by numerical simulations; the results are shown in Fig. 10. The total experimental[simulated] success probability of 76(2)%[77(1)%] is larger than for the FT implementation on logical qubits, indicating that the quantum algorithm executed on logical qubits is currently still operated above the break-even point with its counterpart realization on physical qubits.Fig. 10. Probabilities for the two-solution Grover search on physical qubits.We implement the circuit shown in Fig. 3a on physical qubits and determine the probabilities for each outcome in the experiment (darker columns) and simulation (lighter columns).
Supplementary information
Supplementary Information Transparent Peer Review file
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Gottesman, D. Stabilizer Codes and Quantum Error Correction. Ph D thesis (California Institute of Technology, 1997).
- 2Nielsen, M. A. & Chuang, I. L. Quantum Computation and Quantum Information, 10th edn (Cambridge University Press, 2010).
- 3Reichardt, B. W. et al. Demonstration of quantum computation and error correction with a tesseract code. Preprint at ar Xiv 10.48550/ar Xiv.2409.04628 (2024).
- 4Google Quantum AI and Collaborators. Quantum error correction below the surface code threshold. Nature 638, 920–926 (2025).10.1038/s 41586-024-08449-y PMC 1186496639653125 · doi ↗ · pubmed ↗
- 5Lacroix, N. et al. Scaling and logic in the color code on a superconducting quantum processor. Nature 645, 614–619 (2025).10.1038/s 41586-025-09061-4PMC 1244361140418964 · doi ↗ · pubmed ↗
- 6Dasu, S. et al. Breaking even with magic: demonstration of a high-fidelity logical non-Clifford gate. Preprint at ar Xiv 10.48550/ar Xiv.2506.14688 (2025).
- 7Bernstein, E. & Vazirani, U. Quantum complexity theory. In Proc. 25th Annu. ACM Symp. Theory Comput. 11–20 (1993).
- 8Reichardt, B. W. et al. Logical computation demonstrated with a neutral atom quantum processor. Preprint at ar Xiv 10.48550/ar Xiv.2411.11822 (2024).