One of the main problems in computational physics is predicting the low-energy behavior of many-body quantum systems. The computational complexity of this problem, however, is relatively poorly understood. A recent major progress in this direction has been the no low-energy trivial states (NLTS) theorem; it gives a family of qubit Hamiltonians whose low-energy states cannot be reached by shallow quantum circuits. In this work we provide a fermionic counterpart to this theorem, constructing local fermionic Hamiltonians with no low-energy trivial states. Distinct from the qubit case, we define trivial states via finite-depth fermionic quantum circuits. We further strengthen the result, allowing free access to (generally, deep) Gaussian fermionic circuits into our notion of a trivial state. The desired fermionic Hamiltonian can be constructed using any qubit Hamiltonian which has the NLTS property via well-spread distributions over bitstrings. We also define a fermionic analog of quantum probabilistically checkable proofs (PCPs) and explore the relation of fermionic PCP class with the qubit version.

Minimizing leakage from computational states is a challenge when using many-level systems like superconducting quantum circuits as qubits. We realize and extend the quantum-hardware-efficient, all-microwave leakage reduction unit (LRU) for transmons in a circuit QED architecture proposed by Battistel et al. This LRU effectively reduces leakage in the second- and third-excited transmon states with up to 99% efficacy in 220 ns, with minimum impact on the qubit subspace. As a first application in the context of quantum error correction, we show how multiple simultaneous LRUs can reduce the error detection rate and suppress leakage buildup within 1% in data and ancilla qubits over 50 cycles of a weight-2 stabilizer measurement.

We consider the problem of approximating the ground state energy of a fermionic Hamiltonian using a Gaussian state. In sharp contrast to the dense case [1, 2], we prove that strictly q-local sparse fermionic Hamiltonians have a constant Gaussian approximation ratio; the result holds for any connectivity and interaction strengths. Sparsity means that each fermion participates in a bounded number of interactions, and strictly q-local means that each term involves exactly q fermionic (Majorana) operators. We extend our proof to give a constant Gaussian approximation ratio for sparse fermionic Hamiltonians with both quartic and quadratic terms. With additional work, we also prove a constant Gaussian approximation ratio for the so-called sparse SYK model with strictly 4-local interactions (sparse SYK-4 model). In each setting we show that the Gaussian state can be efficiently determined. Finally, we prove that the O(n⁻¹/²) Gaussian approximation ratio for the normal (dense) SYK-4 model extends to SYK-q for even q > 4, with an approximation ratio of O(n¹/^2−q/4). Our results identify non-sparseness as the prime reason that the SYK-4 model can fail to have a constant approximation ratio [1, 2].

We consider the task of spectral estimation of local quantum Hamiltonians. The spectral estimation is performed by estimating the oscillation frequencies or decay rates of signals representing the time evolution of states. We present a classical Monte Carlo (MC) scheme which efficiently estimates an imaginary-time, decaying signal for stoquastic (i.e. sign-problem-free) local Hamiltonians. The decay rates in this signal correspond to Hamiltonian eigenvalues (with associated eigenstates present in an input state) and can be extracted using a classical signal processing method like ESPRIT. We compare the efficiency of this MC scheme to its quantum counterpart in which one extracts eigenvalues of a general local Hamiltonian from a real-time, oscillatory signal obtained through quantum phase estimation circuits, again using the ESPRIT method. We prove that the ESPRIT method can resolve S = poly(n) eigenvalues, assuming a 1/poly(n) gap between them, with poly(n) quantum and classical effort through the quantum phase estimation (QPE) circuits, assuming efficient preparation of the input state. We prove that our MC scheme plus the ESPRIT method can resolve S = O(1) eigenvalues, assuming a 1/poly(n) gap between them, with poly(n) purely classical effort for stoquastic Hamiltonians, requiring some access structure to the input state. However, we also show that under these assumptions, i.e. S = O(1) eigenvalues, assuming a 1/poly(n) gap between them and some access structure to the input state, one can achieve this with poly(n) purely classical effort for general local Hamiltonians. These results thus quantify some opportunities and limitations of MC methods for spectral estimation of Hamiltonians. We numerically compare the MC eigenvalue estimation scheme (for stoquastic Hamiltonians) and the quantum-phase-estimation-based eigenvalue estimation scheme by implementing them for an archetypal stoquastic Hamiltonian system: the transverse field Ising chain.

Quantum phase estimation is a cornerstone in quantum algorithm design, allowing for the inference of eigenvalues of exponentially-large sparse matrices. The maximum rate at which these eigenvalues may be learned, –known as the Heisenberg limit–, is constrained by bounds on the circuit complexity required to simulate an arbitrary Hamiltonian. Single-control qubit variants of quantum phase estimation that do not require coherence between experiments have garnered interest in recent years due to lower circuit depth and minimal qubit overhead. In this work we show that these methods can achieve the Heisenberg limit, also when one is unable to prepare eigenstates of the system. Given a quantum subroutine which provides samples of a ‘phase function’ g(k) = ^P_j A_je^iφjk with unknown eigenphases φ_j and overlaps A_j at quantum cost O(k), we show how to estimate the phases {φ_j} with (root-mean-square) error δ for total quantum cost T = O(δ⁻¹). Our scheme combines the idea of Heisenberg-limited multi-order quantum phase estimation for a single eigenvalue phase [1, 2] with subroutines with so-called dense quantum phase estimation which uses classical processing via time-series analysis for the QEEP problem [3] or the matrix pencil method. For our algorithm which adaptively fixes the choice for k in g(k) we prove Heisenberg-limited scaling when we use the time-series/QEEP subroutine. We present numerical evidence that using the matrix pencil technique the algorithm can achieve Heisenberg-limited scaling as well.

Solid-state spin qubits is a promising platform for quantum computation and quantum networks^1,2. Recent experiments have demonstrated high-quality control over multi-qubit systems^3–8, elementary quantum algorithms^8–11 and non-fault-tolerant error correction^12–14. Large-scale systems will require using error-corrected logical qubits that are operated fault tolerantly, so that reliable computation becomes possible despite noisy operations^15–18. Overcoming imperfections in this way remains an important outstanding challenge for quantum science^15,19–27. Here, we demonstrate fault-tolerant operations on a logical qubit using spin qubits in diamond. Our approach is based on the five-qubit code with a recently discovered flag protocol that enables fault tolerance using a total of seven qubits^28–30. We encode the logical qubit using a new protocol based on repeated multi-qubit measurements and show that it outperforms non-fault-tolerant encoding schemes. We then fault-tolerantly manipulate the logical qubit through a complete set of single-qubit Clifford gates. Finally, we demonstrate flagged stabilizer measurements with real-time processing of the outcomes. Such measurements are a primitive for fault-tolerant quantum error correction. Although future improvements in fidelity and the number of qubits will be required to suppress logical error rates below the physical error rates, our realization of fault-tolerant protocols on the logical-qubit level is a key step towards quantum information processing based on solid-state spins.

The fault-tolerant operation of logical qubits is an important requirement for realizing a universal quantum computer. Spin qubits based on quantum dots have great potential to be scaled to large numbers because of their compatibility with standard semiconductor manufacturing. Here, we show that a quantum error correction code can be implemented using a four-qubit array in germanium. We demonstrate a resonant SWAP gate and by combining controlled-Z and controlled-S⁻¹ gates we construct a Toffoli-like three-qubit gate. We execute a two-qubit phase flip code and find that we can preserve the state of the data qubit by applying a refocusing pulse to the ancilla qubit. In addition, we implement a phase flip code on three qubits, making use of a Toffoli-like gate for the final correction step. Both the quality and quantity of the qubits will require significant improvement to achieve fault-tolerance. However, the capability to implement quantum error correction codes enables co-design development of quantum hardware and software, where codes tailored to the properties of spin qubits and advances in fabrication and operation can now come together to advance semiconductor quantum technology.

We propose and analyze two types of microwave-activated gates between a fluxonium and a transmon qubit, namely a cross-resonance (CR) and a CPHASE gate. The large frequency difference between a transmon and a fluxonium makes the realization of a two-qubit gate challenging. For a medium-frequency fluxonium qubit, the transmon-fluxonium system allows for a cross-resonance effect mediated by the higher levels of the fluxonium over a wide range of transmon frequencies. This allows one to realize the cross-resonance gate by driving the fluxonium at the transmon frequency, mitigating typical problems of the cross-resonance gate in transmon-transmon chips related to frequency targeting and residual ZZ coupling. However, when the fundamental frequency of the fluxonium enters the low-frequency regime below 100MHz, the cross-resonance effect decreases leading to long gate times. For this range of parameters, a fast microwave CPHASE gate can be implemented using the higher levels of the fluxonium. In both cases, we perform numerical simulations of the gate showing that a gate fidelity above 99% can be obtained with gate times between 100 and 300ns. Next to a detailed gate analysis, we perform a study of chip yield for a surface code lattice of fluxonia and transmons interacting via the proposed cross-resonance gate. We find a much better yield as compared to a transmon-only architecture with the cross-resonance gate as native two-qubit gate.

Leakage outside of the qubit computational subspace poses a threatening challenge to quantum error correction (QEC). We propose a scheme using two leakage-reduction units (LRUs) that mitigate these issues for a transmon-based surface code, without requiring an overhead in terms of hardware or QEC-cycle time as in previous proposals. For data qubits, we consider a microwave drive to transfer leakage to the readout resonator, where it quickly decays, ensuring that this negligibly disturbs the computational states for realistic system parameters. For ancilla qubits, we apply a |1↔|2π pulse conditioned on the measurement outcome. Using density-matrix simulations of the distance-3 surface code, we show that the average leakage lifetime is reduced to almost one QEC cycle, even when the LRUs are implemented with limited fidelity. Furthermore, we show that this leads to a significant reduction of the logical error rate. This LRU scheme opens the prospect for near-term scalable QEC demonstrations.

We analyze whether circuit QED Hamiltonians are stoquastic, focusing on systems of coupled flux qubits. We show that scalable sign-problem-free path integral Monte Carlo simulations can typically be performed for such systems. Despite this, we corroborate the recent finding [I. Ozfidan, Phys. Rev. Appl. 13, 034037 (2020)10.1103/PhysRevApplied.13.034037] that an effective, nonstoquastic qubit Hamiltonian can emerge in a system of capacitively coupled flux qubits. We find that if the capacitive coupling is sufficiently small, this nonstoquasticity of the effective qubit Hamiltonian can be avoided if we perform a canonical transformation prior to projecting onto an effective qubit Hamiltonian. Our results shed light on the power of circuit QED Hamiltonians for the use of quantum adiabatic computation and the subtlety of finding a representation which cures the sign problem in these systems.

Future fault-tolerant quantum computers will require storing and processing quantum data in logical qubits. Here we realize a suite of logical operations on a distance-2 surface code qubit built from seven physical qubits and stabilized using repeated error-detection cycles. Logical operations include initialization into arbitrary states, measurement in the cardinal bases of the Bloch sphere and a universal set of single-qubit gates. For each type of operation, we observe higher performance for fault-tolerant variants over non-fault-tolerant variants, and quantify the difference. In particular, we demonstrate process tomography of logical gates, using the notion of a logical Pauli transfer matrix. This integration of high-fidelity logical operations with a scalable scheme for repeated stabilization is a milestone on the road to quantum error correction with higher-distance superconducting surface codes.

We present a Dicke state preparation scheme which uses global control of N spin qubits: our scheme is based on the standard phase estimation algorithm, which estimates the eigenvalue of a unitary operator. The scheme prepares a Dicke state nondeterministically by collectively coupling the spins to an ancilla qubit via a ZZ interaction, using log2N+1 ancilla qubit measurements. The preparation of such Dicke states can be useful if the spins in the ensemble are used for magnetic sensing: we discuss a possible realization using an ensemble of electronic spins located at diamond nitrogen-vacancy centers coupled to a single superconducting flux qubit. We also analyze the effect of noise and limitations in our scheme.

We examine the problem of determining whether a multiqubit two-local Hamiltonian can be made stoquastic by single-qubit unitary transformations. We prove that when such a Hamiltonian contains one-local terms, then this task can be NP-hard. This is shown by constructing a class of Hamiltonians for which performing this task is equivalent to deciding 3-SAT. In contrast, we show that when such a Hamiltonian contains no one-local terms then this task is easy; namely, we present an algorithm which decides, in a number of arithmetic operations over R which is polynomial in the number of qubits, whether the sign problem of the Hamiltonian can be cured by single-qubit rotations.

We review some of the recent efforts in devising and engineering bosonic qubits for superconducting devices, with emphasis on the Gottesman-Kitaev-Preskill (GKP) qubit. We present some new results on decoding repeated GKP error correction using finitely-squeezed GKP ancilla qubits, exhibiting differences with previously studied stochastic error models. We discuss circuit-QED ways to realize CZ gates between GKP qubits and we discuss different scenarios for using GKP and regular qubits as building blocks in a scalable superconducting surface code architecture.

Leakage outside of the qubit computational subspace, present in many leading experimental platforms, constitutes a threatening error for quantum error correction (QEC) for qubits. We develop a leakage-detection scheme via Hidden Markov models (HMMs) for transmon-based implementations of the surface code. By performing realistic density-matrix simulations of the distance-3 surface code (Surface-17), we observe that leakage is sharply projected and leads to an increase in the surface-code defect probability of neighboring stabilizers. Together with the analog readout of the ancilla qubits, this increase enables the accurate detection of the time and location of leakage. We restore the logical error rate below the memory break-even point by post-selecting out leakage, discarding less than half of the data for the given noise parameters. Leakage detection via HMMs opens the prospect for near-term QEC demonstrations, targeted leakage reduction and leakage-aware decoding and is applicable to other experimental platforms.

One of the most direct preparations of a Gottesman-Kitaev-Preskill (GKP) qubit in an oscillator uses a tunable photon-pressure (also called optomechanical) coupling of the form qˆbˆ†bˆ, enabling us to imprint the modular value of the position qˆ of one oscillator onto the state of an ancilla oscillator. We analyze the practical feasibility of executing such modular quadrature measurements in a parametric circuit-QED realization of this coupling. We provide estimates for the expected GKP squeezing induced by the protocol and discuss the effect of photon loss and other errors on the resulting squeezing.