Correcting detection error in quantum computation and state engineering through data processing

# Correcting detection error in quantum computation and state engineering through data processing

C. Shen and L.-M. Duan Department of Physics and MCTP, University of Michigan, Ann Arbor, Michigan 48109, USA Center for Quantum Information, IIIS, Tsinghua University, Beijing, China
###### Abstract

Quantum error correction in general is experimentally challenging as it requires significant expansion of the size of quantum circuits and accurate performance of quantum gates to fulfill the error threshold requirement. Here we propose a method much simpler for experimental implementation to correct arbitrary detection errors. The method is based on processing of data from repetitive experiments and can correct detection error of any magnitude, as long as the error magnitude is calibrated. The method is illustrated with its application to detection of multipartite entanglement from quantum state engineering.

###### pacs:
03.67.Mn, 42.50.Dv, 03.67.Bg, 03.65.Ud

A celebrated achievement in quantum information science is establishment of the error threshold theorem, which assures that any experimental error below a certain threshold value can be corrected in quantum computation or communication 1 (); 1' (). This theorem is important as it indicates there is no fundamental obstacle to realization of quantum computation with imperfect experimental devices, as long as the imperfection is small. The experimental realization of fault-tolerant quantum error correction, unfortunately, is still challenging. When all the experiment devices are subject to errors as it is the case for real applications, realization of fault-tolerance requires significant expansion of the size of quantum circuits for complicated encoding and accurate performance of quantum gates to fulfill the error threshold requirement, which is still beyond the capability of current experimental technology.

In this paper, we propose an alternative method to correct a special but practically important type of error — the detection error — in quantum information processing. The detection error is a common source of noise in many quantum information experiments. In particular, it is a significant obstacle to observation of multipartite entanglement in quantum state engineering 2 (). We show here that this type of error can be corrected at any magnitude as long as the error magnitude has been calibrated (for instance, through prior test experiments). The detection error distorts the experimental data by a transformation that depends on the magnitudes of various error possibilities. When the relevant error magnitudes have been calibrated for the detectors by the prior test experiments, the form of the distortion transformation induced by the detection error is known, and then we can find a way to inverse this transformation to reconstruct the original signal. In this way, we can use imperfect detectors to simulate perfect detectors as long as their imperfection has been calibrated. The proposed method is straightforward for experimental implementation as it is based on data processing and removes the difficulty associated with fault-tolerant quantum encoding. To correct the detection error, we only require to repeat the same experiments by some additional rounds to have small statistical error for the inverse transformation. To illustrate its applications, the method is used to significantly improve the detection of multi-qubit entanglement and spin squeezing. In many cases, the signal of multipartite entanglement only becomes visible after the proposed correction of the detection error, in particular when the number of qubits is large.

Any measurements in quantum information can be reduced to population measurements in certain bases (including possibly several complementary bases). If we want to measure properties associated with a state (generally mixed) of qubits, in each chosen measurement basis, there are possible measurement outcomes. By measurements we determine the probability associated with each outcome (). For instance, if we repeat the same experiment times and get the th outcome  times, we estimate the probability by and its standard deviation (the error bar) by using the binomial distribution. If the detectors are perfect, the measured probabilities just give the distribution of the state in the measurement basis . In reality, however, the detectors always have errors, which distort the distribution , making the measured distribution significantly different from . The purpose of this paper is to show how to reconstruct the real distribution from the measured distorted probabilities .

We first consider the case where the measurements haven individual addressing, and each qubit is measured by an independent detector. For detection on a qubit, the most general error model is characterized by a matrix

 D=[1−p0p1p01−p1], (1)

where () denotes respectively the error probability that the detector gives outcome () for the input signal of (). For simplicity of notation, we assume the error matrix has the same form for detection of each qubit (it is straightforward to generalize the formalism to the case where the error rates and in the matrix are qubit-dependent). Furthermore, we assume and have been well calibrated by a prior test experiment. For instance, we may input a known state to the detector and can calibrate and easily from the measurement data.

For qubits, the error model for the detection is then characterized by a matrix , with the element corresponding to the probability of recording the outcome with the input signal . Assume the detection error rates on different qubits are independent to each other and the binary string has zeros and ones. If we need flips from to and flips from to to change the string from to , the matrix element is given by

 Mji=(1−p0)n0−α(1−p0)n1−βpα0pβ1. (2)

The measured probabilities are connected with the real distribution through the distortion transformation To reconstruct the real signal from the measured distribution , in principle we only need to inverse the matrix . However, as is a huge matrix, it is not clear how to inverse this matrix (it is even a question whether the inverse exists).

Our key observation is that the matrix , with the elements given by Eq. (2), has a simple tensor product structure. It is straightforward to show by mathematical induction that

 M=n⨂k=1Dk, (3)

where all the are identical and given by in Eq. (1). Therefore, the inverse can be easily done in an analytic form with

 M−1=n⨂k=1D−1k=n⨂k=1[1−p′0p′1p′01−p′1]i, (4)

where the parameters and are given by

 p′0 = p0/(p0+p1−1), p′1 = p1/(p0+p1−1). (5)

Note that with the substitution in Eq. (5), and have the same form except that and can not be interpreted as error rates any more since in general they are not in the range . The formula also shows that the inverse transformation always exists except for the special case with .

In some experimental systems we do not have the ability to resolve individual qubits. Instead, we perform collective measurements on qubits by detecting how many qubits (denoted by , ) are in the state in a chosen detection basis (this is equivalent to measurement of the collective spin operator along a certain direction). In this case, the detection only has outcomes for an -qubit system. For collective measurements on qubits, the detection error matrix is represented by an matrix . The matrix element corresponds to the probability to register outcome when qubits are in the state. If the detection error matrix for an individual qubit is still given by in Eq. (1), we can directly calculate from : from signal to , if qubits flip from to and qubits flip from to , with the constraints and , is given by

 Lij = ∑n01,n10B(j,p1,n01)B(n−j,p0,n10) (6) = ∑qB(j,1−p1,q)B(n−j,p0,i−q),

where and we have let , and hence satisfies the constraint . As the dimension of the matrix depends linearly on the qubit number , it is typically not difficult to numerically calculate its inverse matrix if is not very large. As will be shown in the appendix, there is also a simple analytic formula for : if we denote the dependence of in Eq. (6) on as , we have

 L−1ij=Lij(p′0,p′1) (7)

where are given by the simple substitution in Eq. (5). With the inverse matrix , the real signal can be similarly reconstructed from the measured data as .

The above formulation can be extended straightforwardly to qudit (-dimensional) systems where the individual detection error matrix in Eq. (1) is replaced by a matrix. For independent detection of -qudits, the overall error matrix still has the tensor-product structure as shown by Eq. (3), which allows easy calculation of from .

With the inverse error matrix , it is straightforward to reconstruct the true distribution from the measured data . The price we need to pay is that compared with , there is an increase of the standard deviation (error bar) in our estimate of by the formula . With some tedious but straightforward calculation, we find

 Δgi=√[∑j(M−1ij)2fj−g2i]/N (8)

As and has matrix element (when , has matrix element which leads to exponential increase of the error bar with the qubit number . To maintain the same error bar , the number of repetitions of the experiment eventually needs to increase exponentially with . For practical applications, this exponential increase of by the factor is typically not a problem for two reasons. First, as the detection error rate is usually at a few percent level, the exponential factor remains moderate even for hundreds of qubits. Second, this exponential increase only applies when we need to measure each element of the distribution . In most of quantum information applications, we only need to measure certain operators which are expressed as tensor products of a constant number of Pauli operators for different qubits. In this case, does not have the exponential increase as we show now.

Suppose we need to measure an operator , which is expressed as , where is a component of the Pauli matrices when or the identity operator when . The number of the Pauli matrices in the tensor product expansion of is called the support of . To measure the operator , we choose the measurement basis to be the eigenbasis of for the th qubit. In this measurement basis, is diagonal with the matrix element diag, where diag for and diag for . Under the distribution , the expectation value of is given by , where denotes the diagonal matrix element of . Therefore, by defining a corrected operator , we can get the true expectation value directly from the experimental data . Using the relation , is expressed as For ,

 diag(σμkk)D−1k = [(1−2p′0)−(1−2p′1)] (9)

and for , diag. For simplicity of notation, we take . In this case, diagdiag for , and the corrected operator is related with the original operator by a simple scaling transformation .  The scaling transformation is independent of the qubit number , so the error bar of does not have exponential increase with when the operator has a constant support .

The scaling transformation also applies to collective operators, but some caution needs to be taken for calculation of their variance. For instance, if we take the collective spin operator , it is easy to see that  as each of the terms of has support . However, as which has non-uniform support for its superposition terms, one finds that . With this transformation, we can correct the distortion to the spin squeezing parameter by the detection error. Assume that the mean value of is along the -direction with and the squeezing is along the -direction. The squeezing parameter is given by 3 (). Using the transformation for and , we find that

 (10)

where is the contribution to by the detection noise. After correction of the detection error, gets significantly smaller compared with in particular when the qubit number is large, and thus can be used to verify a much bigger entanglement depth using the criterion in Ref. 4 (). From Eq. (10), we find that the variation . As typically , the error bar for after correction of the detection error gets significantly larger, and we need to correspondingly increase the rounds of the experiment to reduce the statistical error.

To illustrate application of the error correction method here, as an example, we apply it to detection of genuine multi-partite entanglement in graph states. For a graph state of qubits associated with a -colorable graph , the genuine -party entanglement can be detected with the following witness operator 5 ()

 WGn=3I−2⎡⎣q∑l=1⎛⎝∏k∈Ql(Sk+I)/2⎞⎠⎤⎦ (11)

where denotes the set of qubits with the th color , is the identity operator, and is the stabilizer operator for the th qubit (which is a tensor product of the Pauli operators   for the th qubit and for all it neighbors in the graph ). A state has genuine -qubit entanglement if . For an ideal graph state, all its stabilizer operators have expect values . With detection error, the value of gets significantly degraded. As an example, Fig. 1 shows the values of all for two particular -colorable graph states: a -qubit GHZ state (GHZ) and a linear cluster state (LC). We assume detection error with for each qubit. With a known magnitude , the detection error can be easily corrected by a scaling transformation , where is the support of the corresponding stabilizer operator . Fig. 1 shows that after error correction, is almost unity. Its error bar increases a bit after the correction, but is still small. To show the influence on the entanglement detection, we assume the experimentally prepared graph state corresponds to the ideal target state distorted by small depolarization noise independently acting on each qubit, so , where the noise super-operator and 6 (). In Fig. 2, we show the witness as a function of the preparation error rate , both before and after correction of the detection error (with an error rate ). For both GHZ and LC states, without correction of the detection error, we cannot detect any -qubit entanglement even for a perfectly prepared state with . After correction of the detection error, we can confirm genuine -qubit entanglement as long as the preparation error . So, correction of the detection error significantly improves the experimental performance, and the improvement gets more dramatic when the qubit number increases.

Finally, we briefly comment on the sensitivity of our error correction method to calibration of the detection error. In this method, the error magnitude (or magnitudes , for general cases) is assumed to be known. If we have a relative error in calibration of the magnitude , .i.e., , the scaling transformation on the detected operator leads to an relative error in the observed quantity . As long as , which is typically the case as , the relative error actually gets reduced and the method here can tolerate some uncertainty in calibration of the error magnitude .

In summary, we have shown a method to correct any detection error through simple processing of the experimental data. The method applies to measurements in general many-particle settings, with or without separate addressing. Different from the conventional quantum error correction, this method does not require encoding or change of the quantum circuits and works under arbitrary magnitudes of the detection noise, as long as the error magnitude has been calibrated. The cost of this method is moderate as it only requires repetition of the same experiment by some additional rounds to gain enough statistics and thus the method can readily apply to many experimental settings.

Appendix:Proof of Eq (7). We can relate the matrix to the matrix defined in Eq. (3). Denote the space of -bit binary strings with bits of as , and has dimension . The matrix element represents the probability of recording a -bit binary string as , and is the probability of recording a signal as any string in the space. As a collective measurement does not distinguish the binary strings in the same space , is related to by The probability is apparently independent of the exact form of , as long as belongs to the space , so we can pick up any in without alternation to the result of summation. From Eq. (4), we know . Let us define

 Njk≡∑μ∈SjM−1μν =∑μ∈SjMμν (p′0,p′1)=Ljk(p′0,p′1),

where is an arbitrary element in . Now we show that gives inverse of the matrix :

 ∑jLijNjk=∑j∑σ∈Si,μ∈SjMσρM−1μν = ∑σ∈Si∑j∑μ∈SjMσμM−1μν=∑σ∈Siδσν=δik

In the second line, we have changed the subscript in to as both belong to . This proves Eq. (7) in the text.

This work was supported by the NBRPC(973 Program) 2011CBA00300 (2011CBA00302), the IARPA MUSIQC program, the DARPA OLE program, the ARO and the AFOSR MURI program.

## References

• (1) M.A. Nielsen, I.L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press (2000).
• (2) H.-J. Briegel, W. Dür, J. I. Cirac1, and P. Zoller, Phys. Rev. Lett. 81, 5932 (1998).
• (3) D. Leibfried et al., Nature 438, 639 (2005); T. Monz et al., Phys. Rev. Lett. 106, 130506 (2011); R. Islam et al., Nature Commun. 2, 377 (2011).
• (4) D. Wineland, J. J. Bollinger, W. M. Itano, and D. J. Heinzen, Phys. Rev. A 50, 67 (1994).
• (5) A. Sorensen and K. Molmer. Phys. Rev. Lett, 86, 4431 (2001).
• (6) O.Gühne and G. Tóth. Physics Reports, 474, 1 (2009).
• (7) W. Dur and H. J. Briegel. Phys. Rev. Lett, 92, 180403 (2004).
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters