Two-qubit pure state tomography by five product orthonormal bases
Wang Yu1, 2, Shang Yun1, 2, 3, 4, †
Institute of Mathematics, Academy of Mathematics and Systems Sciences, Chinese Academy of Sciences, Beijing 100190, China
University of Chinese Academy of Sciences, Beijing 100049, China
National Center for Mathematics and Interdisciplinary Sciences, Chinese Academy of Sciences, Beijing 100190, China
MDIS, Academy of Mathematics and Systems Sciences, Chinese Academy of Sciences, Beijing 100190, China

 

† Corresponding author. E-mail: shangyun602@163.com

Project supported partially by the National Key Research and Development Program of China (Grant No. 2016YFB1000902), the National Natural Science Foundation of China (Grant No. 61472412), and the Program for Creative Research Group of the National Natural Science Foundation of China (Grant No. 61621003).

Abstract

In this paper, we focus on two-qubit pure state tomography. For an arbitrary unknown two-qubit pure state, separable or entangled, it has been found that the measurement probabilities of 16 projections onto the tensor products of Pauli eigenstates are enough to uniquely determine the state. Moreover, these corresponding product states are arranged into five orthonormal bases. We design five quantum circuits, which are decomposed into the common gates in universal quantum computation, to simulate the five projective measurements onto these bases. At the end of each circuit, we measure each qubit with the projective measurement {|0⟩⟨0|,|1⟩⟨1|}. Then, we consider the open problem whether three orthonormal bases are enough to distinguish all two-qubit pure states. A necessary condition is given. Suppose that there are three orthonormal bases . Denote the unitary transition matrices from to as U1 and U2. All 32 elements of matrices U1 and U2 should not be zero. If not, these three bases cannot distinguish all two-qubit pure states.

1. Introduction

A quantum state is the carrier of quantum information. One of the central topics in quantum science is the ability to characterize and manipulate quantum states, which would be applied to quantum communications, quantum cryptography, and quantum computation. Suppose that there is an unknown quantum state. It is impossible to create an identical copy of this state by the no-cloning theorem.[1] However, the information of the unknown quantum state can be transmitted to a remote location by the process of quantum teleportation,[2] with the help of classical communication and previously shared quantum entanglement. Now consider the following problem. There is a machine which can produce a large number of identical unknown quantum states. We are permitted to design any device to repeatedly measure this state and read out the probability distributions of the measurement outcomes. Can we infer all the information of the unknown state?

This is the problem of quantum state tomography. Two aspects are essential for this problem: the measurements and the method to reconstruct the unknown state via the measurement probabilities. Various measurements are considered for achieving quantum state tomography, such as Pauli measurements,[3,4] projective measurement onto mutually unbiased bases,[58] or onto general orthogonal bases,[912] and symmetric informationally complete measurements.[1315] The data of measurement probabilities is collected after we take measurements on the ensemble of identical unknown quantum states. There are different methods applied on the data in order to calculate all the coefficients of the unknown state, such as linear inversion,[16,17] maximum-likelihood estimation,[9,18,19] Bayesian mean estimation,[20] and the convex optimization program.[12,21]

In this paper, we focus on the two-qubit pure state tomography. From the result of papers,[9,22] the density matrix of an unknown two-qubit state can be determined with the expectation values of 16 Pauli measurements. It was proved that a two-qubit pure state |ϕ⟩ is uniquely determined among all states (UDA) with 11 Pauli measurements.[4] That is to say, there does not exist any other state, pure or mixed, which has the same measurement results as |ϕ⟩. The projective measurements of five mutually unbiased bases have the same property.[6] In Ref. [12], another set of five orthonormal bases are constructed by a sequence of orthogonal polynomials. The method of the convex optimization program can be used to search for a density matrix, whose measurement probabilities match the data obtained experimentally. It has been found that several projections onto the tensor products of Pauli eigenstates are enough for two-qubit pure states tomography. These tensor product states are arranged into five orthonormal bases. Furthermore, the projective measurements onto these bases is simulated by five quantum circuits. Each circuit is decomposed by the common quantum gates used in universal quantum computing,[22] which is the footstone in building a universal quantum computer and has been researched on various materials. The measurement parts are the projective measurement {|0⟩⟨0|,|1⟩⟨1|} on each qubit. We can directly calculate the unknown two-qubit pure state with partial measurement probabilities of these bases.

Moreover, a necessary condition is given for the unsolved problem of whether three orthonormal bases can distinguish all two-qubit states. This problem dates back to the famous Pauli problem, whether or not the position and the momentum distributions can uniquely determine the wave function of a particle (up to a global phase).[23] The answer to the Pauli problem is negative and many related topics have been studied. The finite version of Pauli problem is considered by Asher Peres.[24] He noted that the probabilities of two orthonormal bases should be sufficient to determine a finite dimensional pure state up to a set of ambiguities. The findings in Ref. [25] suggest that the set is not a set of measure zero. A general question is the following. How many different orthonormal bases are needed at least in order to uniquely determine a pure quantum state in d-dimensional Hilbert space ? The minimal number of orthonormal bases is gradually clear except for d = 4.[10,2629] There exist four orthonormal bases which can distinguish any pair of different pure states in , but it is not clear whether three bases can do this. There is a partial answer for this problem. It is proved that three product orthonormal bases cannot distinguish all two-qubit states, not even four. The product basis here is restricted to the form of {ϕiψj:i, j = 1, 2}, where {ϕ1, ϕ2} and {ψ1, ψ2} are two orthonormal bases of . Given any three orthonormal bases of , the unitary transition matrices from to can be written as U1 and U2. We show that all 32 elements of U1 and U2 should not be zero. If not, these three bases are not enough. For example, if we choose the first basis as the standard basis and choose the second basis as the basis of Bell states, no matter how we choose a third basis, these three bases cannot distinguish all two-qubit pure states.

The rest of this paper is organized as follows. In Section 2, we describe the five product bases and construct five quantum circuits to simulate the projective measurements onto these bases. Then, we show the method to directly calculate the unknown pure state with partial measurement probabilities. In Section 3, we give a necessary condition that three orthonormal bases must satisfy in order to distinguish all two-qubit pure states. We will end up with a conclusion in Section 4.

2. Five product bases for two-qubit pure state tomography

As a warm up, let us consider the one-qubit state tomography. It is known in many textbooks that three orthonormal bases uniquely determine an unknown qubit state,[22] which can be the eignenbases of Pauli operators σx, σy, and σz. Denote the eignenbases as {|+⟩, |−⟩}, {|L⟩, |R⟩}, and {|0⟩,|1⟩} correspondingly. Here

The Hadamard gate H can be described as | + ⟩⟨0| + | − ⟩⟨1|. Denote gate as | + ⟩⟨0| − i| − ⟩⟨1|.

Let U be a d × d unitary operation. If we want to measure a d-dimensional pure state |ϕ⟩ with the projective measurements onto the orthonormal bases {U|0⟩,..., U|d − 1⟩}, we can take the unitary operation U on the state |ϕ⟩, then measure U|ϕ⟩ with the projective measurements onto the orthonormal bases {|0⟩,..., |d − 1⟩}. As the probabilities of the j-th outcome for both sides are equal, |⟨ψ|U|j⟩|2 = |⟨j|U|ψ⟩|2, where j = 0,..., d − 1.

Let V1 and V2 be the unitary transition from basis {|0⟩, |1⟩} to basis {|+⟩, |−⟩} and {|L⟩, |R⟩}. Then V1 = H, ; V2 = |L⟩⟨0| + |R⟩⟨1|, . We can simulate the projective measurements onto these Pauli eigenbases with the circuits in Fig. 1.

Fig. 1. Three circuits for one qubit tomography. The measurement parts are the projective measurement {|0⟩⟨0|, |1⟩⟨1|}. Before measurement, we take the unitary operation H, , and I on the unknown qubit state respectively.

Now, we construct five product orthonormal bases to directly determine an unknown two-qubit pure state, which are as follows:

It is easy to check that each basis above is orthonormal. The first basis is the standard basis in 4-dimensional Hilbert space . The projective measurement onto the states of is simulated by the circuit in Fig. 2.

Fig. 2. Circuit for the projective measurement onto . Let the unknown two-qubit pure state be |ϕ⟩. We directly measure each qubit with the projective measurement {|0⟩⟨0|, |1⟩⟨1|}. The state would collapse to |00⟩, |01⟩, |10⟩, or |11⟩.

Next we use the circuits to simulate the projective measurements onto the bases , which are depicted in Fig. 3. Let Uj be the unitary transition from the basis to the basis , j = 1, 2, 3, 4. The task is to decompose the gates by the common gates in universal quantum computation. For example, . If we take the operation on the states of the basis , we will obtain the states of the basis .

Fig. 3. (color online) Circuits for the projective measurements onto bases . The circuits at top left, top right, bottom left, and bottom right simulate the projective measurements onto , and respectively. There are two two-qubit operations in these circuits. The first one is the swap operation, which is depicted in Fig. 4. The other one is a controlled unitary operation, . The one-qubit operations are the Hadamard gate H and the gate that appeared in the circuit for one qubit tomography. At the end of each circuit, we measure each qubit with the projective measurement {|0⟩⟨0|, |1⟩⟨1|}.
3. Method to calculate the unknown two-qubit pure state

An arbitrary unknown two-qubit pure state, separable or entangled, can be expressed with the vector form a|00⟩ + b|01⟩ + c|10⟩ + d|11⟩, where |a|2 + |b2| + |c|2 + |d|2 = 1. Any complex number g can be expressed as |g|eg. So we can write the unknown state with a new form

where θa, θb, θc, θd ∈ [0,2π). In this section, we demonstrate that all the unknown coefficients can be calculated by partial measurement probabilities of the five bases.

Fig. 4. Simulation of swap operation. The swap operation is: |00⟩⟨00| + |10⟩⟨ 01| + |01⟩⟨10| + |11⟩⟨11|. If we take the swap operation on input state |jk⟩, the output state will be |kj⟩, j, k = 0, 1. It can be decomposed into three control-not operations. The bases states of and would be transformed into the bases states of and by this operation.

If we measure |ϕ⟩ with the projective measurement onto the basis , the state will collapse to |00⟩, |01⟩, |10⟩, or |11⟩. Denote the measurement probabilities as , , , and respectively. The coefficients of {|a|,|b|,|c|,|d|} can be calculated by the probabilities of the basis . The following equations hold:

So, , , , . The unknown state will be clear until the phases {θa, θb, θc, θd} are calculated. We can let one of these phases be zero. Two pure states |ϕ1⟩ and |ϕ2⟩ correspond to the same one if there is a nonzero complex number g such that |ϕ1⟩ = g|ϕ2⟩, where |g| = 1. As |ϕ1⟩⟨ϕ1| = gg|ϕ2⟩⟨ϕ2| = |ϕ2⟩⟨ϕ2|, their density matrices are the same. So the global phase can be extracted and neglected.

At least one probability in is not equal to zero. Next we divide it into four cases and calculate the remaining phases {θa, θb, θc, θd} at each case. Record the probability as if the state collapses to the k-th element of basis .

For example, the probabilities of the basis are the following: , , , and . Then |a| = 1, b = c = d = 0. We do not need to calculate the phases {θb, θc, θd} and the unknown state is |00⟩ for the freedom choice of the global phase. The nonzero probability may be another one. We list the nonzero probability and the corresponding state in the following table (Table 1).

Table 1.

The nonzero probability and the corresponding state.

.

There are subcases totally. For example, . The unknown state can be written as . We only need to calculate the phase θb. This task can be fulfilled by considering the probabilities that the state collapses to |0 + ⟩ and |0L⟩. The following equations hold:

So the values of cosθb and sinθb can be determined by the probabilities . Then we assign a unique value of θb at the range of [0,2π).

When , the unknown state is |a||00⟩ + |c|eiθc|10⟩. In order to calculate the phase θc, we can consider the probabilities that the state collapses to | + 0⟩ and |L0⟩.

When , the unknown state is |a||00⟩ + |d|eiθd|11⟩. In order to calculate the phase θd, we can consider the probabilities that the state collapses to |L − ⟩ and |LR⟩.

When , the unknown state is |b||01⟩ + |c|eiθc|10⟩. In order to calculate the phase θd, we can consider the probabilities that the state collapses to | − L⟩ and |RL⟩.

When , the unknown state is |b||01⟩ + |d|eiθd|11⟩. In order to calculate the phase θd, we can consider the probabilities that the state collapses to | + 1⟩ and |L1⟩.

When , the unknown state is |c||10⟩ + |d|eiθd|11⟩. In order to calculate the phase θd, we can consider the probabilities that the state collapses to |1 + ⟩ and |1L⟩.

There are four subcases in total. For example, . The unknown state can be written as |a||00⟩ + |b|eb|01⟩ + |c|ec|10⟩. The phases θb and θc can be determined by Eqs. (5) and (6).

When , the unknown state is |a||00⟩ + |b|eiθb|01⟩ + |d|eiθd|11〉. The phases θb can be determined by Eq. (5). The phase difference θdθb can be determined by considering the probabilities that the state collapses to | + 1⟩ and |L1⟩.

When , the unknown state is |a||00⟩ + |c|eiθc|10⟩ + |d|eiθd|11⟩. The phase θc can be determined by Eq. (6). The phase difference θdθc can be determined by considering the probabilities that the state collapses to |1 + ⟩ and |1L⟩.

When , the unknown state is |b||01⟩ + |c|eiθc|10⟩ + |d|eiθd|11⟩. The phase θd can be determined by Eq. (9). The phase difference θdθc can be determined by considering the probabilities that the state collapses to |1 + ⟩ and |1L⟩.

We write the unknown two-qubit pure state as |a||00⟩ + |b|eiθb|01⟩ + |c|eiθc|10⟩ + |d|eiθd|11⟩. We can use Eqs. (5) and (6) to determine the phases θb and θc. The phase difference θdθb can be determined by Eq. (11). Then all the phases are clear.

Now we draw a conclusion. The coefficients {|a|, |b|, |c|, |d|} can be calculated by the probabilities , which are obtained from the projective measurement onto the basis . The remaining task is to calculate the phases for the nonzero amplitudes. The number of nonzero amplitudes may be k as listed in the case k, where k ∈ {1, 2, 3, 4}. We can let one phase be zero for the freedom choice of the global phase. For the cases 2, 3, and 4, we use two, four, and six measurement probabilities respectively to determine the phases. These probabilities are some parts obtained from the projective measurements onto the bases , where each basis state is the tensor product of Pauli eigenstates. Until we determine the amplitudes and the phases, the unknown two-qubit state will be clear.

4. A necessary condition

Let be m orthonormal bases of the d-dimensional Hilbert space , where . Given any pair of different pure states ρ1 and ρ2, these bases can distinguish them if for some basis states . The minimal number of orthonormal bases needed to distinguish all d-dimensional pure states is clear except for d = 4.[10] Four orthonormal bases can distinguish all two-qubit pure states; however, the problem remains unsolved whether three orthonormal bases are enough. Here we give a necessary condition that three orthonormal bases should satisfy. If not, three bases cannot distinguish all two-qubit pure states.

Denote three orthonormal bases as . We record the bases , , and as {|1⟩, |2⟩, |3⟩, |4⟩}, {|ϕ1⟩,..., |ϕ4⟩}, {|ψ1⟩,..., |ϕ4⟩} respectively. The unitary transition matrices from to are written as U1 and U2. U1 = (xi, j), U2 = (yi, j), i, j = 1,..., 4. In other words, the bases of and are the following:

All 32 elements of matrices U1 and U2 should not be zero. If zero appears, we can always construct a pair of different pure states, whose measurement probabilities under the three orthonormal bases are the same.

Here we give the proof. Without loss of generality, the zero matrix element is x11. We may as well let the first rows of matrix U1 and U2 be real numbers. Namely, xk1 and yk1 are all real, k = 1, 2, 3, 4. This is because the probabilities that state |ϕ⟩ collapses to the basis states |ϕk⟩ or c|φk⟩ are equal, |c| = 1. So the bases {|ϕ1⟩, |φ2⟩, |ϕ3⟩, |ϕ4⟩}, and {c1|ϕ1⟩, c2|ϕ2⟩, c3|ϕ3⟩, c4|ϕ4⟩} have the same power of discrimination, where |ck| = 1 for k = 1, 2, 3, 4.

We can construct a pair of unnormalized vector states:

where Re(ck) and Im(ck) are the real and imaginary parts of the complex number ck. Now we prove that there is a nontrivial solution of {ck: k = 2, 3, 4}, so that the probability distributions of |Φ1⟩ and |Φ2⟩ are the same under the projective measurements onto the bases , and .

It is easy to prove that their probability distributions onto the first basis are the same.

For the second basis ,

where , . If Xk = 0 for all k, we have |⟨Φ1|ϕk⟩|2 = |⟨Φ2|ϕk⟩|2. Their probability distributions onto the basis are the same. From the assumption, x11 = 0. We have |⟨Φ1|φ1⟩|2 = |⟨Φ2|φ1⟩|2 even X1≠ 0.

Define . With a similar analysis, if for k = 1, 2, 3, 4, their probability distributions onto the basis are the same one.

Now we consider the following linear equations:

There are five linear equations about six unknown real variables Re(cj),Im(cj): j = 2, 3, 4}. There must be a nontrivial solution. Then for such a pair of two different unnormalized vector states |Φ1⟩ and |Φ2⟩, we have |⟨Φ1|ϕk⟩|2 = |⟨Φ2|φk⟩|2, |⟨Φ1|ψk⟩|2 = |⟨Φ2|ψk⟩|2, k = 1, 2, 3. As the bases and are orthonormal, we have .

Then,

So the following equations are satisfied:

Let N be the norm of |Φ1⟩ and |Φ2⟩. The probability distributions of pure states |Φ1⟩/N and |Φ2⟩/N are the same under the projective measurements onto the bases , , and .

The zero element may appear at the j-th row of the matrices U1 or U2. Then the two vector states can be constructed as: |Φ1⟩ = |j⟩ + ∑kj [Re(ck) + i· Im(ck)]|k⟩, |Φ2⟩ = |j⟩ − ∑kj[Re(ck) + i · Im(ck)]|k⟩. With a similar analysis, there exists a nontrivial solution so that the probability distributions onto the three orthonormal bases are equal. This completes the proof.

This necessary condition can give us some insights to the minimal number of orthonormal bases which can distinguish all two-qubit pure states. If the minimal number is actually three, we should find them with the constraint of the condition. For example, if we choose basis as the standard basis and choose basis as the basis of Bell states, no matter how we choose a third basis, these three bases cannot distinguish all two-qubit pure states. As there are zero elements in matrix U1, if any three orthonormal bases which satisfy this condition cannot distinguish all two-qubit pure states, the minimal number is four.

5. Conclusion

In this paper, we show that any unknown two-qubit pure state can be directly reconstructed by the measurement probabilities of projections onto 16 tensor products of Pauli eigenstates. These projections are arranged into five orthonormal bases. Moreover, we construct five quantum circuits to simulate the projective measurements onto the states of these bases. In the end, we focus on the unsolved problem of whether three orthonormal bases can distinguish all two-qubit pure states. A necessary condition is given. If the condition is not satisfied, three orthonormal bases are not sufficient to distinguish all two-qubit pure states.

Reference
[1] Wootters W K Zurek W H 1982 Nature 299 802 https://www.nature.com/articles/299802a0
[2] Bennett C H Brassard G Crépeau C Jozsa R Peres A Wootters W K 1993 Phys. Rev. Lett. 70 1895
[3] Gross D Liu Y K Flammia S T Becker S Eisert J 2010 Phys. Rev. Lett. 105 150401
[4] Ma X Jackson T Zhou H Chen J X Lu D W Mazurek M D Fisher K A G Peng X H Kribs D Resch K J Ji Z F Zeng B Laflamme R 2016 Phys. Rev. 93 032140
[5] Yan F Yang M Cao Z L 2010 Phys. Rev. 82 044102
[6] Adamson R B A Steinberg A M 2010 Phys. Rev. Lett. 105 030406
[7] Lima G Neves L Guzmán R Gómez E S Nogueira W A T Delgado A Vargas A Saavedra C 2011 Opt. Express 19 3542 https://www.osapublishing.org/abstract.cfm?uri=oe-19-4-3542
[8] Giovannini D Romero J Leach J Dudley A Forbes A Padgett M J 2013 Phys. Rev. Lett. 110 143601
[9] Daniel F V J Paul G K William J M Andrew G W 2001 Phys. Rev. 64 052312
[10] Carmeli C Heinosaari T Schultz J Toigo A 2015 Euro. Phys. J. 69 179
[11] Goyeneche D Caňas G Etcheverry S Gómez E S Xavier G B Lima G Delgado A 2015 Phys. Rev. Lett. 115 090401
[12] Carmeli C Heinosaari T Kech M Schultz J Toigo A 2016 Europhys. Lett. 115 30001
[13] Renes J M Blume-Kohout R Scott A J Caves C M 2004 J. Math. Phys. 45 2171
[14] Řeháček J Englert B G Kaszlikowski D 2004 Phys. Rev. 70 052321
[15] Kalev A Shang J W Englert B G 2012 Phys. Rev. 85 052115
[16] Band W Park J 1970 Found. Phys. 1 133
[17] Park J Band W 1971 Found. Phys. 1 211
[18] Lvovsky A I Raymer M G 2009 Rev. Mod. Phys. 81 299
[19] Yin Q Xiang G Y Li Ch F Guo G C 2017 Chin. Phys. Lett. 34 030301
[20] Blume-Kohout R 2010 New J. Phys. 12 043034
[21] Baldwin C H Deutsch I H Kalev A 2016 Phys. Rev. 93 052105
[22] Nielsen M A Chuang I L 2000 Quantum Computation and Quantum Information Cambridge Cambridge University 172 389, 398
[23] Pauli W 1933 Die allgemeinen Prinzipen der Wellenmechanik 24 Berlin Springer-Verlag
[24] Peres A 1993 Quantum Theory: Concepts and Methods Dordrecht Kluwer Academic, The Netherlands POVMs are discussed in Sections 9-5 and 9-6, and PSI-complete measurements in Section 3-5
[25] Flammia S T Silberfarb A Caves C M 2005 Found. Phys. 35 1985
[26] Moroz B Z 1983 Int. J. Theor. Phys. 22 329
[27] Moroz B Z 1994 Theor. Math. Phys. 101 1200
[28] Heinosaari T Mazzarella L Wolf M M 2013 Commun. Math. Phys. 318 355
[29] Jaming P 2014 Appl. Comput. Harmon. A 37 413