Benchmarks of Nonclassicality for Qubit Arrays

TL;DR

This paper introduces practical, resource-efficient benchmarks for testing multi-qubit nonclassicality and entanglement in qubit arrays, using fewer measurement settings than full tomography, and verifies their effectiveness through numerical simulations.

Contribution

It develops measurement-efficient benchmarks based on geometric structures from the $N$-qubit Pauli group, enabling verification of nonclassicality with minimal measurement settings.

Findings

Benchmarks require at most N+1 measurement settings.

Simulations confirm the bounds are tight and detect nonclassicality.

Results demonstrate the benchmarks' applicability to realistic superconducting qubit arrays.

Abstract

We present a set of practical benchmarks for $N$-qubit arrays that economically test the fidelity of achieving multi-qubit nonclassicality. The benchmarks are measurable correlators similar to 2-qubit Bell correlators, and are derived from a particular set of geometric structures from the $N$-qubit Pauli group. These structures prove the Greenberger-Horne-Zeilinger (GHZ) theorem, while the derived correlators witness genuine $N$-partite entanglement and establish a tight lower bound on the fidelity of particular stabilizer state preparations. The correlators need only $M \leq N+1$ distinct measurement settings, as opposed to the $2^{2N}-1$ settings that would normally be required to tomographically verify their associated stabilizer states. We optimize the measurements of these correlators for a physical array of qubits that can be nearest-neighbor-coupled using a circuit of controlled-$Z$ gates with constant gate depth to form $N$-qubit linear cluster states. We numerically simulate the provided circuits for a realistic scenario with $N=3,...,9$ qubits, using ranges of $T_1$ energy relaxation times, $T_2$ dephasing times, and controlled-$Z$ gate-fidelities consistent with Google's 9-qubit superconducting chip. The simulations verify the tightness of the fidelity bounds and witness nonclassicality for all nine qubits, while also showing ample room for improvement in chip performance.

Figures

Equations

$$\alpha =\sum _{i=1}^M{\lambda }_i{O}_i,$$

$$⟨\alpha ⟩=\sum _{i=1}^M{\lambda }_i\text{Tr}(\rho O_i)\le {\beta }_{\text{QM}}=M.$$

$$⟨\alpha ⟩\le {\beta }_{\text{LHVT}}=M-2.$$

$$⟨\alpha ⟩\le {\beta }_{\text{bisep}}=M-2,$$

$$\mathcal{B}=\frac{⟨\alpha {⟩}_{\text{exp}}-M+2}{2},$$

$$F\ge F_{\text{ID}}=\frac{⟨\alpha {⟩}_{\text{exp}}-M+4}{4}=\frac{\mathcal{B}+1}{2},$$

$$|\psi\rangle=\sum_{j=1}^{2^{N-M+1}}\Bigg{[}a_{j}|\kappa_{j1}\rangle+\sum_{k=2}^{M-1}b_{jk}|\kappa_{jk}\rangle\Bigg{]},$$

$$\langle\alpha_{1}\rangle_{\textrm{exp}}=\langle\psi|\alpha_{1}|\psi\rangle=\sum_{j=1}^{2^{N-M+1}}\Bigg{[}|a_{j}|^{2}\langle\kappa_{j1}|\alpha_{1}|\kappa_{j1}\rangle+\sum_{k=2}^{2^{M-1}}|b_{jk}|^{2}\langle\kappa_{jk}|\alpha_{1}|\kappa_{jk}\rangle\Bigg{]}.$$

$$\langle\alpha_{1}\rangle_{\rm exp}\leq\sum_{j=1}^{2^{N-M+1}}\Bigg{[}|a_{j}|^{2}M+\sum_{k=2}^{2^{M-1}}|b_{jk}|^{2}(M-4)\Bigg{]}=4F+M-4,$$

$$F\equiv ⟨\psi |{\Pi }_1|\psi ⟩\ge \frac{⟨{\alpha }_1{⟩}_{\exp }-M+4}{4}\equiv F_{\text{ID}}.$$

$$\rho =\sum _{j=1}^m{c}_l|{\psi }_l⟩⟨{\psi }_l|,$$

$$⟨{\alpha }_1{⟩}_{\text{exp}}\le \sum _{l=1}^m{c}_l(4F_l+M-4),$$

$$F\equiv \text{Tr}(\rho {\Pi }_1)=\sum _{j=1}^m{c}_l{F}_l\ge \frac{⟨{\alpha }_1{⟩}_{\text{exp}}-M+4}{4}\equiv F_{\text{ID}}.$$

$$\rho =\frac{1}{d}\sum _{i=1}^b{\lambda }_i{S}_i,$$

$$\rho_{AB}=\Big{(}\frac{1}{d^{A}}\sum_{i=1}^{b^{A}}\lambda^{A}_{i}S^{A}_{i}\Big{)}\otimes\Big{(}\frac{1}{d^{B}}\sum_{j=1}^{b^{B}}\lambda_{j}^{B}S^{B}_{j}\Big{)}=\frac{1}{d^{AB}}\sum_{k=1}^{b^{AB}}\lambda^{AB}_{k}S^{AB}_{k}.$$

$$D=\left(\begin{array}[]{cc}1&e^{-\Delta t/T_{2}}\\ e^{-\Delta t/T_{2}}&1\\ \end{array}\right)^{\otimes N},$$

$$\rho\rightarrow\int e^{-i\zeta_{i}(\pi/2+\delta\varphi)/2}\,\rho\,e^{i\zeta_{i}(\pi/2+\delta\varphi)/2}\,dP(\delta\varphi)=\frac{1}{2}\bigg{[}\rho+\zeta_{i}\rho\zeta_{i}+i(\zeta_{i}\rho-\rho\zeta_{i})\bigg{(}\frac{\sin w}{w}-\frac{\sin w}{2(w+\pi)}-\frac{\sin w}{2(w-\pi)}\bigg{)}\bigg{]},$$

Full text

Benchmarks of Nonclassicality for Qubit Arrays

Mordecai Waegell

Institute for Quantum Studies, Chapman University, Orange, CA, USA

  

Justin Dressel

Institute for Quantum Studies, Chapman University, Orange, CA, USA

Schmid College of Science and Technology, Chapman University, Orange, CA, USA

dressel@chapman.edu

Keck Center for Science and Engineering, Chapman University, One University Drive, Orange, CA 92866; Phone: (714) 516-5949

Keywords: Benchmark, Quantum Computing, Entanglement Witness, Bell Inequality, Quantum Circuits

Abstract: We present a set of practical benchmarks for $N$-qubit arrays that economically test the fidelity of achieving multi-qubit nonclassicality. The benchmarks are measurable correlators similar to 2-qubit Bell correlators, and are derived from a particular set of geometric structures from the $N$-qubit Pauli group. These structures prove the Greenberger-Horne-Zeilinger (GHZ) theorem, while the derived correlators witness genuine $N$-partite entanglement and establish a tight lower bound on the fidelity of particular stabilizer state preparations. The correlators need only $M\leq N+1$ distinct measurement settings, as opposed to the $2^{2N}-1$ settings that would normally be required to tomographically verify their associated stabilizer states. We optimize the measurements of these correlators for a physical array of qubits that can be nearest-neighbor-coupled using a circuit of controlled-$Z$ gates with constant gate depth to form $N$-qubit linear cluster states. We numerically simulate the provided circuits for a realistic scenario with $N=3,...,9$ qubits, using ranges of $T_{1}$ energy relaxation times, $T_{2}$ dephasing times, and controlled-$Z$ gate-fidelities consistent with Google’s 9-qubit superconducting chip. The simulations verify the tightness of the fidelity bounds and witness nonclassicality for all nine qubits, while also showing ample room for improvement in chip performance.

I Introduction

As hardware is developed to implement quantum circuits on increasing numbers of qubits, it will be valuable to have economical benchmarks of fully quantum behavior. From the outset of quantum computing it has been clear that the advantage of a quantum computer lies somewhere in its ability to readily perform tasks that are physically challenging or impossible for a classical system. Therefore, ideal hardware benchmarks should certify the ability of the hardware to generate such nonclassical behavior. Indeed, a wide variety of benchmarking techniques have been developed recently Friis2018 ; Gheorghiu2019 , including gate-fidelity benchmarks using randomized gate sequences that avoid the state-preparation and measurement errors, and state-preparation benchmarks that certify particular states while avoiding the exponential scaling of state tomography.

Despite these recent achievements, quantifying the specific nonclassical resources that lead to quantum computational advantage has remained an elusive goal vedral2010elusive . Several earlier proposals for suitable measures like entanglement hill1997entanglement ; wootters1998entanglement ; wootters2001entanglement ; wong2001potential , Bell-nonlocality EPR ; Bell1 ; Bell2 ; greenberger1989going ; greenberger1990bell ; Bravyi2018 , or quantum discord and its variations lanyon2008experimental ; ferraro2010almost ; Modi2012 , proved to be insufficient on their own due to the discovery of algorithmic counter-examples Bennett1999 ; Bennett1999b ; Meyer2000 ; dakic2010necessary ; Bera2017 . Recent advances suggest a strong connection between quantum advantage and contextuality KS ; galvao2005discrete ; howard2014contextuality ; wigner1970hidden ; klyachko2008simple , which is a general structural feature of quantum mechanics that subsumes nonlocality. The most pragmatic metric of nonclassical behavior in quantum devices, however, has been the violation of two-qubit Bell inequalities, or similar entanglement witnesses that can apply to few-qubit subsets of a multi-qubit device Friis2018b .

In this article we provide a set of practical hardware benchmarks that naturally generalize two-qubit Bell inequality tests to $N$ qubits, based on the Greenberger-Horne-Zeilinger (GHZ) theorem. As with Bell inequalities, our nonclassicality benchmarks use the experimental violation of a classical bound to quantify the nonclassical behavior of the circuit. Beyond quantifying nonclassicality via a bound-violation, these benchmarks also provide tight lower bounds on the fidelities with which particular stabilizer subspaces have been prepared, and thus witness genuine $N$-qubit entanglement for all states that lie within the targeted subspaces. These benchmarks are optimized for testing controllable qubit arrays with nearest-neighbor coupling. As such, we provide efficient circuits for preparing cluster states that maximally violate these benchmarks with controlled-$Z$ entangling gates, using a constant gate depth of 4 (up to hardware-specific decompositions of the controlled-$Z$ gate Chow2011 ; Ghosh2013 ; Martinis2014 ; Barends2014 ; kelly2015state ; Chow2014 ). Though our benchmarks efficiently verify genuine $N$-qubit entanglement using cluster states, many of the benchmarks may be applied to other stabilizer states and we expect similar benchmarks to exist for all stabilizer states.

The benchmarks we present here generalize earlier work that was experimentally tested with $N=3,4$ photons greganti2015practical , where they were compared to previously proposed state-dependent methods for efficiently verifying the fidelity of particular entangled $N$-qubit preparations Guehne2007 ; Wunderlich2011 . These prior methods have already been used to verify multi-qubit entanglement in state-of-the-art experiments with 12 qubits Gong2019 and 18 qubits Wang2018 , since the exponential scaling required for traditional state tomography is increasingly prohibitive. Notably, for large $N$ our GHZ-based benchmarks produce a tighter preparation-fidelity bound than these existing methods and similarly produce entanglement witnesses with better scaling.

II Results

Nonclassicality Benchmarks:— Our benchmarks consist of measurable correlators that are compared to derived upper bounds; violation of these bounds characterizes nonclassicality. Each such benchmark corresponds to a specific prepare-and-measure circuit on $N$-qubits with $M\leq N+1$ different measurement settings. The $M$ observables form a structure called an ID (also called an identity product W_Primitive ), which is a set of mutually commuting $N$-qubit Pauli operators whose overall product is the $N$-qubit identity, up to a sign. We express an ID as an $M\times N$ table of single-qubit Pauli operators and the identity $\{Z,X,Y,I\}$, labeled $O_{ij}$ with $i=1,...,M$ and $j=1,...,N$. We also define the shortened label $O_{i}=\bigotimes_{j=1}^{N}O_{ij}$ to indicate the $N$-qubit observable obtained as the product of the $i$th row of an ID. We omit tensor product symbols for compactness.

To obtain the Bell inequality for each ID greganti2015practical , we choose a particular eigenspace $\Pi$ represented by a projector of rank $2^{N-M+1}$, which is specified by the set of $N$-qubit Pauli observables $\{O_{i}\}$ that form the $M$ rows of the ID (see Figs. 1 and 2), and a specific choice of their respective eigenvalues $\{\lambda_{i}\}$. We then define the correlator observable for this chosen eigenspace,

[TABLE]

such that its expectation value in a state $\rho$ has an upper bound of $\beta_{\textrm{QM}}=M$, saturated by the chosen eigenspace $\rho=\Pi$,

[TABLE]

For example, we could prepare the joint eigenstate of the ID of Fig. 1(a), with negative eigenvalue $\lambda_{1}=-1$ for the 3-qubit Pauli observable $O_{1}=YXY$, and positive eigenvalues $\lambda_{2}=\lambda_{3}=\lambda_{4}=+1$ for the remaining observables $O_{2}=YYZ$, $O_{3}=ZXZ$, and $O_{4}=ZYY$. Then, $\langle\alpha\rangle=\textrm{Tr}(\Pi\alpha)=4$ since each term in the sum becomes +1.

In the spirit of Bell Bell1 ; Bell2 , if one tries to explain the observed correlation by choosing a complete set of local hidden variables $v_{Zj},v_{Xj},v_{Yj}\in\{+1,-1\}$ that predict the outcomes of the single-qubit Pauli measurements, then at least one of the terms in the correlator sum becomes -1, resulting in a smaller upper bound,

[TABLE]

Experimental violation of this bound thus indicates nonclassicality in the form of a violation of local realism. Though the locality loophole is always open for neighboring qubits on a chip, this violation is still a useful witness for nonclassical states prepared by the chip, much like for Bell inequalities or Bell-Leggett-Garg inequalities white2016preserving . The derivation of this bound is reviewed in the Methods Section.

As an independent result, maximizing the expectation value of the correlator over all biseparable quantum states in the $N$-qubit Hilbert space produces the upper bound,

[TABLE]

which happens to coincide with the bound for local hidden variable theories. Experimental violation of the bound thus also witnesses genuine $N$-partite entanglement. In the Methods section, we provide the proof that the joint eigenspaces of the IDs in this article are maximally entangled, as well as the derivation of this bound.

In light of the convenient fact that $\beta_{\textrm{bisep}}=\beta_{\textrm{LHVT}}$, we define the nonclassicality benchmark score for a given physical $N$-qubit device as the experimentally determined value,

[TABLE]

such that $\mathcal{B}\leq 0$ fails to witness either entanglement or the violation of local realism, while $0<\mathcal{B}\leq 1$ witnesses nonlocal $N$-partite-entangled states. The nonclassicality benchmark score thus serves as a metric of uniquely quantum behavior, with $\mathcal{B}=1$ indicating maximum nonclassicality that saturates the correlator bound. Each $N$-qubit ID provides a benchmark corresponding to a distinct nonclassical eigenspace of an $N$-qubit physical device, and thus the hierarchy of IDs presented in Fig. 1 provides a corresponding hierarchy of benchmarks.

Lower Bounding the Fidelity:— The correlator also serves to bound the fidelity from below greganti2015practical ,

[TABLE]

where $F=\textrm{Tr}(\rho_{\textrm{exp}}\Pi)\in[0,1]$ is the fidelity that the experimentally prepared state $\rho_{\textrm{exp}}$ lies within the eigenspace $\Pi$ stabilized by the chosen ID. We provide a general derivation of this bound in the Methods section. Importantly, in the limit $\langle\alpha\rangle_{\textrm{exp}}\to\beta_{\textrm{QM}}=M$, we have $F_{\textrm{ID}}\to 1$, and thus as the fidelity of the preparation is improved, this lower bound obviates the need for full tomography of these preparations.

Taken together, the inequalities of Eqs. 3, 4, and 6 provide a practical and efficient characterization of the prepared $N$-qubit state, as well as a robust benchmark of its nonclassical behavior, using only $M\leq N+1$ measurement settings. We present minimal benchmark IDs in Fig. 1 for $N=3,\ldots,9$, and detail minimal IDs up to $N=33$ qubits in Supplementary Figures 1 through 5. These minimal IDs saturate the conjectured bound $N\leq(M-2)(M-1)/2$. We also present a family of maximal benchmark IDs in Fig. 2 for all $N\geq 10$ that saturate the bound $M-1\leq N$.

Benchmark Circuits and Simulation:— The IDs in this article have been specially chosen so that the prepare-and-measure circuit for each measurement setting requires a gate depth of 4 on any array of $N$ physical qubits with only nearest-neighbor controlled-$Z$ couplings, making them a scalable and uniform set of benchmarks for implementations of this type. Figure 3 shows the circuits for $N=4,5$, from which the generalization to all $N$ should be straightforward. In general, each circuit prepares an $N$-qubit linear cluster state, which is contained within the maximally entangled subspace of the corresponding ID.

In order to evaluate the usefulness of these benchmarks in real-world physical implementations, we simulated the performance of these circuits for each of the IDs in Fig. 1. We simulated each circuit over a range of $T_{1}$ energy relaxation times, $T_{2}$ dephasing times, and angular jitter for the controlled-$Z$ gate rotations, using the ranges given in Figs. 4 and 5. We also considered the effect of initialization and readout error for each qubit. The ranges of values were chosen to match the reported values of the 9-qubit Google chip Barends2014 ; kelly2015state , with the experimental values roughly in the center of each simulated range. We ran one version of the simulation using a nominal initialization error for each qubit of $P_{\textrm{e}}=2\%$, and another version where we used the observed initialization errors for each of the nine qubits on the Google chip. Final readout error has been neglected as correctable for ensemble statistics. Selected plots from the simulations are shown in Fig. 4, while scatter plots of the lower fidelity bound, $F_{\textrm{ID}}$, are shown in Fig. 5 for the full ranges of simulated values. Note that in order to minimize the effect of the two worst qubits on the chip (boldface values in Figs. 4 and 5), we always used the last $N$ qubits on the chip to form our $N$-qubit IDs in the simulation. See the Methods section for additional details about how the numerical simulations were performed.

Judging by our simulated data shown in Figs. 4 and 5, we expect the 9-qubit Google chip to be able to violate the classicality bounds for all nine qubits. We can see clearly that the qubit initialization error is the dominant source of error as we try to move to larger $N$. This shows that our benchmarking scheme is immediately relevant, since it appears that similar hardware fidelity would only violate the bound for one or two more qubits — but certainly not all 72 on the Bristlecone chip Note1 — once suitable IDs have been found beyond the 9 presented here.

III Discussion

The IDs and implementation circuits presented in this article are good benchmark tests for any physical implementation of qubits in a nearest-neighbor-connected array. They work naturally on a chip with more connectivity than this as well. While our simulations targeted a particular recent chip implementation for concreteness, this does not constrain the general usefulness of this protocol for other multi-qubit systems.

Although some other families of IDs with the same properties as those in Figs. 1 and 2 are known W_Primitive ; waegell2013nonclassical , the minimal IDs, with the largest possible value of $N$ for a given $M$, are not known in general (see the Supplementary Discussion and Supplementary Figures 1 through 5 for the best known cases). Because of their geometric nature, enumerating all of the representative IDs for given values of $N$ and $M$ is a highly nontrivial problem, related to solving the graph isomorphism problem on $N\times M$ colored vertices, and it is thus limited by computational resources. Furthermore, not every ID can be constructed using only nearest-neighbor couplings in linear circuits as in Fig. 3. The increased connectivity of more modern chips, like the Bristlecone chip from Google, should allow the use of more general IDs, although the circuit depth will likely increase by one or two gates.

Each of the IDs presented here also gives rise to a complete proof of the Kochen-Specker (KS) theorem for contextuality KS ; waegell2012proofs ; waegell2013proofs , which can be implemented for any initial state with a few alternative circuits for the different measurement contexts. In general, IDs are the natural building blocks of proofs of the KS theorem in the $N$-qubit Pauli group. This is a slightly more complicated setup, which could inspire different contextuality-based benchmarks in future work.

Finally, maximally entangled IDs with $M<N+1$ give rise to maximally entangled eigenspaces, each of dimension $2^{N-M+1}$, which generalize the codespaces of error correcting codes knill1997theory ; nielsen2010quantum , and $L=N-M+1$ is the number of logical qubits (where $N$ is the number of physical qubits). All $N$-qubit-stabilizer-based error correcting codes (including the toric code kitaev2003fault ) belong to the family of IDs, and while all IDs of this type are error-detecting codes, they cannot all be used to diagnose the syndrome of an error in order to correct it. Many of the well-known error correcting codes generate an ID which proves the GHZ theorem, and all can be used as entanglement witnesses in the manner of this article divincenzo1997quantum . Nevertheless, these more general maximally entangled subspaces may be of significant interest for other applications in quantum information processing, which warrants further investigation. One straightforward application for these subspaces is to perform benchmarks that measure the physical qubits as described in this paper, while simultaneously benchmarking the performance of the logical qubits in some additional way. The two tests may be performed simultaneously because any general logical $L$-qubit state can be prepared for each benchmark, although the circuit is likely to be longer and more complex than Fig. 3, and the performance will be commensurately worse. It is remarkable to note that if the conjectured bound $N\leq(M-2)(M-1)/2$ can be saturated, then the number of logical qubits is bounded by $L\leq((M-2)(M-1)/2-M+1$, and thus the ratio $L/N\rightarrow 1$ in the limit $M\rightarrow\infty$.

IV Methods

Proving the GHZ Theorem: All of the IDs in Fig. 1 have sign -1, and for each qubit $j$, the number of entries $O_{ij}=Z$ in the ID is even, as is the number of entries with $O_{ij}=X$ and with $O_{ij}=Y$. These properties indicate that these IDs give rise to proofs of the GHZ theorem greenberger1989going , which is a logical version of Bell’s nonlocality theorem Bell1 ; Bell2 , without any inequalities. To see this, suppose that a joint eigenstate (i.e., any state in a joint eigenspace) of these observables is prepared. This eigenstate has $M$ eigenvalues $\lambda_{i}$ corresponding to the $M$ observables, and $\prod_{i=1}^{M}\lambda_{i}=-1$, since the product of these $M$ observables is $-I^{\otimes N}$. Suppose that each of the $N$ qubits are now mutually space-like separated, and each is subjected to random local Pauli measurements, and label their outcomes $\lambda_{ij}$, when all $N$ local measurement settings happen to correspond to observable $i$ of the ID. The entanglement correlations that are obeyed by this state are $\prod_{j=1}^{N}\lambda_{ij}=\lambda_{i}$. Putting these relations together we have $\prod_{i=1}^{M}\prod_{j=1}^{N}\lambda_{ij}=-1$. Now, in order for a local hidden variable theory (LHVT) to explain these entanglement correlations, each qubit $j$ must carry local hidden variables $v_{Zj},v_{Xj},v_{Yj}\in\{+1,-1\}$ which predict the outcomes $\lambda_{ij}$, and are pre-arranged to satisfy the entanglement constraints. However, for such hidden variables we would have $\prod_{i=1}^{M}\prod_{j=1}^{N}\lambda_{ij}=\prod_{j=1}^{N}v_{Zj}^{n_{j}}v_{Xj}^{m_{j}}v_{Yj}^{l_{j}}=+1$, since $n_{j}$, $m_{j}$, and $l_{j}$ are all even for the IDs of this article, and thus is is impossible to choose local hidden variables which can satisfy the entanglement correlations of this state. This logical proof without inequalities can be converted into a Bell inequality for use as a benchmark of $N$-qubit nonlocality, as shown in the main text, by noting that for any complete assignment of local hidden variables $v_{Zj},v_{Xj},v_{Yj}\in\{+1,-1\}$ to the ID, at least one of the observables has the wrong eigenvalue.

In general, proving the GHZ theorem does not prove that nonlocal correlations exist between more than just a single pair of qubits among the $N$ collins2002bell ; svetlichny1987distinguishing ; seevinck2002bell ; mitchell2004conditions , nor does it generally witness genuine $N$-qubit entanglement. In contrast, the benchmark IDs we present in this article prove the GHZ theorem and are constructed to be $N$-partite entanglement witnesses toth2005detecting ; toth2005entanglement , such that their corresponding Bell inequalities can only be violated by genuinely $N$-qubit-entangled states. To go further than the results we present here and prove nonlocal correlations exist between every pair of qubits among the $N$, one must violate the corresponding Svetlichny inequalities collins2002bell ; lavoie2009experimental instead, but with the cost that the number of required measurement settings grows exponentially with $N$ collins2002bell .

Bounding the Fidelity: An $N$-qubit ID with $M$ observables $\{O_{i}\}$ has a complete set of eigenspaces $\{\Pi_{k}\}$ satisfying $\sum_{k}\Pi_{k}=I$, each of which can be identified by a unique set of distinct eigenvalues $\{\lambda_{ik}\}$ of $\{O_{i}\}$. Only $M-1$ of the observables in an ID are independent, and if $M-1<N$ the eigenspaces $\Pi_{k}$ are degenerate, and each contains $2^{N-M+1}$ mutually orthogonal vectors $|\kappa_{jk}\rangle$ which share the eigenvalue $\lambda_{ik}$, with $j=1,\ldots,2^{N-M+1}$, such that $\{|\kappa_{jk}\rangle\}$ is a complete orthonormal eigenbasis of the ID. Each of the $2^{M-1}$ eigenspaces $\Pi_{k}$ corresponds to a unique correlator $\alpha_{k}=\sum_{i=1}^{M}\lambda_{ik}O_{i}$. Each experimentally obtained quantity $\langle\alpha_{k}\rangle$ enables us to put a lower bound on the fidelity that an experimentally prepared pure state $|\psi\rangle$ lies within the eigenspace $\Pi_{k}$ greganti2015practical .

With no loss of generality, we will henceforth use correlator $\alpha_{1}$ and the target eigenspace $\Pi_{1}$. We begin by expanding $|\psi\rangle$ in this eigenbasis as,

[TABLE]

such that \sum_{j}\big{(}|a_{j}|^{2}+\sum_{k=2}^{2^{M-1}}|b_{jk}|^{2}\big{)}=1.

Since the expansion is in an eigenbasis of $\alpha_{1}$, we find

[TABLE]

Note that $\langle\kappa_{j1}|\alpha_{1}|\kappa_{j1}\rangle=\sum_{k=1}^{M}\lambda_{k}^{2}=M$, since all eigenvalues of $|\kappa_{1}\rangle$ match those in the correlator $\alpha_{1}$ by construction, and thus square to 1. However, any other $|\kappa_{jk}\rangle$ does not lie within $\Pi_{1}$, so is characterized by eigenvalues distinct from those characterizing $\Pi_{1}$. Moreover, since the product of all eigenvalues for the observables of a given ID is fixed for any eigenstate, only even numbers of eigenvalues can differ from those characterizing $\Pi_{1}$, which necessarily causes at least two terms of $\langle\kappa_{jk}|\alpha_{1}|\kappa_{jk}\rangle$ to become $-1$, resulting in an upper bound of $\langle\kappa_{jk}|\alpha_{1}|\kappa_{jk}\rangle\leq M-4$ for those eigenstates. Using these two observations we obtain,

[TABLE]

where $F=\sum_{j}|a_{j}|^{2}$, and we have used \sum_{j}\big{(}|a_{j}|^{2}+\sum_{k=2}^{2^{M-1}}|b_{jk}|^{2}\big{)}=1. We can rewrite this relation as

[TABLE]

Noting that the left hand side of this equation is the fidelity $F$ for the preparation $|\psi\rangle$ to lie within the eigenspace $\Pi_{1}$, the right hand side $F_{\textrm{ID}}$ gives a lower bound $F\geq F_{\textrm{ID}}$ for the fidelity. For IDs with $M=N+1$, the target subspace $\Pi_{1}$ contains only one eigenvector, so the fidelity $F$ is also a state preparation fidelity for the particular target eigenstate $|\kappa_{1}\rangle$. For IDs with $M<N+1$, the target subspace $\Pi_{1}$ is degenerate, so the fidelity $F$ is the fidelity for $|\psi\rangle$ to lie within that subspace.

Next we generalize the above derivation to the case of mixed states. For a general convex combination of $m$ pure states,

[TABLE]

where $\sum c_{l}=1$, we can expand each $|\psi_{l}\rangle$ using appropriate eigenbases of the ID as in Eq. (7) and follow the same arguments to obtain

[TABLE]

where $F_{l}\equiv\langle\psi_{l}|\Pi_{1}|\psi_{l}\rangle$. We can rewrite this as,

[TABLE]

As in the pure state case, the left hand side is the fidelity $F$ for the mixed state $\rho$ to lie within the target subspace $\Pi_{1}$, while the same expression for the right hand side $F_{\textrm{ID}}$ places a lower bound on this fidelity.

Witnessing Genuine $N$-Partite Entanglement: An $N$-qubit ID provides an entanglement witness if it is maximally entangled W_Primitive ; waegell2014bonding . Entanglement is usually discussed in reference to the separability of states. However, there is a way to reason about the entanglement of a set of observables directly without reference to states. We define a maximally entangled set of $N$-qubit observables as one with the property that there exists no bipartition of the $N$ qubits into subsets of $R$ and $N-R$, such that all of the observables in each subset $\bigotimes_{k\in[1,R]}O_{ik}$ mutually commute. It follows from this definition that the joint eigenstates of this set are maximally entangled $N$-qubit stabilizer states.

To see this, consider that every stabilizer state (space) of $N$ qubits has a stabilizer group of $b=2^{g}$ mutually commuting Pauli observables $\{S_{i}\}$ and corresponding eigenvalues $\{\lambda_{i}\}$, and its density operator can be written as,

[TABLE]

where $g$ is the number of independent generators in the set, and $d=2^{N}$ is the dimension of the Hilbert space. Note that if $g<N$, then $\rho$ projects onto a subspace of rank $r=2^{N-g}>1$, and that $g=M-1$ for a minimal ID, which is just a specific subset of one or more complete stabilizer groups. If a stabilizer state is the tensor product of two smaller stabilizer states on subsystems $A$ and $B$, it follows that its density operator can be written as,

[TABLE]

For the bipartition of the system into $A$ and $B$, all of the stabilizer operators $S^{A}_{i}=\bigotimes_{k\in A}O_{ik}$ mutually commute by definition. It follows that one can find such a mutually commuting bipartition for any separable state, and therefore if no such bipartition exists, then the set of observables is maximally entangled. All of the IDs presented in this article are maximally entangled in this way, which results in a witness inequality with the same bound as the Bell inequality.

All states within a maximally entangled eigenspace of an ID are maximally entangled, meaning that for all of them, the maximum squared-Schmid-coefficient across all bipartitions is $1/2$. For such an eigenstate $|\psi\rangle$, a standard entanglement witness is $\mathcal{W}=\mathbb{1}/2-|\psi\rangle\langle\psi|$, and an experimental measurement of $\langle\mathcal{W}\rangle<0$ is a witness of genuine $N$-partite entanglement toth2005entanglement . Noting that a superposition state $a|\psi\rangle+b|\psi_{\bot}\rangle$ can only violate this bound for $F=|a|^{2}>1/2$, we obtain $F_{\textrm{ID}}\leq F\leq 1/2$ for all biseparable states. Plugging this into $F_{\textrm{ID}}=(\langle\alpha\rangle_{\textrm{exp}}-M+4)/4$ yields $\langle\alpha\rangle_{\textrm{bisep}}\leq M-2$, which is Eq. (4).

Numerical Simulation Details: In the simulation, the state is first degraded by initialization error. That is, ideally the $N$ qubits are prepared in an initial ground state $\otimes_{i=1}^{N}|0\rangle$. However, each qubit has an error probability $P_{\rm e}^{(i)}$ of being initially excited, which produces a mixed initial bit state $(1-P_{\rm e}^{(i)})|0\rangle\langle 0|+P_{\rm e}^{(i)}|1\rangle\langle 1|=(1-2P_{\rm e}^{(i)})|0\rangle\langle 0|+P_{\rm e}^{(i)}I$, and thus a degraded initial state $\rho=\otimes_{i=1}^{n}[(1-2P_{\rm e}^{(i)})|0\rangle\langle 0|+P_{\rm e}^{(i)}I]$ with ground state fidelity $\prod_{i=1}^{n}(1-2P_{\rm e}^{(i)})$. The final readout error for an ensemble average can be corrected if the readout misidentification probabilities $P_{\rm e}^{i}$ are known, and thus we have neglected the role of the readout error.

Each gate in Fig. 3 is then applied to the initial state $\rho$. For the Hadamard gate, it is sufficient to use a $Y_{90}$ rotation, $\exp(-iY\pi/4)$. We decompose the controlled-$Z$ gate into an implementable $ZZ_{90}$ entangling gate and single-qubit corrections: $\exp(i\pi/4)[\exp(iZ\pi/4)\otimes\exp(iZ\pi/4)]\exp(-iZZ\pi/4)$. We degraded each gate by $T_{1}$ energy relaxation and $T_{2}$ dephasing processes for the corresponding gate times $\Delta t$. For the energy relaxation time $T_{1}$, the first-order corrections for each individual qubit are accumulated and then applied to $\rho$. For each qubit \Delta\rho_{i}=\big{(}a^{\dagger}_{i}\rho a_{i}-\frac{1}{2}\{\rho,a^{\dagger}_{i}a_{i}\}\big{)}\Delta t/T_{1}^{i}, where $a_{i}$ is the lowering operator of the $i$th qubit tensored with identity for the other qubits, and $\rho\rightarrow\rho+\sum_{i}^{N}\Delta\rho_{i}$. This linear-order Lindblad-form update is sufficient since $\Delta t/T_{1}^{i}\ll 1$. For the dephasing time $T_{2}$, we directly construct the matrix,

[TABLE]

for efficiency and apply gate dephasing using element-wise multiplication (MATLAB syntax .), as $\rho\to\rho$ . $D$.

For simulating gate infidelity, we assume that the single-qubit gate fidelities are high enough for their errors to be neglected, and so simulate only a range of fidelities for the 2-qubit controlled-$Z$ gates. As a crude model for infidelity of a controlled-$Z$ gate, we add a random angular jitter $\delta\varphi$ only to the $ZZ$ rotation step, $\exp[-iZZ(\pi/2+\delta\varphi)/2]$, and average over the effect of this jitter using a raised cosine distribution with a width $w$, $dP(\delta\varphi)=d(\delta\varphi)[1+\cos(\pi\delta\varphi/w)]/(2w)$, where $\delta\varphi\in[-w,w]$ has compact angular support. This yields the averaged state update,

[TABLE]

where $\zeta_{i}$ is the tensor product of Pauli $Z$ for the two qubits the controlled-$Z$ is acting on, and identity for all of the other qubits. The limit as $w\to 0$ restores the unperturbed gate. This crude error model includes only one possible physical mechanism of infidelity for the controlled-$Z$ gate, but gives an indication of the gate sensitivity to imprecise angular control. Since the initialization error dominates the infidelity, the effect of the angular jitter is small.

Code availability: The MATLAB code used to generate our data is available from the authors on reasonable request.

Data availability: The data that support our findings are available from the authors on reasonable request.

Acknowledgments: We thank Josh Mutus and Daniel Sank for helpful commentary, as well as Eric Freda for helping to create some of the figures in this paper. MW was partially supported by the Fetzer Franklin Fund of the John E. Fetzer Memorial Trust. JD was partially supported by the Army Research Office (ARO) grant No. W911NF-15-1-0496, as well as No. W911NF-18-1-0178.

Competing Interests: There are no competing interests.

Author Contributions: MW developed the IDs for the benchmarks in this article, developed the benchmark inequalities, coded the simulations, and wrote the manuscript. JD developed the theory for the simulations, and co-wrote the manuscript.